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Abstract 

We outline an ultraviolet renormalization procedure for hamiltonians acting in the light- 
front Fock space. The hamiltonians are defined and calculated using creation and annihila¬ 
tion operators with no limitation of the space of states. Both, the regularization of the initial 
hamiltonian and the definition of the renormalized effective hamiltonians, preserve the light- 
front frame kinematical symmetries. The general equations for the effective hamiltonians are 
illustrated by second order calculations of the self-energy and two-particle interaction terms 
in Yukawa theory, QED and QCD. Infrared singularities are regulated but not renormalized. 
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1. INTRODUCTION 

This paper describes elements of a renormalization theory for hamiltonians. The two 
main elements are the similarity renormalization group for deriving the hamiltonians and 
the Schrodinger equation for their eigenstates. The hamiltonians are defined in the Fock 
space. 

The Fock space is constructed using light-front creation operators acting on the vacuum 
state. The vacuum is annihilated by the corresponding annihilation operators. We use this 
representation of states because, due to the boost symmetry, it is easier than the usual 
equal-time representation to work with in a relativistic theory. 

The similarity renormalization group is defined in terms of effective creation and anni¬ 
hilation operators which appear in the effective hamiltonians. The effective hamiltonians 
can be calculated using the method of successive approximations and perturbation theory. 
Eventually, one obtains a large number of terms and one must study their contribution to the 
effective Schrbdinger equation. The space of solutions of the Schrodinger equation includes 
bound states. 

The goal of this paper is to present the general framework. In addition, the general 
presentation is illustrated by second order perturbative calculations of effective hamiltonians 
in Yukawa theory, QED and QCD. The effective hamiltonians lead to the on-shell scattering 
matrices which agree with standard results. The bound state eigenvalue problems for the 
effective hamiltonians are only briefly described because little is known about them beyond 
the nonrelativistic limit. This article is meant to help in facilitating further studies by 
providing a comprehensive introduction and showing examples of a systematic derivation of 
the effective hamiltonians. From this point of view, details of the approach are less important 
than the general framework and they are forseen to change in future. 

One reason for renormalization to play an important role in the hamiltonian approach is 
divergences. Initial expressions for bare hamiltonians are divergent because they are obtained 
from formal expressions in local field theory. Nevertheless, such hamiltonians are known to 
closely reproduce experimental data through tree diagrams. The most prominent example 
is the canonical hamiltonian of QED. Renormalization problems appear only when one is 
forced to sum over intermediate states and the sum diverges. The divergences correspond to 
diverging loop integrals in Feynman diagrams. However, the hamiltonian approach greately 
differs from the lagrangian diagrammatic approaches. 

One apparent difference is that the hamiltonian sums over intermediate states involve 
integrals over the three-dimensional momentum space while the integrals in the lagrangian 
calculus are four-dimensional. Although a connection exists for finite integrals which are not 
sensitive to cutoffs, when the integrals diverge the connection is broken. In the hamiltonian 
approach, we introduce the effective Fock space basis and we construct the hamiltonians using 
the three-dimensional regularization and renormalization. In the lagrangian approach, one 
directly calculates Green functions using four-dimensional regularization and renormalization 
techniques. Equivalence of the two approaches in diverging cases remains to be shown, 
especially when the bound states are taken into account. 

Besides removing divergences, the renormalization group is useful in the hamiltonian 
approach because it introduces a hierarchy of scales. Phenomena of different scales are 
dealt with in a certain order. This enables us to solve problems involving many scales. 
Particle theories contain many, possibly even infinitely many different scales. Using the 
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renormalization group approach, we can start from the hamiltonian of a basic theory that 
couples all degrees of freedom of all scales and we can reduce the initial hamiltonian to a 
renormalized hamiltonian in which couplings between degrees of freedom of vastly different 
scales, vanish. Eventually, the goal is to obtain an effective hamiltonian for such a small 
range of scales that the spectrum of solutions can be found in practice. 

In the hamiltonian renormalization group approach, one introduces different scales of 
momentum. We define the momentum scales for particles using their relative momenta. For 
example, consider the case where some dressing mechanisms for particles are expected to 
work at different scales. One can think of the constituent quarks as built from the smaller 
current quarks and gluons. The relative momenta of the latter determine the resolution 
scale which is required to see them. For addressing such issues, the renormalization group 
evolution in the relative momentum scales is a natural choice. But, so far, it has not been 
exactly understood in a hamiltonian approach how the constituent structure appears in 
QCD. 

A number of model subspaces of the Fock space need to be considered when solving for 
the spectrum of a field-theoretic effective hamiltonian because the full space of states is too 
large for computations. Also, different physical problems require different model subspaces. 
Nevertheless, working within a subspace of interest, one should secure that the results for 
physical quantities are independent of the renormalization group cutoff which limits the 
relative energy transfers. The cutoff independence can appear only in a certain range of 
cutoffs which corresponds to the model subspace. However, once the cutoff independence in 
the hnite range is achieved, one expects to have solved the theory in this range. 

The following diagram will illustrate the situation. 


H{oo, n, 6, A) 

i 

i 

full i 

i 

H (A,n,(5, A) 


1 

-^ H{oo,h,6,A) 

i 

I limited 

2 i 

-- H {X,n,6,A) 


In this diagram, the vertical arrows indicate evolution in the renormalization group pa¬ 
rameter A. A limits the relative energy transfers. It will be explained in detail in next 
Sections how the limits are imposed. A ranges, in principle, from inhnity to zero. The 
hamiltonians depend on additional parameters n, S and A. This notation is used here for 
the purpose of the introductory discussion and it will be partly altered in the next Sections. 

n stands for the cutoff on the change of the particle number. It dehnes the limits on the 
numbers of creation and annihilation operators that can appear in a single hamiltonian term. 
For example, the canonical expressions for light-front hamiltonians in local held theories of 
physical interest limit the number of creation and annihilation operators in a single term 
to 4, and the particle number can change at most by n = 2. One can set similar or less 
stringent restrictions on the effective hamiltonians with hnite width A. 

S stands for the infrared cutoh. For example, it may dehne the lower bound on the 
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longitudinal momentum carried by particles that can appear or disappear as a result of a 
single interaction. 

A stands for the ultraviolet cutoff which defines the upper limit on the relative transverse 
momentum of particles which can appear or disappear in a single interaction. 

The left branch of the diagram is marked “full” because it corresponds to the renormal¬ 
ization group evolution in the full Fock space. The parameters n, S and A do not limit the 
particle number in the space of states. 

Counterterms in the initial H{oo,n,S, A) are constructed to make the physical results 
have well dehned limits when the cutoffs n, S and A are relaxed. The construction of 
counterterms in perturbation theory will be discussed in detail in the next Sections. Once 
the counterterms remove the regularization dependence from the effective dynamics the 
arguments n, 6 and A in H{X, n, S, A) in the lower left corner of the diagram are equivalent 
to their limiting values, n = cx), 5 = 0 and A = oo. Then, A) = H{\) which 

is independent of n, S and A for their extreme values. The infrared regulator S may still 
appear in the effective hamiltonian if there are massless vector particles in the theory. 

The horizontal arrow marked 1, denotes simplihcations which can be made in the initial 
hamiltonian if some of its matrix elements are small and/or one can drop certain couplings 
without losing accuracy in representing solutions of the theory. The model parameters 
h, S and A are introduced to limit dynamics to the selected part. For example, h can 
limit the number of bare particles in the initial model space. Also, 5 can limit the bare 
particle momenta from below when it is known that only momenta in a certain range matter. 
Similarly, A can limit from above the free energies of bare particle states which are taken 
into account in a calculation when it is known that some range of the free energies dominates 
the dynamics and the introduction of the upper limit does not lead to signihcant errors. 

Accuracy of the step denoted by the arrow 1 has to be checked by relaxing the model 
cutoffs and measuring the resulting changes in the spectrum of H{\,h,6, A). Naturally, 
these cutoffs may have to be varied in a big range because they are introduced along the 
arrow 1 for bare particles. 

The right branch of the hgure is marked “limited” because it describes the renormalization 
group flow in the limited bare model space selected by parameters h, S and A. 

The effective hamiltonians at the bottom of the diagram, namely, H{X,n,6, A) and 
H{X, h, S, A) with energy transfers limited by hnite A, are connected by the arrow marked 2. 
This arrow denotes the procedure of introducing helpful cutoffs which enable us to approxi¬ 
mately solve for the spectrum of the effective hamiltonian H{X, n, S, A). This time, however, 
the final computation cutoffs are introduced at the level of the effective particles, not at the 
level of the initial bare particles. 

The arrow 2 denotes the replacement of the whole effective hamiltonian matrix by a 
limited matrix which is obtained from the full matrix by introducing the model cutoffs h, 6 
and A. The procedure of obtaining the small matrix will be discussed below. The spectrum 
of the small matrix may be very close to the corresponding part of the spectrum of the full 
matrix because A is small (this will become clear later). The accuracy must be verified by 
relaxing the cutoffs h, S and A and observing convergence of results as in the case of the 
arrow 1. However, again, it is natural to expect that the cutoffs h, S and A may have to be 
varied only in a small range which corresponds to A. Thus, a hnite dynamical problem to 
solve is dehned. 

The “full” renormalization group evolution can be calculated using the approach of the 
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present paper. The “limited” evolution can be calculated using the matrix element techniques 
introduced earlier by Glazek and Wilson |]^ who drew on the work of Wilson . The 
matrix elements techniques were introduced for application to QCD. |Q Alternatively, one 
can use Wegner’s flow equations for hamiltonian matrix elements in cases soluble with the 
energy-independent width, p The present approach can be viewed as a special case of 
the earlier procedures because the hamiltonians which we consider transform by the same 
unitary transformation as a single creation or annihilation operator. However, by having 
defined the renormalization group transformation for creation and annihiltion operators, we 
remove the need to consider the model space dependence of the renormalization procedure. 

On the other hand, to make the transition to the bound state degrees of freedom, such 
as mesons and barions built from quarks and gluons, the matrix elements techniques are 
unavoidable. The transition from constituents to their bound states in the renormalization 
group flow, although achievable in principle, is not discussed in the present paper. In the 
present approach, quantum numbers of effective particles such as spin, flavor or color, cannot 
change in the flow. 

Both ways in the diagram which start from the initial hamiltonian if(cx), n, 6, A) and go 
through (a) the arrow 1 and the arrows “limited” and (b) the arrows “full” and the arrow 
2, lead to the numerically manageable H{X,n,S,A). When calculations of some matrix 
elements of the effective hamiltonian are done in perturbation theory by expanding in powers 
of a coupling constant at a selected scale A up to a limited order in the expansion, both ways 
of going through the diagram are, in principle, equivalent, because only a finite range of 
particle numbers and momenta can be reached in a limited number of steps of the size A 
starting from the finite values selected by the states used in evaluation of the matrix elements 
in question. 

Differences arise only when one solves for the spectrum of an effective hamiltonian and 
when one attempts to vary the model space parameters n, 6 and A. In the “full” calculation, 
one obtains a single effective hamiltonian which one can solve in successively enlargeable 
model spaces. In the “limited” calculation, the model space restrictions are imposed at the 
beginning and they lead to an effective hamiltonian whose action cannot be considered in 
a larger model space without repeating the renormalization group calculation in the larger 
space. 

A priori, it is not clear how important different effective interaction terms are. Pro¬ 
portionality to different powers of the coupling constant is certainly helpful in qualitative 
estimates. But it is not known precisely how to choose the basis states for evaluation of the 
effective hamiltonian matrix elements. Finding the basis which can span a good approxima¬ 
tion to the full solution will certainly require extensive trial and error studies. This feature 
can be illustrated in an elementary example of the following 2x2 matrix. 


a + bg^ gv 
gv c + dg^ 

This matrix is a model of a full effective hamiltonian calculated to second order in g including 
all couplings between all effective Fock sectors as given by the “full” calculation. Thus, we 
have the hamiltonian terms order 1, order g and order g"^. In a perturbative calculation using 
matrix elements, which is focused on the upper sectors, one would calculate only the terms 
a and bg^. 
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Assume that a and c are of the same order, and b and d are of the same order. For 
arbitrarily small g, the eigenvalues are given by a and terms quadratic in g. No terms linear 
in g arise. The quadratic corrections include contributions due to the term gv which couples 
different sectors. The role of this coupling needs to be estimated. The presence of dg^ seems 
to be irrelevant because it couples to the upper sector through the off-diagonal terms order 
g. Hence, it seems to contribute only in order g"^ to the eigenvalues. 

It is well known that the above analysis is wrong in the case with degenerate diagonal 
matrix elements. It does not matter how the degeneracy arises. For example, consider the 
case of a finite g such that a + bg'^ = c + dg"^. The eigenvalues are equal a + bg"^ ±gv. They are 
linear in g instead of being quadratic, for arbitrary v. The lowest eigenvalue eigenstate is a 
superposition of the upper and the lower sector instead of being dominated by the upper one. 
In this example, the degeneracy is not visible until the term dg^ is included in the calculation. 
As a second example consider the case with degenerate matrix elements c = a -|- bg^ and 
d not included. The simple nondegenerate perturbative expansion is again not applicable. 
But the addition of the term dg^ can lift the degeneracy and make the simple perturbative 
expansion work. 

Corrections originating from the interaction terms such as gv may be suppressed in the 
power expansion in the coupling constant for small A despite the degeneracy. The argument 
is following. The range of v in momentum variables is given by A. If A is reduced in the 
renormalization group flow down to a number which implies restriction on the energy trans¬ 
fers on the order of some positive power of g then, the resulting interaction can contribute 
to the eigenvalues only in the order implied by g and A together which is higher than g^. In 
addition, the effective interaction v may contain small factors. For example, in the effective 
e’''e“ sector of positronium, the emission of photons is proportional to the velocity of elec¬ 
trons which is order a, on average. The interaction term gv which couples states with an 
additional photon, plays no role in the eigenvalue in order if the width A restricts energy 
changes to order rueiectron (c.f. 0). 

Terms such as dg'^ have been discussed in the light-front approach to QCD. [§] In heavy 
quarkonia which are dominated by the QQ state, terms such as dg"^ in other Fock sectors than 
the QQ, may be expected to lift up energies of effective gluons to a sufficiently high value 
so that the effective hamiltonian part a + bg‘^ in the QQ sector alone may have eigenstates 
which are a good approximation of the full solution for the heavy mesons. 

But it is not known yet what happens in the eigenvalue equations of hamiltonians cal¬ 
culated to higher order than second. Unfortunately, the space of states with fast moving 
fermions and bosons is highly degenerate. A large number of matrix elements must be 
calculated and examined. One needs to choose the states required for evaluation of the 
hamiltonian matrix elements before it is known which of those states or matrix elements are 
important. Since the states are built using creation and annihilation operators, the operator 
calculus for hamiltonians provides a relatively small number of operators which universally 
produce the very many matrix elements in question, such as gv, bg^ and dg"^ and similar 
higher order terms. In particular, the operator calculus becomes essential when the number 
of states involved in the renormalization group flow is large. 

In general, when one attempts to perform a multistep renormalization group evolution 
which is indicated by the small vertical arrows, at each step using an expansion in powers 
of the corresponding coupling constants, the degree of dressing of bare particles grows with 
the number of steps. Eventually, the number of steps can approach infinity when the ratio 


6 


of a single step in the cutoff reduction to the initial cutoff tends to zero. In this case, it is 
not clear that both branches (a) and (b) in the diagram above lead to the same hamiltonian 
H{X, n, S, A) for effective particles. 

The reason is that the number of bare particles in the initial model space may have to be 
enormously enlarged to support the effective particle picture in the small hnal model space. 
The same concerns the momentum cutoffs. The required enlargement of the initial space of 
states may grow with the number of steps in the renormalization group flow and even become 
inhnite in the vicinity of a hxed point. Then, proceeding along the branch (a) by introducing 
the model limitations h, S and A at the beginning of the renormalization group evolution 
(arrow 1) and working with matrix elements of the hamiltonians in the limited space (arrows 
“limited”), may give different results from those that come out from the calculation along 
the branch (b). In the latter one, the model space cutoffs h, S and A are introduced only in 
the hnal small energy transfer hamiltonian (arrow 2) when solving for its spectrum, with a 
numerically verihable accuracy. 

An explicit example of a difference between proceeding along the two branches of the 
diagram above is provided in Ref. 0 which discusses an approach analogous to the branch 1 
and “limited” (c.f. Ref. |]^). In the Tamm-Dancoff (TD) approach, there are restrictions on 
the particle number which naturally lead to the sector-dependent counterterms as described 
in Ref. P], for example, for masses. On the other hand, in the procedure of the branch 
“full” and 2 no such sector dependent counterterms arise. The present paper will provide 
examples of sector-independent mass counterterms. 

The above diagram and the 2x2 matrix model illustrate the general structure of our 
light-front hamiltonian approach to quantum held theory. We summarize the steps here. 

The hrst step is the calculation of the effective hamiltonian. 


H{X) = Sl^^sA H{oo,n,6,A) Sx,n,SA • (1-1) 

S denotes the similarity transformation. Eq. (1.1) corresponds to the arrows marked full 
in the diagram. This equation represents the ideal situation for massive particles where the 
ultraviolet renormalization removes dependence of the effective hamiltonian matrix elements 
on n, 6 and A when the regularization is being removed. The similarity renormalization 
scheme is presented in the next Sections. 

The second step is the solution of the effective Schrodinger equation 


H{X) 1^) = E IV^) . (1.2) 

H{X) has the same dynamical content and eigenvalues as H{oo, n, 6, A) and E is independent 
of A. But Eq. (1.2) greately differs from the eigenvalue equation for H (cxo, n, 5, A). The major 
difference is that the dynamics of H{X) has a limited range on the energy scale. Therefore, 
one can hope it is possible to solve the eigenvalue problem scale by scale. Scattering processes 
are described by the same hamiltonian. Next Sections will give 2nd order examples. 

Solutions to Eq. (1.2) provide renormalization conditions for the hnite parts of the 
counterterms. When one has to go beyond perturbation theory in the solution, as it is the 
case in the bound state dynamics, a general method is necessary for reducing the whole 
problem to a manageable one. This step is marked by the arrow 2 in the diagram. In 
the case of the matrix model, this step corresponds to the calculation of the model space 


7 



hamiltonian in the upper-left corner of the matrix. The similarity renormalization scheme 
guarantees that this step is free from ultraviolet divergences because the width A is hnite. 


In the general case, one can apply the well known Bloch technique of calculating 
the model space hamiltonians (see Ref. ^ for the application of this technique in the 
usual renormalization group approach). Suppose we want to evaluate the model two-body 
hamiltonian using H{X) with A < m, where m is the one-body mass. We can introduce 
the projection operator P on the effective two-particle sector with a limited center-of-mass 
energy. We also introduce the operator R which generates the multi-particle and high energy 
components of eigenstates from their limited mass two-body part. By assumption, R satishes 
the conditions (1 — P)R = RP = R and PR = R{1 — P) = 0 and the equation P + R — 
1) H{X) {P + R) = 0. Then, the model two-body dynamics is described by the hamiltonian 

H2 = (P + R^R)-X^ (P + P^) H{X) (P + P) (P + P+P)-i/2 _ 

The same approach can be used for larger model spaces. The model space is characterized by 
the parameters h, S and A in the diagram. So, the operation P depends on these parameters. 
But the resulting spectrum in the range of interest should not depend on the model space 
boundary once it encloses the dominant basis states in the selected range of scales. 

The heuristic Eq. (1.3) can be applied in perturbation theory even for sizable coupling 
constants since the interaction strength is considerably reduced by the similarity factors. 
These factors limit the interaction range in the relative momentum space of the interacting 
particles to be of the order of A. Thus, the perturbative expansion relies on the form factor 
suppression of the interaction vertices which change the particle number. In principle, no 
A-dependence is generated in the model dynamics because the form factors are generated 
through the transformation S. This feature will be illustarted in the lowest order examples 
in Section 3. Note, that the operation P is ultraviolet hnite. 

The general scheme outlined above is still prone to the infrared regularization depen¬ 
dence which may appear in Eq. (1.3), in particular, in gauge theories. Extensive studies are 
required to hnd out if the model dynamics for gauge invariant quantities is infrared conver¬ 
gent. One may discover new efective interactions coming from the infrared region. 

m 
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Construction of a general hamiltonian approach to particle dynamics by necessity touches 
upon fundamental issues. Therefore, these issues are discussed in this paper as they ap¬ 
pear. Some issues, however, especially those related to the ground state formation, are not 
discussed extensively because such discussion would have to include description of states 
containing very large numbers of bare particles. Dynamics of such states is not understood. 
For example, it is not known if going through the branch “full” of the diagram above one 
can obtain interactions which couple few-particle states with practically inhnitely many wee- 
particle states. Similarly, issues concerning gauge invariance and rotational symmetry are 
not extensively dicussed because they require better understanding. For example, current 
conservation leads to important cancelations in the on-energy-shell hamiltonian matrix el¬ 
ements. Such cancelations are absent in the off-energy-shell hamiltonian matrix elements. 
The on-energy-shell hamiltonian matrix elements exhibit rotational invariance which is ab¬ 
sent in the off-energy-shell hamiltonian matrix elements. 

Granted all the reservations, the calculus we develop in this paper is, in a way, similar 
to the transformation discussed by Melosh [T^. The important differences are that we 


provide a renormalized dynamical theory applicable to particles of different kinds. If one 
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restricts attention to QCD, the calculus is expected to help in establishing a connection 
between the constituent quark model, Feynman parton model, and perturbative quantum 
chromodynamics. 

This paper is organized as follows. Section 2 describes the general formalism. Section 3 
contains examples of lowest order applications. Section 4 concludes the paper. The list of 
references is focused on the similarity renormalization group approach to light-front hamil- 
tonian dynamics in the Fock space. The reader should be aware of this limtation. 

2. EFFECTIVE HAMILTONIANS 

This Section is divided into three parts. The hrst part describes a general method of 
calculating effective hamiltonians in the Fock space. The second part presents our regular¬ 
ization scheme for bare hamiltonians. The last part discusses renormalization conditions and 
the effective eigenvalue problem. 


2.a Similarity transformation 

We construct a family of effective hamiltonians in the light-front Fock space. The family 
is parametrized by a scale parameter A which ranges from inhnity to a hnite value. A limits 
energy transfers in the interaction terms. 

The hamiltonians are written in terms of sums of ordered products of creation and an¬ 
nihilation operators. The hamiltonian labeled by A is expressed in terms of creation and 
annihilation operators which correspond to A. We commonly denote these operators by 
qx- In addition, the creation and annihilation operators carry labels of quantum numbers 
such as momentum, spin, flavor or color. We will not indicate those numbers in the initial 
presentation, unless it is necessary. 

All hamiltonians in the family are assumed to be equal. Thus, 

Hx,{q\i) = Hx^{qx^). ( 2 . 1 ) 

For A = cx), the hamiltonians ifoo are expressed in terms of operators creating and 
annihilating bare particles, goo- Hamiltonians H^o can be constructed from the canonical 
held theoretic expressions for the energy-momentum density tensors. 

Unfortunately, expressions for in local held theories are divergent. They need to 
be regularized by introducing a bare ultraviolet cutoh which we shall denote by e. The 
ultraviolet cutoh A from the previous Section corresponds to A^/e where A is an arbitrary 
hnite constant which carries the necessary dimension of a mass. The limit of removing the 
bare ultraviolet cutoh will correspond to e —0. 

Hx=oo = for all values of e. For the limit e —0 to exist the hamiltonians 
must include a number of additional terms (called counterterms) whose structure will be 
determined later. 

Hoo may have to include an infrared regulator, generically denoted by 5. For example, 
this is required in QED with massless photons and in QCD with massless gluons. 5 —0 
when the infrared regularization is removed. The parameter 5 is indicated explicitly only if 
needed. 

Our key assumption is that the particle degrees of freedom for all diherent scales A are 
unitarily equivalent to the bare particle degrees of freedom: 
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q\ = UxqooUl. ( 2 . 2 ) 

This assumption says that the quantum numbers of bare and effective particles are the same 
for all values of A. A few examples will illustrate the physical origin of this assumption. 
(1) Constituent quarks have the same quantum numbers as current quarks. (2) We use the 
same quantum numbers for photons and electrons independently of the kind of processes 
we consider in QED or in related effective theories such as the nonrelativistic Schrodinger 
equation with Coulomb potentials between charges. (3) Pions and nucleons in nuclear physics 
have the same quantum numbers quite independently of what kind of models, pion-nucleon 
vertex form factors or dynamical assumptions one uses. 

It follows from Eq. (2.2) that creation and annihilation operators for Ai and A 2 are uni- 
tarily equivalent and connected by transformations of the form The transformations 

Ux^ or Ux 2 will depend on the bare cutoffs but the transformation Ux^Ul^ for hnite Ai and 
A 2 will have to be hnite in the limit e —>■ 0. On the other hand, these transformations may 
depend on the infrared cutoff in the case of massless particles, unless there is a mechanism 
removing this dependence. 

We shall dehne transformations Ux using equations of the similarity renormalization 
scheme of Glazek and Wilson i i- That scheme was originally developed for application 
to QCD §. 

The transformation Ux is dehned indirectly through a differential equation for effective 
hamiltonians. Hamiltonians with labels Ai and A 2 are connected by integration of the dif¬ 
ferential equation from Ai to A 2 . Our guiding principle in writing the differential equation 
is that the transitions between states of effective particles with considerably different relative 
momenta should be suppressed. 

A few examples of effective theories illustrate why we adopt this principle. (1) Emission 
or absorption of short wavelength photons are not essential in the formation of atoms. (2) 
Emission of hard pions from a nucleon is not important in nonrelativistic nuclear physics. 
(3) Constituent quarks have moderate momenta and their effective dynamics seems to be 
independent of the very hard gluon emissions. Similarly, (4) the high momentum transfer 
phenomena are independent of the small momentum transfer effects such as binding. A 
standard way of achieving this kind of picture in theoretical models is to include form factors 
in the interaction vertices. The form factors quickly tend to zero when momenta change by 
more then the size of a specihc cutoff parameter. 

An alternative mechanism for decoupling of different momentum scales is the quick rise 
of the free energy as a function of the momentum and, effectively, appearance of the energy 
gap between different momentum scales. For example, this second mechanism works in 
atomic physics where the relativistic expression for a free electron energy is replaced by the 
nonrelativistic one. The linear dependence of the energy on a large momentum is replaced 
by a quadratic one. The quadratic rise of the energy suppresses contributions of high energy 
intermediate states. This suppression stabilizes the nonrelativistic bound state theory with 
a small coupling constant in the small relative momentum range. In a relativistic case, we 
have to include form factors to suppress transitions between different momentum scales. 

The cutoff parameter in the form factors sets the scale for allowed changes of momenta. 
It determines the range or width of the interaction in momentum space. Precisely this 
notion of the interaction width in momentum space is the origin of our scale A which labels 
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renormalized effective hamiltonians. The effective hamiltonians will contain factors that are 
analogous to the vertex form factors which are commonly used in nonlocal models (see also 
Ref. PI). 

The transformation U\ transforms the regularized hamiltonian of a local theory 
into a nonlocal effective hamiltonian Hx in which the large momentum transfer dynamics is 
integrated out. 

The momenta of individual particles are not restricted. This is required to obtain a boost 
invariant spectrum of solutions. The reason is that boosts change the momenta unlimitedly. 
Restricting the individual particle momenta would exclude boost invariance beyond the width 
scale. This would be a problem because the width needs to be small for the hamiltonian 
eigenvalue equation to be soluble in practice. 

Also, the larger is a relative momentum the larger change is generated by a boost. There¬ 
fore, when the free energy of interacting particles in their center-of-mass frame (i.e. the free 
light-front invariant mass) is larger than A then, the immediate change of energy due to the 
interactions will be limited by the large energy itself instead of A. This construction also 
reflects the property of quantum systems that strong interference occurs between waves of 
similar wavelengths within a range of wavelengths on the order of the wavelengths them¬ 
selves. 0 

The infinitesimal transformations are constructed in such a way that small energy de¬ 
nominators cannot appear in perturbative calculations of the effective hamiltonians. 0 0 
Only large energy changes are integrated out. Therefore, the dynamical interference effects 
for states of similar free energies are not included in the derivation of an effective hamilto¬ 
nian. The calculation of the strong coherence effects is deferred to the later step of solving 
for the eigenstates of the effective hamiltonian. This second step may be nonperturbative. 
For example, the Coulomb potential of QED is formally of the hrst order in a but it leads 
to an overwhelmingly rich atomic structures beyond perturbation theory. 

Our differential equations will require a separation of the changes in the creation and 
annihilation operators from the changes in the coefficients in front of their products. It 
is also initially assumed that terms with large numbers of operators in a product will not 
dominate or mediate the dynamics of interest. If the latter assumption turns out to be invalid 
then, the present formalism may only provide a way to approach the resulting problems, but 
that case is not in the focus of the present paper. The only comment due here is that, if the 
ground state dynamics leads to spontaneous symmetry breaking, or condensates, in a well 
dehned renormalized hamiltonian theory then, we will have a tool to study those phenomena 

[i> 


theoretically in great detail, c.f. Refs, p, |jT2| and |^. Nothing more can be said about it 
at this point. All hamiltonians we consider can, by assumption, be very well approximated 
by a hnite sum of hnite products of creation and annihilation operators. 

The unitary equivalence of the creation and annihilation operators for the scale A and the 
creation and annihilation operators at the inhnite scale, i.e. those appearing in H^o = 
and the equality of the hamiltonians at all scales, imply that 


Hx{qx) = UxHx{qoo)Ul = H^{q^). 
We denote Hx^q^o) ='Hx and obtain 


(2.3) 


Wa = UlH^U>.. 


(2.4) 
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Thus, the effective hamiltonian Hx for effective particles at the scale A is obtained from the 
hamiltonian Tix by replacing the creation and annihilation operators for the bare particles 
by the creation and annihilation operators for the effective particles with the same quantum 
numbers. The bare creation and annihilation operators are independent of A. We only need 
to calculate the A-dependent coefficients in front of the products of q^a in Hx. 

The differential equation for Hx is proposed in the following form 

= (2.5) 

where 

Tx = ( 2 . 6 ) 

Instead of calculating Ux we shall be using Tx in our calculations. With the help of Eq. 
(2.5), all calculations can be carried out using the bare creation and annihilation operators. 
Subsequently, at the end of the calculations, one replaces the bare operators by the dressed, 
effective ones and one obtains the desired effective hamiltonian Hx{qx) for the scale A. 

Description of the key elements of the calculation of Hx starts here. Hx has the following 
structure 


Hx = Fx[Gx]. (2.7) 

Fx[Q>] denotes an operation on the operator Qx which inserts special numerical factors and 
makes Hx have the properties we want Hx to have to represent dynamics of effective particles 
at the scale A after q^o in Qx is replaced by qx- Using the unitary equivalence, we also have 

Hx{qx) = Fx[Gx{qx)i ( 2 . 8 ) 

where 


Gxiqx) = UxQxUl 


(2.9) 


A similar relation holds for Tx{qx) and Tx since the latter is expressed in terms of goo- 

The operation Fx will be dehned now in a number of steps. We begin with a dehnition 
of operators which Fx acts on. The operator Qx, or any other operator we will consider, is 
a superposition of terms each of which is an ordered product of creation and annihilation 
operators. The ordering is arbitrary but needs to be dehned and the sign in front of each 
term needs to be determined. For example, we will adopt this order: creators of fermions, 
creators of bosons, creators of antifermions, annihilators of antifermions, annihilators of 
bosons, annihilators of fermions. So, each term contains a product of creation operators 
standing to the left of a product of annihilation operators. At least two operators must 
appear in a product and at least one creation and one annihilation operators must appear, 
i.e. no product contains only creation or only annihilation operators, by dehnition. This is a 
universal property of light-front hamiltonians. Other forms of hamiltonian dynamics do not 
have this property and immediately lead to the problem of understanding the ground state 
formation before any of its excitations can be considered beyond perturbation theory. We 
shall not further elaborate on this issue here (see e.g. [Q and ||13|| ). 
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The operator Q\ is divided into two parts, Qi\ and Q 2 \. Qi\ is the part of Q\ which 
is a superposition of the terms of the form a^a for a equal q^o of any kind. In principle, 
one could also include in Qix some terms with a larger number of creation and annihilation 
operators, e.g. terms containing two creation and two annihilation operators. This option 
may be useful in cases where the effective particles interact strongly and the formation of 
resonances or bound states needs to be taken into account already when one evaluates the 
effective hamiltonians. However, in general, the plane-wave Fock space basis states are not 
eigenstates of such operators and one is forced to consider more complicated basis states. 
Such cases are beyond the scope of the present article. Here, we limit Qw to terms a^a. 

Qi\ becomes the effective free part of Ga, denoted Gia, after q^o is replaced by q\. The 
effective free hamiltonian part will be equal to Gia because Gia will not be changed by the 
operation F. Eigenvalues of Gia will be also called free energies. 

The remaining part of Qx, he. Q 2 X = G\ — Gix, contains all terms besides terms of the 
form a^a. To obtain the corresponding part of the effective hamiltonian Hx, one first replaces 
qoo by qx and obtains G 2 a- Then, one applies the operation Fx to G 2 a- Ta inserts the vertex 
form factors. 

Let the momentum labels of all creation operators in a single product in an interaction 
term be ki, /c 2 ,..., kj and the momentum labels of all annihilation operators be k[, /C 2 ,..., kf. 
Each momentum has three components, k~^ ranging from 0 to cx) and two transverse com¬ 
ponents k-^ = {k^, k"^), both ranging from —00 two -|-cx). The 2 :-axis is distinguished by our 
choice of the light-front. The three light-front momenta are conserved. The sum of mo¬ 
mentum labels of all creation operators, ki, equals the sum of momentum labels of all 
annihilation operators in this term, I]/=i kj. We denote these sums by P = (P’*', P"*"), (P’*’ 
is positive). Each momentum label has three components, p = Each k^ or k'^ is 

a positive fraction of P"*"; Xi = kf /P"*", 1 > > 0 and x) = A;'^/P+, 1 > x) >0. We have 

Xi = J2j=i x'j = 1. We also define 

= kn - XnP^ ( 2 - 10 ) 

for all momenta in the hamiltonian term. = Z]/=i = 0. 

The above dehnitions of momentum variables are standard for wave functions in the 
light-front Fock space. But the way they are used here is not standard. We introduce this 
notation for all translationally invariant operators. The dehnition of P is not connected with 
any particular state. P is not changed as a whole by the action of the hamiltonian term. 
It is only redistributed from the set of momenta of the annihilated particles to the set of 
momenta of the created particles. 

Thus, each term in the hamiltonian is characterized by P and two sets of variables, 
X/ = {{xi, for creation operators and X'j = {{x'p K'^)yjZ'[ for annihilation operators, 

in that term. For example, if a product contains two creation operators and one annihilation 
operator, we have xi = x, X 2 = 1 — x and x[ = 1. Also, nf = —nf = and = 0. But 
P can be arbitrary and only a part of the total momentum of some state in which action of 
the product in question replaces one particle of momentum P by two particles of momenta 
xP + K and {1 — x)P — K (the latter expressions are valid for -|- and T components, = 0). 
It is convenient to speak of P as a parent momentum and about the individual particle 
momenta as daughter momenta. The parent momentum in a hamiltonian term equals one 
half of the sum of momenta labeling all creation and annihilation operators in the term. 
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Each daughter particle carries a fraction of the parent momentum. The parent momentum 
may be carried by one or more particles. 

The operation on a product of creation and annihilation operators produces 


F, 


i=l 


J 

ki n ^k'. 




= h{X,,X'j)X{alX{au' 


( 2 . 11 ) 


Z=1 


j = l 


The function f\{Xi,Xj) can be any suitable function of the daughter momentum vari¬ 
ables which represents our physical intuition and suits our calculations. The arguments of 
fx are explicitly invariant with respect to seven kinematical Poincare transformations of the 
light-front frame. This feature leads to explicit symmetry of our theory with respect to 
these transformations. There are three additional conditions that the function fx will have 
to satisfy. 

The hrst condition is that fx should be expressible through the eigenvalues of Gix cor¬ 
responding to the sets Xi and X'j so that fx equals 1 for small differences between the 
eigenvaules and it quickly goes to zero when the differences become large. This is an ap¬ 
plication of the basic condition of the similarity renormalization scheme for hamiltonians. [Q] 
1^ The width of fx is set by A. One can consider functions fx which depend on Xj and X'j 
in a more general way than through the eigenvalues of the free hamiltonian but that option 
will not be discussed here. 

The hrst condition dehnes the effective nature of the hamiltonian labeled by A. Namely, 
the effective particle states which are separated by the free energy gap which is much larger 
than A, are not directly coupled by the interactions. In other words, A limits the free energy 
changes induced by the effective interaction. Moreover, as a consequence of /a ~ 1 for similar 
energies, 1 — fx is close to zero for the similar energies and it vanishes proportionally to a 
power of the energy difference. The higher is the power the smaller is the role of states of 
similar energies in the calculation of the effective hamiltonian. This will be explained later. 
Consequently, the higher is the power the smaller is the role of nonperturbative phenomena 
due to energy changes below the scale A, in the calculation of the effective hamiltonian. 
Thus, there exists a chance for that the full hamiltonian diagonalization process can be 
divided into two parts: a perturbative calculation of the effective renormalized hamiltonian 
and a nonperturbative diagonalization of that effective hamiltonian. This is our factorization 
hypothesis in the hamiltonian approach. 

The second condition is that both, 1 — fx and dfx/dX, must vanish faster than linearly in 
the free energy difference. This condition is required to exclude small energy denominators 
in perturbation theory. This will become clear below. The second condition implies that 
1 — fx vanishes as at least second power of the energy difference near zero. 

The third condition for fx is most difficult to satisfy and a trial and error approach will 
be required to verify if this condition is satisfied. The third condition is defined by saying 
that multi-particle interactions (especially interactions that change the number of effective 
particles by many) should not be important in the effective hamiltonian dynamics which 
is characterized by changes of energies below the scale A. This may be possible if fx as a 
function of the daughter variables approximates the shape of one particle irreducible vertices 
which is characteristic to the theory under consideration. Structure of Gx depends on the 
choice of the function fx- Some choices will lead to more complicated effective interactions 
than others. The best choices for the most efficient description of physical phenomena at 
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some scale A are such that the effective particles interact in a way that is most easy to 
understand and which can be parametrized with the least possible number of parameters 
over the range of scales of physical interest. Therefore, one can conceive variational estimates 
for the best choices of fx that minimize complexity of the effective hamiltonians. Nothing 
more can be said in this work about the accuracy one can achieve in satisfying the third 
condition in general. But it is clear that creation of effective particles will be suppressed 
when the width A becomes comparable to the effective masses of those particles. 

To satisfy the first condition above we define a boost invariant gap between free energy 
eigenvalues for effective particles which is to be compared with the running cutoff parameter 
A. The reason to avoid dependence of fx on the parent momentum P is that otherwise the 
effective dynamics could depend on the frame of reference, even within the class of frames 
reachable by kinematical symmetries of light-front dynamics. The difficulties may emerge 
despite the fact that no physical result depends on the effective cutoff scale A, by construction, 
and no dependence on the frame of reference can emerge through dependence of fx on A. 
However, when approximations are made, the symmetry is no longer exact. The cutoff 
independence may appear only in a window corresponding to the approximations made. 
The problem with boost invariance is that boosts can change momentum arbitrarily and in 
order to obtain fully boost invariant results the window would have to be infinite. Otherwise, 
one would obtain frame-dependent results for, a priori, frame-independent quantities. 

If the explicit boost invariance is not preserved by the procedure of calculating the ef¬ 
fective hamiltonian then, the boost invariance condition for bound state masses provides a 
good check on the theory. But no clue is provided within the perturbative calculation of the 
effective hamiltonian for what to do if the nonperturbative diagonalization of the effective 
hamiltonian does not produce boost invariant eigenmasses. Therefore, we choose to preserve 
kinematical symmetries of light-front dynamics explicitly. The price we have to pay for this 
great simplification is the length of definition of our procedure. 

As a side remark, we should mention that rotational symmetry of the spectrum, if ob¬ 
tained, depends on interactions. Our procedure does not explicitly preserve rotational sym¬ 
metry. Nevertheless, it has been shown that if counterterms provide enough freedom through 
their finite parts and multi-particle effects are suppressed then, there exist reasons to believe 
that the rotationally invariant results may be obtained. [0 

The free energy eigenvalues relevant to a particular hamiltonian term with daughter 
variables Xj and Xj are 


and 


where 


and 


I 

E 

2 = 1 


kt 


P^2 ^ 

P+ 


^ 3 


P+ 


M] = Y. 

2=1 


+ m?(A) 


Xj 


( 2 . 12 ) 


(2.13) 


(2.14) 
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(2.15) 


M 


l2 


J 

= E 

i=i 


+ m^(A) 


x'- 


The individual effective particle masses are allowed to depend on the effective hamiltonian 
width parameter A. We dehne the mass difference for a hamiltonian term to be 


AM^ = M'j - Ml 

and the mass sum to be 


(2.16) 


j:M‘^ = M'l + Mj. (2.17) 

Now we will describe details of the function f\{Xj,Xj). We introduce a parameter za- 
Different ways of defining this parameter are useful for different purposes. It is important to 
keep in mind the physical requirement of the renormalization group approach and quantum 
mechanics that large free energy states will be allowed to interact with considerably larger 
free energy changes than small free energy states are. Consequently, the mass differences 
for masses smaller than or on the order of A will be limited by A and the mass differences 
for masses much larger than A will be limited by the size of the large masses. Following the 
similarity renormalization scheme [Q |Q, za can be chosen in such a way that za is close to 
zero for AM^ small in comparison to A^ or T^M^ and \zx\ is close to 1 for AM^ large in 
comparison to A^ or comparable to HM^. For example, 


AM^ 

EM^ + X^' 


( 2 . 18 ) 


This type of definition (including EWl^) is useful for estimates that are necessary in high 
order perturbation theory, especially in the analysis of overlapping divergences. |I[] The new 
feature here is that the introduction of YiM^ does not violate the light-front boost invariance 
and basic cluster decomposition properties. Other options exist too. For example, one can 
consider functions of the ratio AMf^/YM'^. Such functions would depend on A as a separate 
parameter so that the ratio would have to be closer to 0 for smaller A’s. There exists also 
an option of replacing Eq. (2.18) by a more complicated function of AM^, YM^ and A 
which could accelerate or slow down the decrease of the effective hamiltonian width with the 
decrease of A or the masses themselves. In principle, one can also use variables x and k in 
other combinations than M^. 

The function fx{Xi, Xj) is defined for the purpose of this article to be a function of zf", 
n > 1, which is analytic in the vicinity of the interval [0,1] on the real axis, equals 1 for 
zx = 0 and quickly approaches 0 for za ~ 1; 


fx{X,,X'j)=f{zT). (2.19) 

For example. 


/(«) 


/ u{i - Mo) y 
Vmo(1 -u) j 


( 2 . 20 ) 


where 1 > mq > 0 and m > 1. The larger the exponent m the closer f{u) approaches 
9{uo — u) for 0 < M < 1. Eq. (2.20) concludes our definition of the operation Fx. 
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The smallest possible value of in Eq. (2.18) is + J2j=i^j{^) ■ 

Thus, ;2 a is small for small positive when AAi^ is small in comparison to particle masses. 
Therefore, uq must be much smaller than 1 to force AAi^ to be small in comparison to 
when A^ is small. 

One can also force AM.^ to be small in comparison to the particle masses by making 
A^ negative so that it subtracts from its minimal value. Then, the mass difference is 

compared to the snm of kinetic energies due to the relative motion only. 

It is also usefnl to limit the small mass differences by choosing an inhnitesimally small 
Uq and introducing A^ = Uq A^. Then, |AA^^| < A^ in the 6'-fnnction limit. In this case, 
the band-diagonal hamiltonian width becomes independent of the mass sum for as long as 
the latter is small in comparison to A^. 

Next, we dehne the inhnitesimal transformation T\ in Eq. (2.5). Eq. (2.5) can be 
rewritten using Eqs. (2.8) and (2.9), into the form 

w = fg + fQ' = [/^i, T] + [/^ 2 , r \. ( 2 . 21 ) 

The prime denotes differentiation with respect to A. We have simplihed the notation of 
FxiQx] to fg. 

In the further development, three nniversally valid relations will be often used without 
saying. Namely, /^i = and (1 - /)^i = fVi = 0. 

Equation (2.21) involves two nnknowns, g and T. Additional conditions are reqnired to 
dehne T. One recalls that if the interaction is absent, i.e. when ^2 = 0, then no evolution 
with A may appear. Therefore, in the limit of negligible interactions, both g' = 0 and 
T = 0 . Thus, we expect that g' should differ from zero if and only if the interactions are 
important. The hrst term on the right-hand side is order T since is order zero in powers 
of the interactions. The second term on the right-hand side is at least of second order in the 
interactions. The hrst term can be used as a seed for dehning T through a series of powers 
of the interaction. 

We associate the derivative of g with the second term on the right-hand side. The hrst 
term on the right-hand side and a part of the second term which is left after the derivative 
of g is dehned, together determine T . Thus, we come to the dehnition of T through the 
commutator [^i,T]. We write 

A = {B}g, (2.22) 

when 

[A,g,]=B. (2.23) 

The subscript of the curly bracket in Eq. (2.22) will be often omitted in later discussion 
nnless this omission might lead to a confusion with the usnal cnrly brackets. The bracket 
operation of Eqs. (2.22) and (2.23) can be written in the explicit form since the most 
general structure of the operators A and B we shall enconnter is a snperposition of prodncts 
of creation and annihilation operators. Snppose B contains a term which involves a product 

Xialilay (2.24) 

i=l j = l 
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Then, contains the same product (as a part of the same expression) with an additional 

factor equal 


^ k'f + m]{\) 

3=1 ^ 3 


i=l 


kt 


-1 


For translationally invariant operators, the sums of individual energies satisfy Eqs. 
and (2.13), respectively, and the factor (2.25) equals 


(2.25) 

( 2 . 12 ) 


AM 2 

P+ 


(2.26) 


where is the parent momentum for the product under consideration and the mass dif¬ 
ference is dehned in Eq. (2.16). When all terms in the operator B are multiplied by the 
corresponding factors (or, equivalently, divided by the energy denominators) one obtains the 
operator A defined in Eq. (2.22). 

The factor (2.25) explodes to inhnity when the denominator approaches zero. Hence, 
for the operator A to be well dehned, the coefficients of products of the form (2.24) in 
the operator B must vanish at least as fast as the energy denominator itself when the 
denominator approaches zero. Therefore, our dehnition of [Qi,T] must be given in terms of 
an operator which has such property. In order to satisfy this condition, Eq. (2.21) is split 
into two equations as follows. 


/^' = /[/^2,r]. 


(2.27) 


[T,g3\ = {l-f)[fQ2,T]-fG. (2.28) 

These equations imply that T does not contain small energy denominators provided the 
functions f\ and /{ satisfy our initial assumptions. Our second condition introduced below 
Eq. (2.11) on the functions 1 — Xj) and fx{Xi, X'j) guarantees that T is well dehned 

and tends to zero in the region of vanishing energy denominators because the right-hand 
side of Eq. (2.28) vanishes at least as fast as the hrst power of the energy diherences. 

Equation (2.27) is a hrst order diherential equation. One has to provide initial conditions 
to dehne the hamiltonian theory. The initial conditions will be set in examples discussed in 
this paper by canonical light-front hamiltonians and the corresponding counterterms. The 
difficulty of the renormalization theory is to determine the counterterms so that the ehective 
hamiltonians have well dehned limits when the bare cutoh is removed. In other words, one 
has to hnd the class of initial conditions at A = cx) which imply e-independent H\A for all 
hnite A’s when e —0. 

A general iterative procedure for calculating the ehective hamiltonians will be described 
now in complete analogy to Refs. and |^. However, instead of the iteration with two 
coupled equations for Hx and T^, we describe here an iteration based on a single equation 
for Plx with an explicit solution for Tx already built in. Simple algebra and substitution of 
Eq. (2.28) into Eq. (2.27), lead to 



f\G2Xi 



(1 - h)Q2x 


Qix, 


(2.29) 
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Equation (2.29) is the main element of the renormalization group formalism in this paper. 
Note that the right-hand side is given in terms of a commutator. Therefore, the effective 
renormalized hamiltonians contain only connected interactions. This is essential for cluster 
decomposition properties of the effective hamiltonians. 

Equation (2.29) is of the form 
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^Qx = T^[Qx]. (2.30) 

The right hand side contains terms which are bilinear in the effective interaction strength. 
The initial condition for Eq. (2.29), or (2.30), is given at A = cx): Q\=oo = Ge- 

With accuracy to the hrst order in powers of the interaction strength, Gx is independent 
of A and Gt is equal to the initial regularized hamiltonian expression one intends to study, 
denoted by Hf‘\ In this initial approximation, = fx^G\ \ where G\^ = and is 
the similarity factor calculated using eigenvalues of Gf\- forms our hrst approximation 
to the similarity renormalization group trajectory of operators Gx parametrized by A. 

Eq. (2.30) can now be written in the iterative form for successive approximations to the 
trajectory Gx- Namely, 


d ^(n+l) 
d\^^ 


= Tr[Gi 


This is a conveniently abbreviated version of 


d- ^(n+l) 
d\^^ 


r{n)p{n) 

JX »2A ) 





2X 


g(n.) 

^lA J 


(2.31) 


(2.32) 


denotes a function of z^x'^ expressed through eigenvalues of Gi^\ such as in Eqs. (2.19) 
and (2.20). The initial condition is set by Thus, the solution is 


/*oo 

= gin+i) _ 

J A 


(2.33) 


By dehnition. Goo contains the e-regulated canonical hamiltonian terms and counterterms. 
The counterterms remove the part of the integral in Eq. (2.33) which diverges for hnite A 
when e —>■ 0. This dehnition is based on the requirement that matrix elements of the 
hamiltonian of the ehective theory have a limit when e is made very small. The condition 
that the necessary Goo exists is the hamiltonian version of renormalizability. It does not 
require the number of counterterms to be hnite, although a hnite number has the clear 
advantage of simplicity. 

Denote the part of the integrand in Eq. (2.33) which leads to the divergence by [G^'^] 

and the remaining part by contains the initial regulated hamiltonian 

L conv 

terms and counterterms. We dehne the counterterms in from the condition that ma¬ 

trix elements of the corresponding ehective hamiltonians with hnite A, i.e. of 
are independent of e when e —0. This condition amounts to the requirement that the co¬ 
efficients of products of creation and annihilation operators are hnite for hnite free energies 
corresponding to the momentum labels of the operators in the product. 

Note that and it is not necessary to know to calculate 

One calculates after is made independent of e when e —>• 0. 


- div^ 
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The way of deriving from divergences in the integral on the right-hand side of Eq. 

(2.33) is following. 

The diverging dependence of the integral on e when e —0, is typically of the form e~^ 
or loge times operator coefficients. The operator coefficients can be found by integrating 
the diverging part of the integrand from some arbitrary hnite value of A, say Aq, to inhnity. 
The divergence originates from the upper limit of the integration and it is independent of 
Aq. The remaining hnite part of the integral is sensitive to the lower limit of integration and 
depends on Aq. The counterterm does not depend on Aq but it contains an arbitrary hnite 
part which emerges in the following way. 

The counterterm subtracts the diverging part of the integral. But subtracting terms with 
diverging functions of e times known operators does not tell us what hnite parts times the 
same operators to keep. Thus, one needs to add arbitrary hnite parts to the numbers 1/e 
and log e in the counterterms. These hnite parts are unknown theoretically and have to be 
htted to data. In particular, observed symmetries may impose powerful constraints on the 
hnite parts. 

The diverging part of the integrand may be such that the lower limit of integration 
produces the same operators as the upper limit but the diverging numbers such as e~^ or 
loge from the upper limit are replaced by hnite numbers at the lower limit. Those hnite 
numbers depend on Aq but, once they are replaced by the required unknown hnite parts, 
one obtains a valid expression for the counterterm. The replacement can be achieved by 
adding to the integral the same operators multiplied by the numbers which are equal to 
the diherences between the unknown numbers and the numbers resulting from the lower 
limit of the integration. Thus, the unknown numbers we need to add to the integral of the 
diverging part of the integrand from Aq to inhnity, depend on Aq. We can write which is 

independent of Aq as + • The free hnite parts of the counterterms 

are contained in Oi'finlei^o) and one can ht them to data using ehective hamiltonians at some 
convenient scales A. 

More than one scale A may become necessary for accurate determination of the free 
parameters when their values have to be of considerably diherent orders of magnitude and 
require knowledge of physical phenomena at diherent scales. For simplicity, let us assume 
here that a single scale A = Aq is sufficient for practical calculations. In this case, the 
renormalization conditions are set through the description of selected physical phenomena 
using One may also consider renormalization conditions for parameters in Hx^ which 
are set at some nearby scale Ai 7 ^ Aq. This will be illustrated in the next Section. 

The complete recursion including construction of counterterms in Eq. (2.33) is given by 

eO” = d/£.l(V) + /f (2.34) 

In the limit n —cx), if the limit exists, one obtains 

p\ POO 

Gx = GEfinite{^o)+ ds[Ts[Gs]]div~ ['^s[Gs]]conV (2.35) 

J Xq J X 

Hx is obtained from Eq. (2.35) through the replacement of q^c by qx (to obtain Gx) and 
action of on Gx- These steps complete our general dehnition of the renormalized ehective 
hamiltonians in the light-front Fock spaces. The structure of the formalism is general enough 
for extension to other applications but those are not discussed in this paper. 
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The renormalized effective hamiltonians may still depend on the infrared regulators in the 
case of theories with massless particles. The limit of the infrared regulators being removed 
(5 —0 ) may lead to infrared singularities which cannot be avoided without solving for 
eigenstates of the effective hamiltonians. It depends on how the effective hamiltonians couple 
Fock sectors with many massless particles in the eigenvalue problem. For example, gauge 
theories contain massless spin-1 bosons with singular dependence of the polarization vectors 
on the small longitudinal momentum. This region contributes infrared cutoff dependent 
interaction terms. To build a bridge between the renormalized effective hamiltonians in the 
light-front Fock space and the formal held theory with a highly symmetric local lagrangian, 
such as for gauge theories, one has to understand what conditions justify calculations of 
the effective hamiltonians in separation from a nonperturbative solution of the eigenvalue 
problem when the infrared cutoff is important. For example, a neutral state of localized 
charged fermions will not radiate bosons with wavelength much larger than the distance 
between the fermions. A precise general answer to the question how the infrared cutoff 
disappears from the hamiltonian spectrum for gauge theories is not known to the author. 
Such infrared factorization property in the hamiltonian approach requires extensive studies 
in gauge theories. The infrared regularization and renormalization conditions in second order 
perturbation theory will be discussed in the next Sections. 

Our description of the perturbative calculus for renormalized effective hamiltonians starts 
here. The key observation is that one can calculate the rate of change of Q\ by expanding 
it into a power series in the effective interaction Q 2 \ at the same running scale A. This is 
obtained by repeated application of Eq. (2.29). One rewrites Eq. (2.29) as 


dGx 

dX 


f\G2X, {-fxG2x}g^^ 


+ 






(2.36) 


Then, one can substitute in the last term on the right-hand side of Eq. (2.36) by the 
preceding terms on the right-hand side. Two successive substitutions produce an expression 
for with four explicit powers of the effective interactions and the remaining terms are of 
higher order [note that (1 — fx)Gix = 0]- 


^ = 1/6, {-/'6}1 + [/6, {(1 - /) 1/5, {-/'6}1}1 
+ 1/6. {(1 - /) 1/6, {(1 - /) 1/6, {-/'6}|}]}] + o(e=). (2.37) 

We have omitted subscripts 2, A and Qix on the right-hand side. All the subscripts should 
appear in the same pattern as in Eq. (2.36). Correspondingly, the inhnite chain of substi¬ 
tutions produces an expression ordered by explicit powers of the effective interactions, to 
inhnity. 


d 

dX 


OO 

6 a = E [/6 . ({(1 - /) [/6, )<”’{-/'6 }(]}) 


72=0 



(2.38) 


The round bracket raised to the n-th power means n consequtive repetitions of the symbols 
from within the round bracket. The subscripts are omitted for clarity as in Eq. (2.37). 

The above expansion in powers of the effective interactions provides a systematic order 
by order algorithm for building an expression for the effective hamiltonian. The energy 
denominators and functions fx are calculated using eigenvalues of Qix- Therefore, in Eqs. 
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(2.37) and (2.38), the derivatives of the function fx contain two kinds of contributions: those 
resulting from differentiating the explicit A dependence in the arguments (for example, 

in Eq. (2.18) ), and those resulting from differentiating the free energy eigenvalues (for 
example, in Eq. (2.18)). Since the free energy terms include effective masses which 

depend on the width A, the derivatives of the effective masses appear in the equations on the 
right-hand side. Moving them to the left-hand side leads to coupled nonlinear differential 
equations for the effective hamiltonians. 

The general iterative approach in Eq. (2.34) or the expansion in Eq. (2.38), can be 
analysed using expansion in the running coupling constants. One can divide Gi\ into two 
parts: one which is independent of the coupling constants and another one which vanishes 
when the coupling constants are put equal to zero. The parts depending on the coupling 
constants are moved over to ^2A and treated as an interaction. After Qix is reduced to the 
part which is independent of the interactions, the derivatives of fx in Eqs. (2.36) to (2.38) do 
not introduce additional powers of the interaction strength and the series is strictly ordered 
in powers of the interactions according to their explicit appearance in the formula (2.38). 
This series then provides the perturbative expansion in terms of the running couplings. 

The simplest case of the perturbative expansions is the case where one expands in powers 
of a single coupling at a single scale. Firstly, one expands the renormalization group equations 
into a series of terms ordered by powers of the bare coupling qq. Secondly, one evaluates the 
effective coupling gi at the chosen scale Ai as a power series in the bare coupling. Thirdly, the 
latter series is inverted and the bare coupling is expressed as a series in the effective coupling 
gi- Then, one can pursue perturbative calculations in terms of the effective coupling. In 
particular, one can reduce the hamiltonian width to A2 < Ai and calculate g 2 as a series in 
gi- Such steps can be repeated. For example, one can reduce the width in each step by a 
factor 2. [Q] ^ steps will reduce the width by the factor 2“^. This way one can build the 
renormalization group flow indicated by the chains of small arrows in the diagram discussed 
in Section 1. 

2.b Regularization 

A canonical bare hamiltonian obtained from a local held theory is divergent. This Section 
describes how the ultraviolet singularities in the canonical hamiltonian are regularized with 
the bare cutoff e. We also discuss the infrared regularization. Our presentation is ordered 
as follows. First, we briehy explain the relation between the ultraviolet and infrared regu¬ 
larizations. Then, we proceed with dehnitions of the canonical hamiltonian terms. For that 
purpose, we have to discuss the fundamental set of scales in the hamiltonian approach and 
explain the role of classical held theoretic lagrangians in the construction. Then, we describe 
details of the ultraviolet and infrared regularizations. 

The ultraviolet and infrared regularizations in the light-front hamiltonians are closely 
related. The infrared structure is inhuenced by the masses in the initial hamiltonian 

hw 

carries the superscript 0 to indicate that it is the initial hamiltonian from the previous 
Section - it does not include the counterterms which need to be calculated. 

An initial mass value is generically denoted by A light-front energy, of a free 

particle with a four-momentum pm = = (p^^ + mf)/p'^^i tends to inhnity when 

m? > 0 and p'^ tends to zero. But may be hnite or even approach zero in this limit if 
m? = 0 and p^ approaches zero too. The limit of small momentum is always a high- 
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energy limit when > 0. But it ceases to be the high-energy limit for the small transverse 
momenta when 0. Thus, in the initial hamiltonian is capable of switching from 

the high-energy regime in the longitudinal direction to the low-energy one when we take the 
limit ^ 0. Conversely, introducing masses turns the infrared low-energy regime into 

the high-energy regime. 

We begin the construction of with an enumeration of momentum scales which are 
involved in it. We distinguish a set of fundamentally different scales related to the boundary 
conditions for helds at spatial inhnity, small momentum cutoffs, phenomenological parame¬ 
ters and large momentum cutoffs. 

The bare hamiltonian is dehned in terms of the operators q^o- We consider quantum 
helds to be built from these creation and annihilation operators. We consider quantization 
rules for classical held theory as secondary to the quantum theory of particles (see also Ref. 

HD- 

The initial basis in the Fock space is built from the vacuum state |0) by successive actions 
of the creation operators goo- We will need more detailed notation. Fermion, antifermion 
and boson creation and annihilation operators are conventionally denoted by d\ a\ b, 
d and a, respectively. For example, \ka >= &l,a.|0) denotes a state of one bare fermion of 
momentum k = (fc’*', k^). Spin z-axis projection, havor, color or other quantum numbers, are 
denoted by a common symbol a. The momentum variables in the subscripts of creation and 
annihilation operators are distinguished in order to set the scales involved in the dehnition 
of the hamiltonian. The momentum scales will be ordered from the smallest to the largest. 
The order of scales in momentum space will be reverse to the order of scales in position 
space. 

The largest scale to be considered in the position space is the quantization volume. 
This volume will never be allowed to enter any physical consideration. The reason is that 
we want to describie phenomena which we believe to be independent of any conditions at 
remote distances. If this belief is not justihed our formulation will be invalid. 

We assume the quantization volume to be so large that the boundary conditions are 
of no importance. Since we dehne states of bare particles by assigning them momenta, 
two different values of momenta are separated by at least one unit of the inverse size of 
the quantization volume. By assumption, no physically interesting question will concern 
wavelengths comparable with that size. 

The Poincare symmetry of our theory will hold only approximately because the space of 
momentum variables in our hamiltonian quantum mechanics will not be continuous. But 
we shall insure by our choices of scales that the granulation in momentum space will never 
be noticeable and the quantization volume will be effectively inhnite for all our purposes. 
Thus, we will universally adopt continuous notation for momentum variables. For example, 
invariance under boosts along 2;-axis will only mean that there is sufficiently many different 
momentum states so that a change of k'^ by a prescribed factor can be approximated using a 
change by a rational factor with better accuracy than any gedanken experiment could verify. 
The same concerns rotations. 

Our assumption that the quantization volume can be large enough for its size to become 
irrelevant to physics, seems to contradict potential relevance of topological features which 
can be related to the boundary conditions in classical held theory. This is not necessarily the 
case. The place where a large volume boundary conditions can appear in the hamiltonian 
theory is in the formation of dynamically stable infrared components in physical states. 
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These might be called “wee”. Expectation values of field operators in states containing 
many wee particles may be responsible for producing the approximate long distance or large 
correlation length classical held dynamics. 

Spontaneous symmetry breaking will be allowed to appear through the wee components 
but it is not known in the quantum theory that such components do emerge in physically 
interesting cases. It is also not known how many different long wavelength structures may 
form. Most importantly, we do not know the saturation mechanism which leads to the 
formation of such stable components. Therefore, their scale, even if present, escapes our 
control at this time and requires more studies. 

We cannot exclude that some topological features will eventually result from our as¬ 
sumptions. But we cannot address this issue before calculating the effective hamiltonians 
and their eigenstates. The issue will have to be considered if it turns out that the large scale 
features do matter in the formation of the eigenstates. 

Some hypotheses about zero-modes and spontaneously broken symmetry in light-front 
quantum held theories were recently discussed in Ref. which quotes important earlier 
literature on the subject. Basically, one may expect that new terms emerge in the effective 
hamiltonians and they account for the large scale dynamical effects. Susskind and co-workers 
have proposed a way to think about the wee parton dynamics in a model. |T^ 


Ref. 13 


discusses the QCD sum rules techniques. However, the original quantum dynamics of the 
vacuum formation and spontaneous symmetry breaking in the present approach are not yet 
understood and cannot be discussed further here. 

The next smaller size in position space is the inverse of the infrared regulator. Three 
main types of the infrared regularization appear in the light-front dynamics. The hrst is a 
straightforward lower bound on /c+, denoted by The second is a lower bound, denoted by 
6, on the parent -|—momentum fraction that can be carried by a particle in an interaction. 
The third is a mass parameter /r^. fis appears as the mass parameter in the initial 
hamiltonian It is introduced for massless bare particles, fis may effectively cut off the 
small longitudinal momentum region at a small scale order /r^/A, where A is some ultraviolet 
momentum upper bound. When the infrared regularization is being removed, 6^, 6 or 
are sent to zero. 

One should note that 5 is a dimensionless number while 6^ and have dimensions. 
In order to make a connection one needs to introduce a number with the dimension of a 
mass. It is assumed that the hnite detector sensitivity and experimental wave packet widths 
are much larger than the infrared regularization scale and the wide gap between the two 
scales lets us choose the dimensionfull number in such a way that physical quantities will be 
independent of the small infrared cutoff scale. 

One should also note that and 6 distinguish the longitudinal direction from the trans¬ 
verse directions on the light-front while the mass regulator combines them. The mass regu¬ 
lator seems to avoid the distinction. However, introducing a mass term for photons or gluons 
which have only two polarizations does explicitly distinguish directions. In the gauge theo¬ 
ries, understanding the mechanism by which the infrared regularization mass scale becomes 
irrelevant to physics is one of the main theoretical challenges in the hamiltonian theory. 
Apparently, the lack of sensitivity to the infrared cutoffs which distinguish directions is re¬ 
lated to the rotational symmetry and to the gauge symmetry. The extraordinary interest 
in understanding the mechanism of this insensitivity, if it exists, stems from the fact that it 
must combine the hamiltonian quantum mechanics with the Lorentz and gauge symmetries. 


24 






The next smaller scale in position space is set by the size of the volume used for prepa¬ 
ration of incoming and detection of outgoing particles (including bound states) and the 
corresponding time scale. Physics to be discussed in our approach will be contained within 
this scale. Wave packets mean values will be contained within distances on that scale and 
the momentum space widths of the wave packets will be limited from below by its inverse. 
Physical observables are allowed to depend on this scale since the preparation and detection 
of states is a part of a physical process. 

The order of magnitude of momenta larger than the experimental wave packet widths will 
be characterized initially in terms of three different scales. Namely, (1) masses of particles, 
(2) the width of the effective hamiltonian (i.e. A), and (3) the bare cutoff scale e~^. Later, 
when solving for the hamiltonian spectrum, a new scale may emerge dynamically, determined 
by the effective coupling constants, masses and width of the effective hamiltonian. Scale 
invariance at large momenta may be violated through a dimensional transmutation even if 
all mass scales are negligible in comparison to the momenta and A. 

The masses of particles appear in the initial bare hamiltonian. Renormalization effects 
lead to different renormalized masses in the effective hamiltonians. The effective masses will 
depend on the width A and the strength of the interactions. 

The width A will be ranging from cx) to some convenient finite values. Efficient description 
of physical phenomena which involve energy-momentum transfers of the order of k will 
require A to be larger than k. It will also be useful to use A not too large in comparison to k 
in order to avoid too much detail in the effective hamiltonian dynamics. Useful values of A 
in nonrelativistic systems are smaller than the effective masses. For example, in QED, the 
convenient effective cutoff A on the electron momentum in a hydrogen atom is much larger 
than the electron binding energy and much smaller than the electron mass. 

The bare cutoff scale e~^ will be related to much larger momentum scales than any scale 
of interest in physics. It should be stressed that the formal limit A —>■ oo is used here 
only to remove A dependence from the hamiltonian regulated by e. In other words, no A 
dependence appears in the hamiltonians with A larger than the scale implied by e~^. No 
physical quantity will be allowed to depend on e when e —> 0. The similarity renormalization 
scheme for hamiltonians is built to achieve this goal to all orders in perturbation theory (cf. 

i 0)- 

One would have to make additional subtractions if dependence on e were recovered in the 
process of Ending eigenstates of the efiective hamiltonians beyond perturbation theory. Such 
unfavorable result would indicate that perturbation theory in the running interactions is not 
reliable in deriving the efiective hamiltonians. Although this possibility will be neglected in 
the present paper, we mention it because we lack understanding of the mechanisms of the 
ground state formation and the insensitivity to the infrared regularization. The point is that 
the order of removal of the ultraviolet and infrared cutoffs may matter. The infrared regular¬ 
ization by fis turns the infrared divergences into the ultraviolet ones. Then, the ultraviolet 
renormalization can remove the divergences by removing the e-dependence from the efiective 
hamiltonians. But the limit of /i^ —0 destroys this correlation of the ultraviolet and infrared 
divergences. Therefore, it is not excluded that the removal of the infrared regularization may 
lead to different results depending on the order of limits one takes. Such considerations are 
beyond the scope of the present article. This remark concludes our description of scales 
involved in the construction of hamiltonians. 

We proceed to the explicit construction of simplest terms in the hamiltonian . The 
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counterterms are not known exactly from the outset. Light-front power counting rules are 
helpful in determining the structure of the counterterms but more detailes are required 
in practice. The similarity renormalization group is to provide the required insight. 

All starting hamiltonians we consider will contain a free part which we denote by 
The free part for fermions and bosons has the form 



k + m, 


( 0)2 


+ d\^dka) + 






(2.39) 


We adopt the following conventions. Summation over a denotes a sum over all quantum 
numbers except the momentum. 


Ibvr^ 


dk'^ 




d^k 


(2.40) 


The creation and annihilation operators in Eq. (2.39) are the bare ones commonly denoted 
in Section 2.a by goo- They satisfy standard commutation or anticommutation relations 


O'kcr-, 




(2.41) 


with all other commutators or anticommutators equal zero as dictated by the spin and 
statistics assignments of Yukawa theory, QED or QCD. 

The initial mass parameters and do not include effects of any interactions and 
are independent of the interaction strength. We may have to consider limits where the mass 
parameters are close to zero, in comparison to all other quantities of relevance to physics. For 
example, may be the infrared regulator mass denoted by Recall that the subscript 
e indicates that the mass parameters stand in the hamiltonian with \ = 00 . 

The initial hamiltonian contains an interaction part, — Qi‘\ For example, 

electrons may emit photons. One writes the corresponding interaction term in QED in the 
form 


E / [ki][k2][k[]lQnH\ki + fca - (2.42) 

CTlO'2 0'i 

We use here conventions to be specihed shortly. The hamiltonian term (2.42) is contained 
in the expression 


x+=0 


h = j dx d^x^ e'lp^wi^x) 4 L{x)iIj^(o){x 
where the helds 'ip^(Q){x) and A'^{x) for = 0 are dehned by writing 


ikx 


(2.43) 


(2.44) 


and 




(2.45) 
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Spinors Umka and Vmka are defined by boosting spinors for fermions at rest, Uma and Vmai 
to the momentum A;, as if the fermion mass were m. This is done using the light-front 
kinematical boost representation for fermions 


S{m, k) = (mk^) -|- A_(m -|- a'^k'^)]. 


(2.46) 


Namely, Umka = S{m, k)uma and Vmka = S{m, k)vma- Solving constraint equations for the 
free fermion helds in canonical held theory amounts to using these spinors. The same boost 
operation dehnes the polarization vectors for photons which are independent of the photon 
mass. We have Ska = (^ta — i^ka = The spin label a denotes the 

spin projection on the z-axis. We adopt a number of conventions from Ref. [|T^. It is useful 
to work with the above spinors and polarization vectors because they provide insight into 
the physical interpretation of the calculated matrix elements. For example, the spinors and 
polarization vectors help in tracing cancelations which result from the current conservation 
(e.g. see Eq. (3.103) etc. in the next Section). 

Equation (2.43) includes 5 well known terms in addition to (2.42). The other 5 terms 
lead to emission of photons by positrons, absorption of photons by electrons or positrons, 
or to transitions between electron-positron pairs and photons. There is no term leading to 
creation of an electron-positron pair and a photon, or to annihilation of such three particles. 
This is the distinguished property of the light-front hamiltonians: conservation of momentum 
k~^ > 0 excludes a possibility that the three momenta sum up to zero. 

Strictly speaking, one has to limit each k~^ from below by a nonzero positive lower bound 
in order to make sure that the three -|—momentum components cannot add up to zero. 
This lower bound is provided by the inverse of the quantization volume. Our approach will 
ensure that this largest of spatial scales in the theory does not need to be invoked in the 
description of physical phenomena. The regularization procedure will cut off such small 
momenta long before they will have a chance to become relevant. If high-order perturbation 
theory subsequently leads to effective hamiltonians which describe universal low momentum 
components in all physical states then, the notion of a nontrivial vacuum will have to be taken 
seriously into account for practical computational reasons. A priori, we cannot exclude this 
will happen. But we postpone considerations of such a situation until it becomes necessary 
in the future work. 

The product "0 denotes a sum of 6 basic interactions. The products of creation and 
annihilation operators are ordered as indicated at the beginning of this Section. However, 
Eq. (2.43) requires additional steps before one can assign it a well dehned meaning because 
operators such as (2.42) can easily produce states of inhnite norm. In the light-front hamil- 
tonian approach, one needs to dehne the individual terms such as (2.42) in order to provide 
meaning to the whole combination of similar terms in Eq. (2.43) 

First of all, there are inverse powers of k~^ in Eq. (2.42) and k~^ may be arbitrarily close to 
0. For example, when k^ and kf in (2.42) are similar (and they are allowed to be arbitrarily 
close to each other no matter what their own size is), the photon momentum k^ = k^ — fcj*" is 
arbitrarily close to zero. The problem is that the photon momentum appears in the photon 
polarization vector in the denominator: . Unless k^ is close to zero, the 

resulting emission strength will approach cxo for k^ —> 0. Therefore, even for a very small 
coupling constant e, the interaction can be arbitrarily strong. There will be reasons for a 
cancelation of this divergence in special circumstances. For example, in the tree diagrams 
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for the S-matrix elements in QED, the cancelation will be a consequence of the presence of 
more terms in the hamiltonian and the energy conservation in physical processes. However, 
for the off-energy-shell matrix elements of the T-matrix, in loop diagrams, or in bound state 
equations, such cancelations will not be ensured automatically and could lead to untractable 
expressions. 

In particular, one has to keep in mind that in the perturbative calculation of the S-matrix 
it is possible to apply energy and momentum conservation laws for incoming and outgoing 
particles on their mass-shells. In contrast, in the bound state calculations, the individual 
particle momenta cannot simultaneously be on the individual mass-shells and still sum up to 
the bound state momentum - the bound state dynamics is always off-shell and perturbative 
mechanisms for cancelations cease to be sufficient. 

In Eq. (2.42), the inverse powers of the longitudinal momentum also appear in the 
fermion spinors. These can be a source of divergences, too. However, the examples we 
describe in this article do not lead to problems with infrared fermion divergences and we will 
not dwell on this subject here. 

Secondly, the spinor matrix elements depend on the transverse momenta of the fermions 
and the boson polarization vector depends on the transverse momentum of the boson. The 
strength of the interaction grows when the relative transverse momenta grow. Again, even 
in the case of a very small charge e, divergences arise and the question is if the interaction 
is hnite in any way. 

For example, the hnitness problem manifests itself clearly when one attempts to evaluate 
the ratio of norms of the states h\ka) and \k(j). This ratio is certainly not hnite and it is 
ill-dehned. 

Therfore, one might ask if it is useful to consider this ill-dehned hamiltonian term. The 
answer is unambiguous: yes. The reason is that the scattering amplitudes calculated using 
this term in combination with two other terms in second order perturbation theory, agree 
very well with observable scattering of electrons and photons. No loop integration appears 
in these calculations to indicate the divergence problem. 

It is well known that the terms one should put into the light-front hamiltonian are 
provided by the formal lagrangian density for electrodynamics C = 

One can rewrite the lagrangian density into a corresponding light-front hamiltonian density 
by using an expression for the energy-momentum tensor density Integrating T~^~ over 
the light-front hyperplane gives the expression one starts from in building the light-front 
hamiltonian for QED. 

The initial hamiltonian for QED results from formal operations on helds ^jJ+ and A-^. 
One uses the gauge A+ = 0 and solves the constraint equations, substitutes expansions 
of the form (2.44) and (2.45) into the formal expression for T~^~, integrates the density over 
the light-front hyperplane and normal-orders all terms. The normal-ordering produces terms 
that involve numerically divergent momentum integrals. The classical held theory does not 
tell us what to do with the divergences resulting from the ordering of operators. 

To deal with all divergences one has to regularize the theory from the outset. The 
naive connection between the classical theory and the quantum theory as given by the 
quantization rules is broken. The regularization turns out to introduce large terms into 
the initial hamiltonian. To gain control on the regularization effects one has to construct 
a renormalization theory for hamiltonians. The diverging terms which result from normal 
ordering can be safely dropped in the form they appear ill-dehned in the canonical approach. 
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because the renormalization procedure introduces other terms of the same operator structure 
to replace them. 

Reasons to adopt the perturbatively and semi-classically motivated strategy of building 
the initial hamiltonian from a classical largrangian are less clear in the Yukawa theory and 
QCD than in QED. 

Application of the pseudoscalar Yukawa interaction in nuclear physics suggests that the 
coupling constant for nucleons emitting pions is order 10. Perturbation theory is question¬ 
able. However, it is important to realize that the vertices for pions and nucleons which are 
used in the relativistic nuclear physics involve pion-nucleon vertex form factors. The vertex 
form factors are analogous to the similarity functions fx- Thus, the well-known pion-nucleon 
interaction could be identihed with an interaction term in an effective hamiltonian of an 
unknown original theory rather than the initial bare hamiltonian term. 

Following this observation, one can try to calculate an effective hamiltonian in a pseu¬ 
doscalar Yukawa theory for A on the order of the cutoff parameters known in nuclear physics 
and ask what terms and what sizes of the coupling constants would have to appear in that 
hamiltonian to match experimental data, if it is possible. [1^ The narrow similarity vertex 


factors can, in principle, weaken the effective strength of the interaction by reducing its 
width in momentum space to a very small range so that the coupling constant itself being 
order 10 may be too small to invalidate perturbation theory in the whole interaction. Then, 
one could study meson-barion physics using perturbation theory with such narrow effective 
hamiltonians. The goal would be to correlate the vertex form factors parameters with the 
coupling strength and masses of mesons and barions which are used in the effective theory. 

Thus, it becomes interesting to see what kind of effective hamiltonians result from an 
abstract initial obtained from the Yukawa lagrangian with a very small coupling constant 
or from similar held-theoretic lagrangians, such as for a-models. For one could see this 
way examples of the small width effective hamiltonian structures implied by the local held 
theory and one could study off-shell mechanisms which may lead to covariant results in the 
renormalized hamiltonian approach. The size of the coupling is of secondary importance 
from that point of view. Moreover, the range of converegence of expansion in the coupling 
constant in the effective theory is not known and it should be studied because it may be 
large in the small momentum width effective hamiltonians. 

One should also keep in mind that the effective hamiltonians may have a number of 
universal chracteristics which are independent of details of the initial theory. For example, 
there exists some resemblance between effective low-energy interactions in theories with 
pseudo-scalar direct and derivative couplings between bosons and fermions. 

In QCD, one proceeds by analogy with QFD assuming that the constituent quark picture 
of hadrons approximates the solution to QCD. 

On the one hand, conhnement invalidates close analogy with QED. On the other hand, 
asymptotic freedom makes the limit of bare cutoffs being removed more plausible than in 
QED. One hopes that the strength of effective interactions never grows too large when the 
scale A is lowered and the effective hamiltonian for the constituent quark picture for hadrons 
may be obtained. Again, the important feature of the similarity factors is that the strength 
of the interactions is not given solely by the size of the coupling constant. The width of the 
function fx is also important. The smaller is the width the smaller is the strength of the 
interactions. Therefore, the running coupling constant may still he within the perturbative 
domain for the evaluation of the effective hamiltonians of reasonably small widths despite 
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the coupling constant growth when the width gets smaller. 

To explain the regularization for light-front hamiltonians which we apply to expressions 
resulting from field theoretic lagrangian densities, we hrst discuss the term (2.42). In that 
term, the parent momentum P equals k[. The spinors and polarization vectors conveniently 
group a number of terms with different momentum dependences into a combination which is 
invariant under the light-front kinematical symmetry transformations. Among those terms, 
there are terms containing masses, terms which in held theory result from derivatives id-^ 
or id~^, or from inverting the operator id~^. All those derivatives are replaced in the term 
(2.42) by momenta of particles created or destroyed by that term. 

We hrst introduce the daughter momentum variables for the created electron and pho¬ 
ton. We have briehy introduced daughter momentum variables in a similar conhguration in 
Section 2.a when dehning the similarity functions fx. Here, we use the daughter momentum 
variables for the purpose of regularization. The variables are 


xi = ki /k'l = X, 

(2.47.a) 

X2 = k2 /k'l = 1 — X, 

(2.47.6) 

x[ = k't/k't = 1, 

(2.47.C) 

K'l; = kj^ — XiP'^ = 

(2.47.d) 

/^2 ^ k^ X^P ~~ K , 

(2.47.e) 

k'I; = k'l; — x'-^P^ = 0. 

(2.47./) 


For each of the creation or annihilation operators in the interaction term (2.42), we dehne a 
daughter energy variable. Namely, 



^ 2-2 -|- 77^(0) 2 

Xi 

X 

42 + ^( 0)2 

_ «^2 + 4)2 

X‘2 

1 — X 


(2.48.a) 


(2.48.6) 


I ™(0)2 

Xi 


(2.48.C) 


For each creation and annihilation operator in the interaction term (2.42), we introduce a 
factor which is a function, r{yi), of the variable yi = ee^/A^, where the snbscript i denotes 
the operator in question. In the no cntoff limit, e —0. A is an arbitrary constant with 
dimension of a mass {h = c = 1). All masses and momenta are measnred in nnits of A. In 
this article, we choose r{y) = (1 -|- y)~^. Thus, the term (2.42) is regulated by introducing 
the factor 


(1 + eei/A2)-i(l + ee2/A2)-i(l + ee^/AY' (2-49) 
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under the integral. The third factor in the above expression can be replaced by 1, since 
is a hnite constant and it cannot compensate the smallness of e. We shall make such 
replacements wherever the parent momentum is carried by a single creation or annihilation 
operator. 

An interesting feature arises in this regularization scheme. Namely, if there is no infrared 
divergences due to the small values of x, one can remove masses from the regularization 
factors. It is the small x divergences that invite keeping the masses in the daughter energy 
so that the latter grows towards small x even for the zero transverse momenta. 

In the case of terms which contain only 1 creation and 1 annihilation operator, i.e. in Qi’\ 
no regularization is introduced. This is important because any restrictions on the particle 
momenta in these terms would violate kinematical symmetries of the light-front hamiltonian 
dynamics. The reason is that the particle momenta in these terms are equal to the parent 
momenta and limiting the parent momenta violates the light-front boost invariance. 

In the initial expressions for hamiltonian densities of Yukawa theory, QED or QCD, only 
terms with products of up to four helds appear. Therefore, we have only two more situations 
to consider in addition to the cases such as and terms of the type (2.42). In the hrst 
situation we have three creation operators and one annihilation operator or vice versa, and 
in the second situation we have two creation and two annihilation operators. 

Independently of the number of creation and annihilation operators in a product, the 
regularization is introduced by multiplying every creation and annihilation operator in the 
product by a function riy) such as in the factor (2.49), where y = ee^/A^ and is the 
corresponding daughter energy variable. Later, after counterterms are calculated, the same 
regularization factors are introduced in the counterterms. 

An additional step is required in the case of hamiltonian terms which originate from the 
products of four helds including inverse powers of id^ acting on a product of two helds. This 
additional step will be described now. 

In the terms containing the products of four helds and inverse powers of id^, we introduce 
two kinds of a hfth daughter momentum and two corresponding daughter energy variables, 
6512 and 6534 . The numbering originates from assigning numbers to the helds in the product 
according to the schematic notation 0 i 02 (* 5 ’'")~” 0304 - One of the hfth daughter energy 
variables is associated with the operators coming from the helds number 1 and 2 , and the 
other one is associated with the operators coming from the helds number 3 and 4. The 
regularized terms will contain an additional product of functions ^( 1 / 512 ) and with the 

arguments 1/512 = 66512 /A^ and 7/534 = 66534 /AT 

The auxiliary daughter energy variables 6512 and 6534 are calculated as if they represented 
daughter energy variables for an intermediate particle, a boson or a fermion, created and 
annihilated in the vertices which contained the products 0 i 02 and (/> 304 , respectively. Those 
vertices are treated as if each of them contained three helds instead of two but the held of 
the intermediate particle was contracted so that the corresponding creation operator and 
the corresponding annihilation operator are absent in the resulting term. This particular 
dehnition of a gedanken intermediate particle does not refer to any particular Fock states 
and remains valid when the operators goo are replaced with q\ by the unitary transformation 
Ux- The dehnition was inspired by Refs. [Q and m where the correspondence between 


the intermediate states with backward moving particles with spin in the inhnite momentum 
frame and the light-front seagull interaction terms is extensively described. 

Mathematically, the dehnition of 6512 and 6534 is introduced in the following way. Every 
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creation and annihilation operator in the helds 0i, 02 , 03 and 04 is assigned a corresponding 
number Sj, i=l,2,3,4. s* equals +1 for a creation operator and Sj equals —1 for an annihilation 
operator. We dehne k'^ = + sji'l\ and S 5 = {—s^.k'^ — s^k'l)/k'^. The gedanken particle 

is thought to be created in the product of helds including 0304 when S 5 = 1 and it is thought 
to be annihilated in that product when S 5 = —1. We dehne the momentum k^ = 
by the relation S 5 /C 5 = —S 3 /C 3 — S 4 /C 4 = Sifci + 52 ^ 2 - We also introduce two auxiliary parent 
momenta, P 34 = ^(^5 + ^3 + k^) and P 12 = ^{k^ + ki + ^ 2 )- Then, we introduce the daughter 
momentum and energy variables 


where m 


( 0 ) 

e5 


2^512 = 

(2.50.a) 

^512 ~ ^5 ~ ^ 512 Pi 2 i 

(2.50.6) 

±2 , (0) 2 
«512 + K 5 

6512 — , 

2^512 

(2.50.C) 

3^534 = k+/P+, 

(2.50.d) 


(2.50.e) 

±2 , (0) 2 
'^534 + ^5 

6534 — , 

2^534 

(2.50./) 

for regularization of the terms involving 

and equals 


This step completes our dehnition of the 


for regularization of the terms involving 
ultraviolet regularization of initial hamiltonians. 

We proceed to the dehnition of the infrared regularization. Inverse powers of id~^ for 
massive particles are already regulated when the ultraviolet regularization is imposed. This 
was explained above. 

For each creation and annihilation operator of an initially massless particle, we introduce 
a factor which limits the daughter momentum fraction x for that operator to be greater than 
5. An example of such a factor is given by (1 + 5/x)~^. Note that our dehnition also implies 
that the same regularization factor is inserted for the gedanken particles with 0:512 dehned 
in Eq. (2.50.a) and 0:534 dehned in Eq. (2.50.d). 

Our procedure of limiting the longitudinal momentum fractions in the interaction terms 
preserves the light-front boost invariance explicitly but may lead to problems in the renor¬ 
malization procedure. The problems may arise when one expands various functions in powers 
of momentum using typical light-front boost invariant combinations of relative momentum 
variables and the expansion leads to new strongly divergent terms, so that the corresponding 
chain of counterterms will be too hard to calculate and the ehective hamiltonians will not 
converge. [^] The argument is not iron-clad: it is not proven that such a situation has to 
arise. 

We suggest that the way out of the potential problem can be looked for in the expansions 
around nonzero values of x so that terms such as /x will not combine transverse singular¬ 
ities at large k with longitudinal singularities at small x. At this stage, we can only assume 
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that the expansions around the right values of the longitudinal mometum fractions, such as 
X = 1/2 for the quarks in heavy quarkonia or electrons in positronium or, such as x = a 
for the gluons or photons, remove the danger of creating new divergences in the high order 
perturbation theory. 

On the other hand, when one limits p'^ from below by 6 ^, instead of the fractions, 
to exclude appearance of the complicated boost invariant combinations then, the boost 
invariance is not explicitly preserved. But we have to remember that the renormalization 
procedure cannot potentially be complete without specifying renormalization conditions. 
Theoretical calculations of observables which can provide the renormalization conditions 
require solving for the eigenvalues of the hamiltonian. This is an intrinsically complicated 
nonperturbative step. The approximations which are necessary in the calculation will most 
probably violate boost invariance and this may exclude predictions for the boost invariant 
quantities of sufficient accuracy to dehne the renormalization conditions. 

Fortunately, one can preserve the explicit light-front boost invariance throughout the 
calculation of the effective interactions by using the limits on the momentum fractions. In 
this case, the renormalization conditions play no role in obtaining the light-front symmetry 
because the effective hamiltonians are explicitly symmetric. One can also make the necessary 
approximations in the diagonalization procedures without violating the symmetries. Thus, 
the exactly symmetric renormalization conditions can be imposed. 

Besides introducing the cutoff 6 on the momentum fractions carried by massless particles, 
we can also introduce for each initially massless particle a hnite regularization mass term 
which is denoted by ms- In other words, in the case of the initially massless particles, 
^( 0 ) _ Such hnite masses in the daughter energies lead to additional suppression of the 
infrared longitudinal momentum region. The reason for allowing the additional mass terms 
follows from anticipation of counterterms. 

The relative momentum cutoffs lead to the counterterms which have unknown hnite parts. 
In particular, the mass counterterms for the initially massless particles contain unknown 
hnite parts. Since, at this point, the hnite parts are arbitrary and they are not known to be 
zero, we introduce the mass terms and investigate their role. 

However, the introduction of masses for photons or gluons in Qi brings the masses into 
the denominators of perturbation theory for ehective hamiltonians and leads to additional 
singularities which do not appear in the case of massless particles. Therefore, additional 
modihcations may be required to remove such diverging ehects of the masses. For example, 
the instantaneous vector boson exchange interaction may be multiplied by the factor k^/{k‘^+ 
pif) with a suitably dehned The operational dehnition for such a term involving four helds 
(two creation operators and two annihilation operators or, three creation operators and one 
annihilation operator or, one creation operator and three annihilation operators) labeled by 
momenta fcj, i = 1,2, 3,4, is to calculate using dehnitions introduced for Eqs. (2.50.a) 
to (2.50.f) and put k equal to \k^\ with the parent transverse momentum for the whole 
interaction term, = ^{k^ + k^ + k^ + fcf), set to 0. 

If the mass parameters are treated as small perturbations, one does not include them in 
Qi- Then, they do not modify energy denominators in perturbation theory and can hopefully 
be taken into account successively in a power expansion, 

In a close analogy to the mass parameter, one can also introduce a lower bound on the 
transverse momenta. This alternative has the advantage of not modifying the free energy 
term and, therefore, does not introduce additional divergences in lowest order calculations. 
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This option is not further discussed in the present article. 

The infrared regularization completes our definition of the initial hamiltonians. 


2.C Renormalization conditions 

The free finite parts of the counterterms are determined by renormalization conditions. 
The renormalization conditions result from the comparison of theoretical predictions with 
data. However, the calculations of observables require solutions to the bound state or scat¬ 
tering problems for the renormalized hamiltonians. In principle, one could work with hamil¬ 
tonians of any width A. In practice, one is limited to consider some subspaces in the Fock 
space. Therefore, the issue of setting renormalization conditions is subtle. 

Various approximations may influence the fit to data and it may be necessary to scan 
a large number of options. If a fixed point is discovered and the bare cutoff can be made 
sufficiently large, one can pursue a search for a full set of fundamental parameters that 
characterize a theory. But we will not achieve this kind of quality in this paper. Only 
second order perturbative examples are analysed here. 

In theories without confinement, one has an option of defining on-mass-shell renormaliza¬ 
tion conditions for single particles. In the present formalism, it means that one can determine 
free parameters in the effective hamiltonians by demanding that single particle eigenstates 
and the S'-matrix have the required properties. The key examples to be discussed in detail 
here are nucleons and pions in Yukawa theory (pseudoscalar coupling) and electrons and 
photons in QED. The on-shell scattering amplitudes can be used for direct determination of 
the coupling constants, but in the second order calculations described in this paper no cou¬ 
pling constant renormalization appears because it requires higher orders of the perturbation 
theory. 

When one wants to include confinement, there is a trouble with deciding how to choose the 
values of masses for colored particles. Quarks and gluons are not supposed to be observable. 
Should or should not they have diverging infrared cutoff dependent mass terms in the initial 
hamiltonian? We do not have any definite answer to this question today. 

For the confined particles, one can adopt the strategy of fitting observables for their bound 
states (masses, form factors, decay rates, scattering cross sections). This is a very demanding 
strategy which must include solutions to a number of problems that are not solved yet. For 
example, in order to obtain a rotationally invariant spectrum for the bound states one has 
to explain how the apparent rotational symmetry violations in the hamiltonian structure 
off the energy-shell combine to produce a degenerate mass spectrum once a suitable choice 
of the finite parts of the counterterms is made. Unfortunately, this problem is not solved 
even in theories without confinement, the chief example being QED. Another condition is 
the demand that a bound state current matrix element has some definite structure which 
is required by the S-matrix theory for electroweak external probes. In fact, the present 
formalism can be described as an attempt to provide a framework in which these types of 
problems can be sufficiently well defined to seek their solutions. 

The key question we have to answer in practice is how many effective particles have to 
be taken into account in the eigenvalue problem for an effective hamiltonian and how many 
can be included in a real calculation. A good example of a theoretical problem one can think 
of is how the momentum of a proton is shared by its constituents. The phenomenology of 
deep inelastic scattering of leptons on nucleons suggests that a considerable number even at 
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moderate momentum transfers. If the number of effective constituents has to be large, one 
may encounter ambiguities in the determination of free parts of counterterms because the 
proton observables will be calculable only through complicated procedures. On the other 
hand, the constituent quark model suggests that the large number of constituents is not 
necessary for an approximate description of the proton properties. Clearly, we do not know 
yet how hard it is to hnd the hnite parts of the counterterms using bound state observables 
but we cannot exclude that it is doable. 

Qx in Eq. (2.35) contains the free hnite parts of the counterterms in the term Qefinitei^o)- 
The hamiltonian Hx^ can be used for calculating the scattering amplitudes and bound state 
properties. The most familiar example is QED. It is almost purely perturbative. One can 
calculate the physical electron energy as dehned by the lowest eigenvalue of the effective 
hamiltonian for eigenstates with the electron quantum numbers. Thanks to the symmetries 
of the light-front dynamics, the eigenvalue will have the form -|- m‘l)/p'^ and rrie will 
have to be equal to the physical electron mass. Note also that the effective mass term for 
the interacting photons will have to be different from zero (and growing with A) in order to 
obtain massless photon eigenstates. Examples of the renormalization conditions for QED 
will be presented in the next Section but no discussion of effective scattering theory beyond 
the second order perturbative results is offered in the present paper. 

When calculating the unknown renormalization parameters, one can also take advantage 
of another aspect of the present theory. The spectrum of physical masses will consist of 
multiplets, as implied by rotational symmetry, if the hamiltonian belongs to a renormalized 
algebra of Poincare generators. The Poincare algebra commutation relations should be 
studied order by order in perturbation theory to hnd out if they impose constraints on the 
free parameters. The initial algebra is spoiled only at the scale of the regularization cutoff. 
The effective algebra can be satished if the dependence on the cutoff scale is removed in the 
renormalization process. 

The general structure of the similarity transformation which is dehned in terms of the 
creation and annihilation operators, provides a natural way to extend the renormalization 
procedure to the whole Poincare algebra. The renormalization group evolution of all gener¬ 
ators is given by Eq. (2.5). The evolution preserves commutation relations. Therefore, the 
ehective generators satisfy the same algebra as the initial ones. 

The renormalization conditions which are related to symmetries can be imposed using 
commutation relations among the symmetry generators without direct evaluation of physical 
quantities. The renormalization of the whole Poincare algebra or other symmetry generators 
is not further analysed in this article. 

3. EXAMPLES OF APPLICATION 

This Section describes a set of examples of lowest order calculations of renormalized 
ehective hamiltonians in the light-front Fock space using the scheme from Section 2. We begin 
by the description of generic rules for calculating the right-hand sides of the renormalization 
group equations, i.e. Eq. (2.29). The rules follow from the commutator structure. Then, we 
discuss examples from Yukawa theory, QED and QCD. 

3.a Evaluation of commutators 
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The right-hand sides of Eqs. (2.29) and (2.38) are commntators. This implies that 
the interactions which appear in the effective renormalized hamiltonians Hx{qx) and in the 
connterterms in are connected. This Section explains how this resnlt comes abont. 

The commntators can be evalnated in a nnmber of eqnivalent ways bnt some are more 
convenient than others. Snppose we are to evalnate 

H=[A,{B](.\. (3.1) 

A = A(X, Y) n'i, 4. n/ii B = B(V, W) nfei < nSi <!», and C = E. E{z)ala.. The 
right-hand side of Eq. (3.1) eqnals 

Ia Ja TJty Ib Jb 

H=A{x,Y) n<n % E -E nn 

i=i j=i k=i i=i 

_mv^ n < li T) n 4. n (3.2) 

where Eyj = J2'i=i E{wi) and Ey = J2k=i By commnting n/=i in the hrst term 

throngh Oa-^i one generates the contracted terms with a nnmber of contractions ranging 
from 1 to the smaller of the nnmbers Ja and Ib, and a term with n/=i o-yj standing to the 
right of 11^=1 nij.- Then, by commnting ^ 10 ; in the latter term throngh one 

obtains new contracted terms with the nnmber of contractions ranging from 1 to the smaller 
of the nnmbers I a and Jb, and a term eqnal to the second term in Eq. (3.2) with an opposite 
sign which thereby is canceled ont leaving only connected terms in the resnlt for H. The same 
resnlt holds for the boson and fermion creation and annihilation operators independently of 
their commntation relations becanse each kind of the operators is commnted an even nnmber 
of times. 

After the second term in Eq. (3.2) is canceled one is left with a nnmber of partially 
contracted terms in which annihilation operators may still stand to the left of creation 
operators. A nnmber of ordering transpositions need to be done before a generic ordering 
of operators adopted in the previous Section is achieved. In fact, the process of commuting 
factors in A through factors in {B}q in the first term on the right-hand side of Eq. (3.2) 
produced above a number of terms with creation operators moved to the right of annihilation 
operators unnecessarily. These transpositions have to be undone to recover hnal asnwers with 
the adopted ordering. Nevertheless, it is visible that disconnected terms cannot appear and 
the following rule simplihes the calculations. 

The right-hand side of Eq. (3.2) equals the sum of the contracted terms which result 
from A{B}^ by moving n/=i 0 % through 11^=1 and, the contracted terms which result 
from —{B}fjA by moving (^wi through nf=i All other terms cancel out. 


3.b Yukawa theory 

The standard procedure from Ref. leads from the lagrangian density Cy = i’ii / 
d — m — g(l))'ip + — /i^0^) to the light-front hamiltonian expression of the form 


Hy = I dx d^x 


- , + rnJ 


+ -(j){-d 


±2 


+ /i ^)0 
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(3.3) 


-=o 


2 id+ 

We replace fields and 0(x) for = 0 by the Fourier superpositions of creation 

and annihilation operators, order the operators in all terms and drop the terms containing 
divergent integrals which result from the contractions. Then, we introduce the regularization 
factors. 

In the course of calculating effective hamiltonians we will also add new terms to Hy 
which ought to be included due to the presence of the regularization, in accord with the 
renormalization theory from the previous Section. For example, we will add a small term 
5rnl = ml — m? to m? in the hrst term and — /i^ to /i^ in the second term. We will 

calculate these additional terms below using the renormalization theory to order More 
additional terms can be calculated when terms with higher powers of the coupling constant 
g than second are included in the calculation of the renormalized effective hamiltonians. But 
the corresponding calculations go beyond the scope of the present illustration of the theory 
and will not be pursued here. 

In order to obtain the theory for particles with the quantum numbers of nucleons and 
pions one needs to include the isospin and replace the scalar coupling by * 75 . However, 


for the purpose of the illustration of the renormalization procedure to second order in the 
coupling g^ we do not have to introduce these factors explicitly. The additional factors merely 
lead to somewhat different algebra which can be traced throughout the whole calculation 
and hnal results including the isospin and *75 can be read from the results in the Yukawa 
theory. In this Section we assume m > /i > 0. 

The lagrangians with chiral symmetry would require massless fermions and additional 
meson helds, such as a. When the spontaneous symmetry breaking is considered in the 
light-front hamiltonian approach, one has to dehne the hamiltonian including the symmetry 
violating terms. A discussion of scalar theories with symmetry breaking, in the context of 
the hamiltonian renormalization theory, is given in Ref. 0 (see also Ref. ||I^). 


Meson mass squared 

The simplest example of a second order expression for a term in an effective hamiltonian 
in the Yukawa theory is provided by the meson mass squared. We hrst describe steps which 
produce this expression. The number of distinct steps in the procedure is 10: dehning the 
regularized initial hamiltonian, calculation of the effective hamiltonian, analysis of the cutoff 
dependence of hnite matrix elements of the calculated terms and extraction of the structure 
of the divergence, evaluation of the counterterm, isolation of the hnite part, calculation 
of the ehective hamiltonian knowing the structure of the counterterm, solving a physical 
problem such as an eigenvalue problem or a scattering problem, adjusting the hnite part of 
the counterterm to match data (including adjustements for the observed symmetries), and 
computing the hnal expression for the ehective hamiltonian with the hnite part determined 
from the ht to data. 

The simplest example is described in full detail of all the steps. Such extensive presen¬ 
tation is not provided in later examples where more complicated expressions would require 
too much space for such extended presentation. The hrst example is discussed in such detail 
despite the fact that in this case one can also proceed without the explicit discussion of 
all the steps and still obtain the correct result. Correct means here in agreement with the 
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similarity renormalization group for hamiltonians. 

For example, one might propose the structure of the counterterm using, as is usually 
done, some scattering amplitude instead of the matrix elements of an effective hamiltonian. 
Note that one can also impose renormalization condition using a scattering amplitude which 
results from a calculation performed without use of the effective hamiltonian. 

However, the systematic approach from Section 2 is the only tool we have for dealing 
with more complicated cases of light-front hamiltonians and their eigenvalue equations. In 
other words, the simplest available case is used to present all the steps in detail because 
it illustrates the procedure in a familiar setting. Then, we proceed to the expressions for 
fermion masses and interactions between fermions. Details are then discussed only in cases 
where a new feature appears in the calculation. 

Equation (2.36) implies to second order in Q 2 \ that 




Qi2X 


d_ 

d\ 


Pi4) 


Gix — E. 


-G21X 


IX 


11 


+ E 

p=3 


Qipx 


d 

dX 


Pizl) 


Gix — E- 


IX 


'Qp\x 


(3.4) 


11 


where the double-digit subscripts refer to the number of creation and annihilation operators 
(in that order) and the bracket subscript denotes the part which contributes to the rate 
of change of Qw with A. E^x is the eigenvalue of Qw which corresponds to the creation 
and annihilation operators indicated by the subscript 11. The reason for that only one free 
energy eigenvalue appears in the denominators is that Qw of Eq. (3.4) is a one-body operator 
and quantum numbers which label creation and annihilation operators in Qw are the same, 
including momentum. Therefore, the free energy eigenvalues are also the same: both are 
equal to E\\. Consequently, all commutators are written on the right-hand side of Eq. (3.4) 
in the simplihed form. The numerator similarity factors reduce to the single derivatives of 
/I (we have chosen n = 1 in Eq. (2.19)). Terms with more than two intermediate particles 
{p > 3) are of order or higher. 

Assuming that g in Eq. (3.3) is extremely small, writing Qw as a series in powers of g 
and keeping only terms order g"^, we obtain the following result from Eq. (3.4) for the meson 
free energy term. 




meson X 


k ^2 , „2 


E Pxj ^ 

(Xj.(Xk‘ 




(3.5) 


A remarkable feature in this result is that no correction arises to the term /k'^ which is 
protected by the kinematical symmetries of light-front frame. Simply, the total transverse 
momentum does not appear in a boost invariant expression. 

The width dependence of px is determined by the equation 


dpi 

dX 




dX - p^ 


r,(x, k), 


(3.6) 


where A4^ = (k^ + m^)/x(l — x). m? and p^ are the original bare mass squared parameters 
from Eq. (3.3). They do not include terms order g^ and higher because such terms would 
lead to higher order corrections than g'^ for the whole expression. The terms order g'^ and 
higher are treated as interactions in the perturbative calculation. 

In terms of graphs for the effective hamiltonian calculus, Eq. (3.6) represents the contri¬ 
bution of a fermion loop on a meson line. However, the graphs are not provided in order to 
avoid confusion with other diagrammatic techniques. 
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(3.7) 


Using Eq. (2.20) in the limit of a 6'-function, f{u) = 9{uo — u), one obtains 


fizl) = 9 


2 1 + 2 
A H-- /i 


1 

V ^'^0 


(3.8.a) 


For example, for uq = \ one has + ?>p — Ai'^]. Therefore, the derivative of 

f\ with respect to A forces the invariant mass of the fermion-antifermion pair, to be 
equal A^ -|- 3p. In general, the derivative selects the range of energies in the integral where 
the similarity function changes most rapidly. The regions where the function approaches 
a constant, i.e. 1 near the diagonal and 0 beyond the hamiltonian width, are strongly 
suppressed. The region that contributes is the edge of the diagonal proximum. The 
derivative of f\ is large and positive in this region and it approaches a 5-function in the limit 
of Eq. (3.8). 

In the limit of an inhnitesimally small Uq, as discussed below Eq. (2.20), one would 
substitute A^ = Then, 



P + p 




(3.8.6) 


The numerator factor in the square bracket in Eq. (3.6) originates from spinors of the 
intermediate fermions, Tr(^m + — fn) with Pm = Pm = The subscript m 

indicates that the — component is calculated from the mass-shell condition knowing -|- and 
T components. -|- and T components of p and p are constrained by the light-front spatial 
momentum conservation law, p + p = fc, where k is the meson momentum. The pseudoscalar 
interaction with zys gives the same result with an additional term -|-8m^ in the numerator. 

According to Eq. (2.49), 


r,(a;, k) 


1 + e 




+ 



-2 


x) 


(3.9) 


No infrared regularization is required in Yukawa theory with massive particles, m > 0 and 

/i > 0. 

We can see some useful flexibility in the regularization factors at this point. If the 
regularization factors in Eq. (2.49) contain e* divided by 1 — Xj in place of e* one obtains 
here 


r^(x, k ) 



(3.10) 


instead of Eq. (3.9). The spinor bracket on the right-hand side of Eq. (3.6) equals 8(p^)^ 
when one changes variables from x and to p with p^ = and p^ is determined from 
the relation AA. = 2pm'^ + p. Thus, with the modihed regularization factor of Eq. (3.10) 
which is a reasonable option, the spinor bracket can be replaced by = ^pAA."^ — rn?) and 
two separate integrations over x and can be reduced to a single integral over AA^. Such 
simplihcations are helpful in a qualitative analysis of the cutoff dependences. 

In the limit of Eq. (3.8.a) for Uq = \ one obtains 
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(3.11) 
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IGtt 



A2 + 2/i2 J 



dz 


[(1 + 2 / a )2 


2 ’ 


where a = e(A^ + 3iJ,‘^)/A‘^ and Zq = \J^ — dm^/(A^ + 3/i^). Note that for A^ < Am? — 3/i^ 
the derivative of the effective meson mass equals zero and the mass stays at the width 
independent value If one uses Eq. (3.10) instead, the corresponding result is 


Eqs. (3.11) and (3.12) are the same for a 1 which is the limit of removing the regularization 
cutoff, e —0, for a hxed hamiltonian width A. In this limit one has 


dA2 


a / /i^ A / 4m^ 

4^ V ^ A2 + 2/i2 ) “ A2 + 3/i2 


3/2 

^(A^ + 3/1^ 


4m^). 


(3.13) 


If one assumes that the meson mass squared parameter in the effective hamiltonian has 
some hnite value, /ig = /Uaq some Aq such that Aq > 4m^ — 3/i^ then, the integration of 
Eq. (3.13) demonstrates that 


2 ,.2 I /\2 \ 2 \ I ^ (,.2 C_ 2 ' ’ A 


U\ — hn -(A^ — An) ~\~ 


(/i - 6m ) log ^ + /ieon^(A, Ao) 


(3.14) 


Aconvi.?) Ao) denotes the result of integrating the convergent part of the integrand. 


/i 
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(A, Ao) 
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4m^ \ ^ 

s + 3/j? j 



(3.14.a) 


hLn?;(A, Ao) has a hnite (i.e. A-independent) limit for large A. It contains the terms which 
vanish for large A as inverse powers of A^. The dependence of /icont>(A, Aq) on m and /x is not 
indicated explicitly because we will not need it in the discussion of counterterms. However, 
one should keep in mind that the mass parameters determine the value of A = \/Am‘^ — 3/i^ 
where the effective mass stabilizes. We simplify our notation assuming that the effective 
cutoffs are always above the point of stabilization. Below the stabilization point, the meson 
mass has the constant value which is independent of A. We will show later that the constant 
value is equal to the physical meson mass. 

The mass squared term in the effective hamiltonian with a non-negligible coupling g 
grows linearly with A^. A logarithmic correction appears with an opposite sign due to the 
factor Zq in Eq. (3.13), as indicated in Eq. (3.14). However, one cannot make contact in Eq. 
(3.14) with the initial hamiltonian by letting A grow to inhnity because one would obtain 
an ill-dehned result. The factors depending on a in Eqs. (3.11) or (3.12) remove the inhnite 
growth of fix when A —> oo. Eq. (3.12) is simpler than Eq. (3.11) and illustrates the same 
idea so we start with Eq. (3.12). 

Equation (3.12) can be integrated over A from inhnity to any hnite value because the 
factor (1 + a)“^ provides convergence for A^ > A^/e. Nevertheless, the integral diverges as 
a function of e when e —> 0. The divergence appears as a single number. Therefore, the 
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counterterm is also a number. We add to jj? in the initial hamiltonian. We also 

write nl as a series in powers of g, jj^l = jj? + 5^1 + oi^g'^) so that (5/i^ ~ g"^. Thus, 




A 


a 


A2 487r 


1 + 




s + 2 g^ 



Am? 
s + S/i^ 


3/2 



s + 3/i^\ 

A2 j 


+ 0(//"). (3.15) 


This is an example of Eq. (2.33) in a perturbative application to second order in powers of 
g in Yukawa theory. The counterterm 5gl will be calculated following the steps described 
below Eq. (2.33). 

The diverging part of the integrand equals Q;/487r [1 + — 6m‘^)/s] and the remaining 

part is convergent. The convergent part of the integrand has the same structure as in 
glonvi^j ^o) hut the integral now extends from to inhnity instead of from Aq to A^. In 
the convergent part, one can replace the regulating factor by 1. Simplihcations occur in the 
limit e —0 and the result of integration in Eq. (3.15) is 

1 A^ 11 

—^7 + ) ( °®^A2 ^ ~6 

+ o(/)> (3-16) 

where the square bracket originates from the diverging part. 

Following the procedure described below Eq. (2.33), we dehne the counterterm dgl as 
the negative of the integral of the diverging integrand for some arbitrarily chosen A = Aq 
plus an unknown hnite piece corresponding to Aq and denoted by 5/ie/mite(-^o)- Namely, 


2 2 I r 2 I ^ 

Aa-A + ^ 



5gl = 


a 

487r 


A'l ,2 
+ a' 

3 e 


6m^) log- (Aq + 3/i^) — (/i^ — 6m^) 


7 Ag 11^ 


-\~ 5 /mite (Ao) + o((7^), (3.17) 

where 


finite ("^o) 


a 


487r 


^ci + {g^ - Qw?)c2 + (Aq + 3/i^) + {g^ 

O 


with the unknown finite numbers Ci and C 2 . So, in fact. 


6m2)(iog^ + ^) 


+ o{g^) 

(3.17.a) 


5gl = 


a 

487r 


A^(- + Cl) + {g^ - 6m^)(log - + C 2 ) 


(3.17.5) 


Since the whole expression on the right-hand side of Eq. (3.17) is merely a number, it is 
not necessary to hnd Ci and C 2 or any other part of it separately at this point. One can easily 
hnd the number 5/ie/inite(-^o) from the knowledge of g\ at some value of A. The resulting 
counterterm will render well defined finite boson mass squared parameter in the effective 
hamiltonians in the limit e —0. Using Eqs. (3.16) and (3.17) one obtains 
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oi ' A^' 

^1 = + ^ - Ag + (/i^ _ 6m2) log ^ - ^il^^{oo,\) + o{g^). (3.18) 

Equation (3.18) is an example of Eq. (2.35). 

The unknown finite term, 5jjil, has to be found, in principle, by comparison with 
data. We shall discuss an example of a renormalization condition later in this Section. 

Let us assume now that at some arbitrarily chosen value of A = Ai the effective meson 
mass squared required in Eq. (3.5) by a £t to data equals /i^, i.e. We can calculate 

^finitei^o) usiug Eq. (3.18) with A = Ai and /i|^ on the left-hand side replaced by the 
number inferred from the experimental data. The result is 

a \ A^l 

X\-Xl + {g^-Qw?)\og^ +/i2^„^(oo,Ai)+o(/). (3.19) 

Note that one has to include the contribution of the convergent terms in the determination of 
the arbitrary constants when using the value of /i^. Knowing 5/ie/inite(-^o) one can calculate 
jji\. Namely, 

/^A =/^i + A^ — A^-h (/i^ — 6 m^) log ^ -^i) + o(fi'"^)- (3.20) 

As expected, Eq. (3.20) is the same as Eq. (3.14) when Ai = Aq and Hi = Hq. One can also 
trace the origin of all the terms from Eq. (3.14); the diverging and converging terms and 
the counterterm in Eq. (3.15). 

Now, we can approach Eq. (3.11) analogously to Eq. (3.12) without calculating all 
integrals explicitly. Integrating both sides of Eq. (3.11), we have 

ds . (3.21) 

By demanding that removes the diverging (i.e. e-dependent) part of the integral, and 
leaving the finite part of /i^ free so that at some A = Aq the effective boson mass squared 
parameter has a desired value /Tq, we obtain 

/^A = ■ (3.22) 

The integrand is given by the right-hand side of Eq. (3.11) with A^ = s. Since s ranges only 
from Aq to A^ and both are finite we can take the limit e —> 0 under the integral sign and 
the integrand becomes equal to the right-hand side of Eq. (3.13) with A^ replaced with s. 
Integration over s produces exactly the same answer as given by Eq. (3.14). Thus, we see 
that the result of Eq. (3.14) is independent of the regularization scheme. It is determined 
by the initial hamiltonian of Yukawa theory as given by Eq. (3.3). The only unknown in 
Eq. (3.14) is the value of /ig. More precisely, we know that /Xq = -|- (5/iQ -|- 0 ( 5 ^^) and the 

unkonwn term is (5/iQ ~ 5 ^^. 

Note that the above calculations can be carried out in a different way using the following 
observation. Equation (3.6) in the lowest order of perturbation theory has a particulary 
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simple structure. Namely, the only dependence on A comes from the similarity function flz^) 
and both sides of the equation are equal to the derivatives with respect to A. Therefore, 
one can directly integrate both sides and obtain a compact integral expression for for 
arbitrary functions /. 

One should remember that such simplihcations do not occur in higher order calculations 
or beyond simple perturbative expansion when the coupling constant and mass parameters 
depend on A themselves. Therefore, we stress that the dominant contribution to the rate of 
change of /i^ with A comes from the edge of the diagonal proximum. This fact remains gen¬ 
erally valid and the procedure applied above represents a generic situation despite simplicity 
of the example. This example does not involve a distinction between the bare coupling g 
and a renormalized coupling because to order there is none. 

In order to determine 5/ig we need to specify a renormalization condition. An example 
of a renormalization condition which determines /Xq will be described now. 

A natural requirement for /Iq is that the effective hamiltonian has one boson eigenstates 
parametrized by momenta p'^ and p-^ with eigenvalues of the form p~ = + where 

fi is equal to the physical boson mass. Our approach preserves kinematical symmetries of 
the light-front frame explicitly and the eigenvalue is guaranteed to appear in that form. 
Therefore, one can calculate a whole spectrum of eigenvalues for eigenstates of different 
momenta by calculating the single mass parameter fl. 

In order to write the effective eigenvalue equation and hnd out /ig which leads to the 
desired value of /2 (if it is possible), the following steps need to be taken. 

One inserts Eq. (3.14) into Eq. (3.5). Then, one replaces the bare operators, al and a^,, 
in the expression for Qx by the effective ones, and axk, to obtain Gx- Next, by applying 
the operation Fx, one obtains the effective hamiltonian with the form factors in the vertices, 
Hx = Fa[Ga]. 

The effective eigenvalue equation for bosons is an example of Eq. (1.2). Here, it is 
sufficient to consider the eigenvalue equation for Hx in the expansion to second order in g 
since we are interested in which is proportional to g^. The only terms that contribute are 
the free energy term including the effective mass squared and the interactions which change 
the particle number by one. The latter equal the canonical interactions with the similarity 
form factors. 

To zeroth order in p, a physical meson state equals a single effective meson state, and 


No terms order g arise in and the next correction is order g^. This correction can 
be calculated using the operation R and the model hamiltonian dehned in a perturbative 
expansion from Eq. (1.3). Using expansion into a series of powers of g to second order, one 
can restrict the model space to the single effective boson sector. The effect of coupling to 
the fermion-antifermion pair states is included in perturbation theory. 

Another method is to arbitrarily limit the number of effective Fock sectors and diagonalize 
the effective hamiltonian in that limited space. Such procedure could be called the effective 
Tamm-Dancoff approach (ETD), c.f. Refs. Q and [0. The term of the second order in 


g in the eigenvalue will determine dp^. In this calculation one can limit the space of states 
to the one effective boson and the effective fermion-antifermion pair. Note that in the ETD 
approach no ultraviolet renormalization problem emerges. The interactions are tempered by 
the similarity factors of width A which is on the order of particle masses. This way the ETD 
approach avoids the old problem of renormalization in the TD approach. 
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The model calculation using the operation R and the ETD calculation have to agree with 
each other for small coupling constants. The reason is that corrections to the zeroth order 
approximation tend to zero when 0 in the presence of a finite energy gap between states 
with different numbers of particles, which is the case here for 0 < /i < 2 m. 

In this paper, we discuss only the second order corrections in g and a very limited number 
of the effective Fock sectors can count (i.e. 2). Therefore, we can focus on a straight-forward 
perturbation theory anyways. Nevertheless, our simple calculations have two interesting 
aspects. 

The hrst one is that no coupling renormalization effects arise to order Therefore, the 
expansion in powers of g up to g^ is equivalent to the expansion in powers of a well dehned 
effective coupling, no matter how the latter is dehned. Thus, it is important to realize that 
the expansion in powers of g is understood here to be the expansion in powers of the effective 
coupling which appears in the effective hamiltonian of width X] gx = g to order g‘^. It is not 
meant to be the expansion in the initial held theory coupling constant. 

The second aspect is following. The perturbative expansions applied in the ehective 
eigenvalue problem are expansions in the interaction which is suppressed in strength by the 
vertex form factor of width A, provided by the similarity function fx- Therefore, the actual 
range in momentum space of the ehective interactions proportional to g in the ehective 
hamiltonian is inhnitely smaller than the momentum range of the analogous interaction in 
the bare hamiltonian. In other words, the ehective strength of the interactions is greatly 
reduced and much smaller than the value of g itself would imply if it stood in the initial 
bare hamiltonian. Our initial expansion in powers of g can now be understood as a shortcut 
to the expansion in powers of the ehective coupling. The latter expansion may have a 
considerable range of rapid convergence because the form factors reduce the size of coefficients 
in the expansion. The ehective coefficients are much smaller than in the case of the initial 
hamiltonian without form factors. 

Thus, the operation R on Hx expanded in powers of the ehective interaction (the cou¬ 
pling constant itself can be sizable), opens new options for studying the ehective eigenvalue 
problem in the whole Fock space using the basis built with the ehective creation and anni¬ 
hilation operators corresponding to the width A. One can estimate contributions of various 
components by making diherent choices of the model spaces and solving model dynamical 
problems numerically. Wave functions are expected to fall oh sharply with growing momenta 
and particle number if g is not too large for the vertices to contain sufficiently narrow simi¬ 
larity factors so that the particle number changing interactions are small perturbations. In 
that case, the result should still match the outcome of the hrst two approaches when all 
quantities are dehned recursively in powers of g. The coupling needs to be set equal to the 
right value. 

The second order expression in perturbation theory implies 




p'' 


p-t 




Fx[G2 ix]\p). (3.23) 


Ip) denotes a single ehective meson state with momentum p+ and p^, {p'\p) = 167r^p'^S^{p' — 
p). The term order g‘^ produces 
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= /^A - / 9f{4) + (3.24) 

where the notation is the same as in Eq. (3.6). Using Eq. (3.8) with mq = | at A = Aq one 
obtains 


For Aq < \/4:m‘^ — 3/i^ the effective meson mass parameter equals the physical meson 
mass, as promised. For fermions with masses order 0.9 GeV, this implies no corrections to 
a light meson mass, such as tt, for cutoffs smaller than 1.8 GeV. But one has to remember 
that the actual expression for the correction is different for the pseudoscalar ttN interaction 
than in Eq. (3.25), i.e. the spin factor has to be enlarged by 8 m^ (see comments above Eq. 
(3.9)). The size of the coupling constant a may invalidate the perturbative result for pions 
coupled to nucleons because the coupling constant on the order of 10 may appear when one 
attempts to make comparisons with data. 

However, the actual measure of the off-shell effects is not given directly by /ig but by the 
sum of fj,Q and the self-energy resulting from the effective interactions. According to Eqs. 
(3.22) and (3.24), the sum of both contributions in the physical pion mass itself is actually 
equal zero to order a . It is not clear at this stage of the development of the effective 
hamiltonian theory if the off-shell corrections of physical interest must be large due to the 
coupling constant order 10. 

Using Eqs. (3.21), (3.22) and (3.25), one can express the meson mass squared term in 
the initial renormalized hamiltonian in terms of the physical meson mass fl and the initial 
mass parameter jj,. Namely, jj? = and 

= /^^ + / / «) + o{g*). (3.26) 


Fermion mass squared 

In complete analogy to Eqs. (3.4) to (3.6), one obtains the fermion energy operator, 

/ _L2 2 

= y/ lk]—^{bLh. + did,,). (3,27) 

Results for fermions and antifermions are identical and we will describe explicitly only 
fermions. We have 


dm\ 

TT 


= g^ 


dp{zl)Umak{^m + m)Umak , x 
XK\ - 77 -^- r^{x, k). 


d\ 






= (m^ -|- K^)/x + -|- K^)/(l — x). The regularization factor of Eq. 

intermediate particles implies 


r^{x, k) = 


1 + ^A1' + 
A2 


e + m? + g, 


AV 


X 


1 — X 


-2 


(3.28) 
(2.49) for the 

(3.29) 
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The spin factor in Eq. (3.28) can be rewritten as 

T '^mcrk\p^ T ^ T iPm (3.30) 

l^rn between spinors is equivalent to m. The term with q"*" is typical in light-front calculations. 
The term in it, proportional to cancels out. The term linear in does not contribute 

because it is odd in and all other factors (especially the regularization factor) depend 
only on the modulus of k^. Regularizations which are not even in would produce a term 
linear in k^. In the present case, the spin factor is reduced to 


'^mak 


[X 


l)m + 


+ {I — x^)rn^ 


2 k+ 


X 


'^mak — 


+ {1 + X^W? 


X 


(3.31) 


The analogous spin factor for a pseudoscalar interaction with *75 is the same as the right- 
hand side of Eq. (3.31) except for the opposite sign in front of x in the numerator. Inclusion 
of the isospin introduces the number of bosons in the theory in front of the integral in Eq. 
(3.28). 

We observe a similar structure in Eq. (3.28) as in the meson mass dependence on 
in Eq. (3.6). Namely, there are terms diverging linearly and logarithmically and there is a 
series of convergent terms. The integrals are slightly more complicated because two different 
masses determine the integrand dependence on x. Otherwise, the procedure is essentially 
the same. 

Instead of evaluating the integrals explicitly we observe that the divergences amount to 
a number which grows when e —0. We integrate both sides of Eq. (3.28) and obtain 


ml = ml- 




+ {1 + xYm? 
x{M.‘^ — m?) 


r,(a;,K), 


(3.32) 


where, according to Eq. (2.18), = (AI^ — + w? + A^). The e-dependent terms 

originate from the first term in the bracket. It is independent of A. The counterterm ~ 
in nil = + Sm"^ + olg"^) removes the divergence. The finite part of the counterterm is left 

to be determined by data. Thus, one can write 


ml = m^ + 5ml 


+ g 



k'^ + {1 + xYm^ 
x(A4^ — m?) 


(3.33) 


and one can determine the value of from the value of m\ which is required by 

some physical condition to be satisfied by the effective hamiltonian at some value of A. If 
we had defined the divergent part by an integral from some value of A, such as Aq in the 
meson case in the previous example, then we would have to take into account that 
depended on Aq to compensate for the Aq dependence of the diverging integral. When we 
define the counterterm to be given by the whole A-independent part of the integral in Eq. 
(3.32), plus a finite constant to be determined by data, then does not depend on 

Aq. Nevertheless, it can be expressed in terms of For example, if for some reason the 
effective fermion mass squared at A = Aq should be = ttIq then, one obtains 


2 2,2 

mx = mQ + g 





+ {1 + xYm^ 

x(M^ — m^) 


(3.34) 
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ttiq = m? + 5m^ + o{g'^) and ~ can be fonnd from a renormalization condition for the 
physical fermion mass. 

A natnral condition for htting is that the effective hamiltonian at the scale Aq has 
fermionic eigenstates with eigenvalnes of the formp“ = + where fh denotes the 

physical fermion mass. In analogy to Eqs. (3.23) and (3.24), one obtains 


So, 


~ 2 2 
m = rrix 



+ (1 + xYrri^ 

— m?) 


9f{zl) + o{g^). 


(3.35) 


tUq = m 




(3.36) 


The initial can be calcnlated in terms of m^, m^, g and e from Eq. (3.32). The effective 
fermion mass parameter in the interacting hamiltonian of width A is 


+ g' 


' J[xK]f 


K‘ 


+ (1 + xYm? 


+ o{g% 


(3.37) 


x{M.^ — w?) 

Analogons eqnation in the case of nncleons conpled to pions is the same as Eq. (3.37) 
except for (1 — x)^ instead of (1 + x)^ in the nnmerator and the isospin factor 3 in front of 
the integral. In the limit of a 6'-fnnction for the similarity factor / one obtains the resnlt 


m 


= mr + 0 ( 5 '^), m = niN, and 


+ ^9^ j + 3m^ — M^) 


+ {1 — xYm? 

x(A4^ — m^) 


+ o(g^ 


(3.38) 


For A^ = (m + n/u.^)^ — 3m^, where n is a nnmber order 1 or a few, one obtains (a = g^/47r) 


n n n 3Qi 


(3.39) 


where c = 4:/3{nfin/mNY + o(/r^). The expansion formnla for c shows how small the 
correction is for small meson masses. Note that A^ mnst be negative for small n, as explained 
below Eq. (2.20). The exact resnlt for c is that, for example, n = 3 gives c ~ 0.03 and n = 4 
gives c ~ 0.12. Thns, even for qnite large conplings the effective mass parameter in the 
hamiltonian deviates only a little from the physical nncleon mass if the hamiltonian width 
allows momentnm changes of the order of a few meson masses only. Moreover, the physically 
relevant off-shell effects are not given directly by these nnmbers bnt by the difference between 
these and the effects of the interactions present in the effective hamiltonian. The combined 
effect is zero for the nncleon mass itself to order g^. 

Note that, if we nsed Eq. (2.20) in the ^-fnnction limit with an ihnitesimal uq and 
A^ = A^/ydfo for n = 1 then, the 6'-fnnction nnder the integral in Eq. (3.38) wonld be 
replaced by ^(A^ -|- — Ad^). 

The present calcnlation does not say what is the size of the off-shell effects in the theory of 
pions conpled to nncleons bnt Eq. (3.39) does not exclnde that the effects can be calcnlable 
in pertnrbation theory. This is enconraging for the program ontlined in Ref. . 


Fermion-fermion interaction 
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Our next example is the second order calculation of the effective hamiltonian term which 
contains products of two creation and two annihilation operators for fermions. This term is 
particulary interesting because it complements the one boson-exchange interactions which 
result from emission and absorption of the effective bosons by the effective fermions. 

The differential eqnations we need to consider are given by Eq. (2.29) for the two-fermion 
terms. We have 


d\ 


Q22\ — 


fxGl 2 \ I “ /a)^ 21 A 


GlX 


' d 
dX 


(1 — fx)Gl 2 X \ fxG 2 lX 


Gix 


(3.40) 


J 22 


The snbscript 22 denotes a term with two creation and two annihilation operators for 
fermions. 21 denotes a term with one annihilation operator and one creation operator for 
fermions and one creation operator for bosons. 12 denotes a term which annihilates a fermion 
and a boson and creates a fermion. For a hermitean hamiltonian, we have G 12 = ^ 21 - 
higher order terms, the subscript notation has to be generalized to represent the meaning 
of the snbscripts more precisely and, mnltiple snbscripts are needed. It is not necessary to 
introduce the more general notation in second order calculations in this paper. 

The right-hand side of Eq. (3.40) does not contain disconnected interactions (it never 
does) and one can isolate the terms with two creation and two annihilation operators for 
fermions by contracting one creation operator and one annihilation operator for bosons. 

The only term which contribntes to the right-hand side of Eq. (3.40) in Yukawa theory, 
is the ordered and regnlarized third term on the right-hand side of Eq. (3.3). Thus, in Eq. 
(3.40), we have 


^2iA = H j[kikk2]lQTi‘^5{ki + k- k2) gu^k^„jU^k2ur{eef /K^)r{eeb/JSG)h\^^^a\hk^„ 

/ [x^A9UmxP+KaJU^Par{eef /r{eei,/A^) , (3.41) 


1 

T+ 


where e/ = {p? -|- wA)/x and e?, = [p? + h^)/(l ~ ^)- This representation illustrates appear¬ 
ance of the parent and danghter momentnm variables in the interaction term. The factor 
gUmkiafUmk 2 <Tr{eef/hG)r{eeh /will be denoted by g 2 i- The analogous factor in Gi 2 X 
will be denoted by 5 ^ 12 . 

Similarly, in the case of the four-fermion interaction term, we have 


G22X= j[P]^ Y. /[a^'^]bP]^22A&Ip+«,i&(l_^)P_«,3&j,P+P-2^(l-J/)P-pa4, (3.42) 

CTiCT2f73Cr4 

where g 22 x is a fnnction of the daughter momentum variables x, p^, y, and the fermion 
spin projections on the z-axis Ui, (J 2 , <73 and < 74 . Details of the notation will become clear 
shortly. 

To order g"^, only fx depends on A on the right-hand side of Eq. (3.40). Eqnation (3.40) 
can be rewritten in a simplihed fashion in terms of the coefficient fnnctions snch as in Eq. 
(3.42). Namely, 
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( 3 . 43 ) 


^fl'22 — [fi'l25'2l]22 • 

The subscript A and arguments of the functions are not indicated, to simplify the formula. 
Expression in the hrst bracket is called the inner similarity factor ioi ^ 22 - The word “inner” 
is used to distinguish this factor from the form factor introduced by the operation Fx when 
this operation is applied to G 22 - The latter form factor can be called the outer similarity 
factor because it depends on the incoming and outgoing invariant masses only, independently 
of the internal structure of the interaction. 

Momentum variables in Eq. (3.43) can be expressed by the daughter momentum variables 
from Eq. (3.42). One needs to express the parent and daughter variables of Q 12 and ^ 21 , 
and the energy denominators, in terms of x, y and from Eq. (3.42). This is done as 
follows. 

The momentum labels of the fermion annihilation operators in Eq. (3.42) are denoted 
by k 2 and A ;4 and the momentum labels of the fermion creation operators in Eq. (3.42) are 
denoted by ki and ks . The numbers assigned to the fermion momenta correspond to the 
numbers labeling their spin projections on the z-axis in Eq. (3.42). In terms of the origin 
of the annihilation and creation operators in Eq. (3.42), Q 21 provides the fermion creation 
operator with momentum ki and the fermion annihilation operator with momentum k 2 - Gu 
provides the fermion creation operator with momentum k^ and the fermion annihilation 
operator with momentum k^. There is a change of sign due to the reordering of the fermion 
operators. The boson operators from Q 12 and Q 21 are contracted and provide factors similar 
to the factors associated with an intermediate particle in the old-fashioned hamiltonian 
calculations for scattering processes. 

It is useful to think about the effective hamiltonian in terms of a scattering amplitude 
with two vertices but the reader should remember that the formula we are describing is 
not for an S'-matrix matrix element. Therefore, the “scattering” language has a limited 
meaning. The hne point is that, after evaluation of Qx, one has to replace the bare creation 
and annihilation operators by the effective ones in order to obtain Gx- The term in Qx is not 
directly related to any scattering process before the replacement is made. The replacement 
prevents confusion between the hamiltonian calculus which uses the bare operators, and 
the S'-matrix calculus which uses the effective hamiltonian and the corresponding incoming, 
outgoing and intermediate states of effective particles. The scattering language becomes 
particulary confusing in theories with gauge symmetry, spontaneous symmetry breaking and 
conhnement. None of these features appear here in the calculation to order Therefore, 
the scattering language is useful in the current example. 

The intermediate boson momentum is dehned for -|- and T components as k^ = k^ — k^ 
and kf^ = {kf‘^ + pi^)/kf. These four components form the four-momentum of the exchanged 
boson. It is denoted by k^^ to indicate the mass which determines the minus component. 
The same result for k^^ is obtained by subtracting ki from /c 2 instead of ^4 from k^,. It is so 
because the translational invariance of the hamiltonian on the light-front implies momentum 
conservation for the -|- and T components. 

Thus, the inner similarity factor in Eq. (3.43) is 


[/{-n - {-/'}/] 



[-/(4a)]^ 

AE'21 


[-/(4a)]^ 

AEi2 



(3.44) 
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The prime denotes differentiation with respect to A. 

In the case of Qi in Eq. (3.4), the whole inner similarity factor of the analogous structure 
was equal to the derivative of Piz^) divided by a single denominator. For both functions / 
in Eq. (3.4) had the same argument z'l and the two corresponding energy denominators were 
the same. In Eq. (3.44), we have two different arguments of the similarity functions / and 
two different energy changes. Namely, Zi 2 x and AE 12 for the vertex of Q 12 with momenta 
ksm, kim and /cs^, and Z 2 i\ and AE 21 for the vertex of Q 21 with momenta /cim, /^ 2 m and k^^. 

The parent momentum for the vertex of Q 12 is P \2 = {k^ii + k^m + k^PP so that for the 
+ and T components we have P 12 = k^. Similarly, the parent momentum for the vertex of 
Q21 is P21 = {k^fj, + kim + k2m)P so that for the + and T components we have P 21 = k2- 
Now, the rules provided by Eqs. (2.12) to (2.20) imply the following formulae for the 
arguments of the similarity functions /. 


AAI^2 ~ + ^4m)^ ~ + k^rn ~ ^3m)-Pl2- 

(3.45) 

SAI 12 = -^12 + 

(3.46) 

AMI 2 

(3.47) 

AAI 21 = ^2m ~ {k^^^ + kirnP = + kim — k2m)P2l- 

(3.48) 

SAI 21 = “-^21 + 

(3.49) 

AMl^ 

“ SAIii + A2' 

(3.50) 

Equations (2.12) to (2.16) imply 


AF 

AE 12 - + 

-'12 

(3.51) 

and 


p+ 

-'^21 

(3.52) 

In terms of the familiar parameters x, y and from Eq. (3.42), i.e. 


P = ki + k‘i = k2 + ki, 

(3.53) 

X = kf/P+, 

(3.54) 

P = k^ — xP-^, 

(3.55) 


and 
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y = kt/P^, 


(3.56) 


= k^ - yP^ , (3.57) 

the mass differences which determine the arguments of the similarity functions read as fol¬ 
lows. 


AMI 2 = {1 — X 


— p-*-)^ + ^ p-*-^ -I- -I- 


AMl^ = y 


y — X ^ — y I — X 

p^^ -I- -I- {k^ — p-*-)^ -I- 


y X 

The corresponding energy denominators are 


y-x 


AE 12 — 


^ 


y-x 


+ 


1-1/ 


1 — X 


/P^ 


and 


(3.58) 

(3.59) 


(3.60) 


rp^^ -I- + w? {k-^ — p^Y + 

ZA^21 “ - — - — - 

\_ y X y — X 

In Eqs. (3.40) to (3.61), always y > x. When evaluating matrix elements of the effective 
interaction between states of indistinguishable fermions, one obtains results in which the 
momentum and spin variables are properly symmetrized (antisymmetrized) as dictated by 
the statistics. 

Evaluation of the second bracket in Eq. (3.43) gives 


/P+. (3.61) 


[//l2P2l]22 9^m{l—x)P—Ka3^m{l—y)P—piJ4^Y(^4/A )T’(6ei2/A ) 

1 




{y - x)P+ 

The arguments of the regularization factors appear in the mass differences. Namely, 


(3.62) 


where 


AAI?2 — 64 + C12 


= 


nfi + + y^ 


m 


X 12 


X 12 


X 12 


1 - 1 / 

1 — x’ 


(3.63.a) 


(3.63.5) 


kY = -P^ + Xi2K^, 


and 


where 


AAI21 = — ^1 + 621 = rn? 


K 21 + rri^ 
X21 


+ y^ 

1-2:21 


(3.63.C) 


(3.64.a) 
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(3.64.6) 


X 

X21 = 

y 


/^ 2 i — ^ ^ 2 iP • (3.64.c) 

The daughter energies defined by Eqs. (2.48) and (2.49) appear in Eqs. (3.62) to (3.64). 
The energy subscripts will be used to distinguish the corresponding regularization factors 
in an abbreviated notation below. Equation (3.43) reads, in the abbreviated notation, as 
follows. 


^ ^ [ y f 12/21 

y — X AAfii 


1-a: f[^j2i 
y-x AMf 2 


y(y - x) gusu^ r4ri2 guiU 2 rir 2 i. 


(3.65) 


In the familiar limit of Eq. (3.8) where the similarity function / approaches the 6'-function, 
be. fi 2 = 0[\^ + 2w? — AAI 12 ] /21 = 6[\^ + 2w? + AA^ 2 i]) the derivatives of the 
similarity functions become ^-functions and one can integrate Eq. (3.65) using the relation 
/“ ds6{s — b)6{s — c) = 9(c — b)6{c — a). The result is 


fl'22A — g22e + 


1/6'21-126'21-A _|_ (b ~ a;)6*i2-A6'i2-2i 


{y-x)\AMli\ {y-x)\AM\2\ 

The symbols 9 with various subscripts denote the following functions: 

^12-21 = 1 — 6*21-12 = 6 *(|AAf 42 l ~ |b^Al2il)) 


9{y - x) gu3U4 uru guiU 2 rir 2 i. 

(3.66.a) 


(3.66.6) 


and 


6*i2-a = d{\AMl2\ - - 2m^) 


92i-x = 9{\AMl3\-\^-2m^) 


(3.66.C) 


(3.66.d) 


The initial value term of g 22 e at A = cx) is absent in the canonical hamiltonian. It should 
be set to 0 unless it is found that matrix elements of H 22 \ = F\ [G 22 A] between finite free 
energy states depend on e when e —0 and a counterterm containing nonzero g 22 e is required 
to remove this dependence. 

Note again that, if we used Eq. (2.20) in the ^-function limit with an ifinitesimal Uq and 
A^ = A^/ydlo for n = 1 then, A^ + 2m^ in the 6*-function arguments above would be replaced 
by A^. This feature will be used later in the case of nonrelativistic bound states. 

The easiest momentum configuration to analyse is the one where the sum of free energies 
for the momenta of creation operators equals the sum of free energies for the momenta of 
annihilation operators. In this case, we have (fcim + ^sm)^ = {k2m + kirnY — W^e will refer 
to this configuration as the energy-diagonal part of the interaction. In the energy-diagonal 
part of the interaction, we have AE 12 = — AE 21 , and 


-I- rn^ 

1 -y 


-I- mf 


X 


-\- mf p-*-^ -I- 


X 


(3.67) 
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Thus, 


1 — X y y — X ’’ 

where = {k-^ — + {y — x)"^. These relations imply 012-21 = 1 —02i-i2 = 0(1—a; — ?/), 

012 -a = 0 [/i^ + g^ - {y-x){\^ + 2w?)/{l-x)] and 02 i-a = 0 [/i^ + g^ - [y - x){\^ + 2w?)/y]. 
Therefore, 0 i2-a = 1 in the same momentum range where 02 i-a = 1- Thus, the energy- 
diagonal part of the fermion-fermion effective interaction order is 


022A — 022e T 


gU'iUi r4ri2 guiU2 rir2i 


+ q 


,72 


9{y — x)9 


2 I -*2 

p +q 


X 


max{y, 1 — x) 


(A^ + 2m^) 


. (3.69) 


Equation (3.69) is helpful because it provides insight into the more complicated interac¬ 
tion from Eq. (3.66). When the momentum transfer is sufficiently large and {y — a;)(A^ -|- 
2rri^)/max{y, 1 — x) is negligible so that the 0 -functions and the regularization functions r 4 , 
ri 2 , ri and r 2 i in Eq. (3.69) equal 1, the second term on the right-hand side of Eq. (3.69) 
is equal to Feynman’s expression for the one boson-exchange scattering amplitude for two 
fermions. Namely, the numerator factors are standard for the Yukawa interaction and the 
denominator equals /i^ —(^ 3 ^ —/^ 4 m)^ = — {,k 2 m — kimY. The necessary antisymmetrization 
for identical fermions requires evaluation of the matrix elements of the hamiltonian term. 
However, the 0-functions and the regularization factors in Eq. (3.69) produce a difference be¬ 
tween the energy diagonal part of the effective hamiltonian matrix elements and the on-shell 
Feynman scattering amplitude. 

The 0-functions force the momentum transfer carried by the intermediate boson, |^, to 
be larger than [{y — a;)(A^ -|- 2m?)/y — p^] ' . The size of this number depends on how large 
A^ and the ratios of the longitudinal momenta are. If A^ is negative and compensates 2m^, 
the lower bound on the momentum transfer is absent. For larger A^, the ratio oi y — x, i.e. 
the intermediate meson longitudinal momentum fraction, to the parent fermion momentum 
fraction y has to be smaller than p^{\^ + 2m?)~^ for the lower bound on the momentum 
transfer to be absent. Otherwise, the momentum transfer is limited from below. This means 
that the effective interaction term does not include the long distance part of the Yukawa 
potential. 

The regularization factors in the limit e —0 converge pointwise to 1. We shall explain 
how it happens in the more general case below. No diverging cutoff dependence is obtained 
when evaluating matrix elements of G 22 A between states of hnite invariant masses so 
that the matrix elements of H 22 \ are also free from divergences. Therefore, g 22 e = 0. We can 
replace the regularization factors in the limit e —0 by 1 . 

We proceed to the analysis of Eq. (3.66). No divergences appear in the hnite matrix 
elements of the effective hamiltonian of width A when e —0. One can see that this is the 
case using Eqs. (3.63.a) and (3.64.a). Namely, the arguments of the regularization factors 
are hnite for hnite AAI 42 and they approach 0 when e —0. One demands that 

the free invariant masses of the states of fermions used to calculate the matrix elements are 
hnite. Then, the only possibility for AA 1^2 or AAf^i fo diverge emerges when x approaches 
y, i.e. when the longitudinal momentum transfer between the fermions approaches zero. In 
such case, ei 2 and 621 approach inhnity even for a vanishing transverse momentum transfer 
because the meson mass squared is greater than zero. 
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Now, the remaining factors of spinors and energy denominators, the latter multiplied by 
the boson phase-space factor oi y — x, are all finite in the limit x ^ y. The regularization 
factors ri 2 and r 2 i deviate from 1 only in the small region in the momentum space where 
\x — y\ < eh ?/(or in a still smaller region for a nonzero meson transverse momentum). 
All other factors in the interaction are hnite in this region. Therefore, for hnite wave packets 
used in the evaluation of the matrix element, this small region produces a contribution which 
is proportional to e. Thus, it vanishes in the limit e —0. Consequently, the matrix elements 
of g 22 e are equal 0 and the regularization factors can be replaced by 1 . 

The full result for the effective fermion-fermion interaction in the limit e —0 can now 
be rewritten as 



where 


H22X — Fx [G22A. 


E / [yp] ^ 22 A /(4a) bi^P+naA{l-x)P-Ka^bxyP+pa2h{l-y)P-pa^ , (3.70.a) 

0'lC2Cr3(74 '' 


fl'22A — 


^12-21^21-A 
- (I 2 I 


+ 


^12-A^21-12 

- Qi2 


9{y - x)gu3UiguiU2, 


(3.70.6) 


gi2 — k^m 


-k 


Ami 


(3.70.C) 


(I 2 I — k2m 


-k 


lmi 


(3.70.d) 


^12-21 — 
9i2-X 


1 — 6 ^ 21-12 




x){y^ - 


- x){y^ - ql^) - y{y^ - qli)\ , 

4 ) - {y - x){\^ + 2m^) \ , 


(3.70.e) 

(3.70./) 


6^21-A = 0 


y{p^-ql,)-{y-x){\^ + 2m^) 


(3.70.(7) 


The argument of the outer similarity factor /(4a) which limits the width of the effective 
interaction in momentum space, is given by 




(3.70.6.) 


The mass difference, AM 22 = Ml, - Mfs, and the mass sum, SMI 2 = Ml, + Mfs, can 
be expressed by the fermion momenta using relations Af 24 = ( 62 m + k,^)^ and Mf^ = 
(klm T k^Yn) • 

Equations (3.70.a) to (3.70.h) explain the structure of the fermion-fermion effective in¬ 
teraction order g"^ in terms of the two four-momentum transfers, qi 2 and ^ 21 - The transfer 
qi 2 appears in the vertex where the intermediate boson is annihilated and the transfer ^21 
appears in the vertex where the boson is created. In the energy-diagonal part the two 
four-momentum transfers are equal. In the general case they are different. The 6 -functions 
exclusively select which momentum transfer appears in the denominator. The lower bounds 
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on the momentum transfers depend on the ratio oi\x — y\ to y and 1 — x and on the masses 
and A^. 

We can now evaluate matrix elements of the T-matrix between the effective two-fermion 
states using the effective hamiltonian of width A to second order in g. We use the formula 

T(E) = Hn + ^ . ( 3 . 71 ) 

We have Hqx = Gix and Hjx = Fx [G 22 A + Gi 2 A + G 21 A] • The first term on the right-hand side 
of Eq. (3.71) contributes solely through H 22 x- In the second term, in Hjx, only Hi 2 x + H 21 X 
contributes. 

The first term in Eq. (3.71) has its matrix element given by the antisymmetrization 
of the right-hand side of Eq. (3.70.b). The multiplication by f{z 22 x) does not matter 
because f{z 22 x) = 1 in the energy-diagonal matrix elements. Only the energy-diagonal part 
contributes to the cross section. The energy-diagonal part of g 22 x is given by 


fi'22A — 


gu3U4 guiU2 
yt? — 


0{y -x)ex, 


(3.72) 


where 6x = 0 [max{y, 1 — x)^^? — q^) — {y — x){\‘^ + 2m^)] and q = qi 2 = q 2 i- The antisym¬ 
metrization of the right-hand side of Eq. (3.72) produces the contribution of the hrst term 
in Eq. (3.71) to the scattering amplitude. 

The second term in Eq. (3.71) provides the one-boson exchange amplitude with form 
factors in the fermion-boson vertices. The form factors are the similarity functions fx- The 
resulting amplitude is given by the antisymmetrization of Eq. (3.72) with dx replaced by the 
product of the form factors. In the ^-function limit, the vertex form factors equal 1 — O 21 -X 
and 1 — 6 * 12 -A- Their product equals 1 — 6 'a. Thus, the second term provides the same 
contribution as the hrst term but the factor 9x is replaced by 1 — 6 *a. 

The sum of both terms in Eq. (3.71) produces the matrix element of the scattering 
matrix on-energy-shell which is independent of A. Our complete on-shell result is equal to 
the well known one-boson-exchange scattering amplitude. 

Higher order calculations require dehnitions of the incoming and outgoing physical fermion 
states which differ from the single effective fermions due to the interactions which dress effec¬ 
tive fermions with effective bosons. Discussion of the corresponding scattering theory goes 
beyond the scope of this paper. 

There is an important property of the second order calculation above which is worth 
a separate note. When the hamiltonian width in the mass difference, as determined by 
A^ -f 2m^ or A^, becomes small then, the effective meson emission can no longer occur. Thus, 
the effective theory describes fermions interacting by potential forces. The potentials are 
given by factors f{z 22 x) 922 X in H 22 x- Note the necessary presence of the similarity form factor 
off-energy shell. g 22 x contains also the inner similarity factors which force the intermediate 
boson to form a sufficiently high invariant mass state but if the width is small enough they 
become identically equal 1 . 

In particular, f{z 22 x) 922 X in the Yukawa theory becomes equal to the relativistic potential 
which approaches the Yukawa potential in the nonrelativistic limit. The nonrelativistic 
theory is defined by the limiting case where the width A is such that the allowed energy 
transfers are much smaller than the effective fermion masses. This condition limits only 
the relative motion of the effective constituents. It does not limit their total momentum. 
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Thus, we have accomplished a derivation of the potential order in the relativistic quantum 
mechanics of fermions in the Yukawa theory. The relativistic one-boson exchange potential is 
a term in the effective hamiltonian, including its off-shell behavior. Note that in the energy 
diagonal part of the effective potential as well as in the on-energy shell scattering amplitude 
the outer similarity factor equals 1 independently of the size of the momentum transfer. In 
other words, one cannot see the onter similarity factor in the physical scattering amplitude 
order and the only trace of the effective nature of the potential are the form factors in 
the interaction vertices. 

If we use the interaction g%[)i''^^f%!)(() instead of g'ljjtljct) in writing the initial hamiltonian 
of Eq. (3.3), the resulting effective potential corresponds to the one-pion exchange between 
nucleons. Since the formalism is not limited to the nonrelativistic domain of the fermion 
momenta or to the lowest order perturbation theory, one can investigate this type of poten¬ 
tials in a wide range of applications in meson-barion and quark-pion physics. However, the 
second order calculation cannot tell us how the coupling constant depends on the width A. 
To correlate the conpling strength with the hamiltonian width one needs to calculate the 
effective hamiltonians including terms order g^. 

The remaining question is how can one explain the size of the conpling constant if it is of 
order 10. We cannot offer an explanation in the second order calculations. However, we can 
suggest that the similarity flow may lead to such large couplings on the basis of the following 
reasoning. 

The outer similarity factor and the vertex similarity factors continne to reduce the 
strength of the effective interactions when A continues to decrease. The only way the ef¬ 
fective hamiltonian of the very small width may have the same eigenvalues as the unitarily 
equivalent large width hamiltonians, is by the rise of the conpling constant. Therefore, we 
forsee that a complete calculation including the conpling constant dependence on A will lead 
to the nucleon-nncleon one-pion-exchange potentials with narrow form factors and large cou¬ 
pling constants. One may hope to derive similar structures for other meson exchanges. The 
general trend may be that the smaller is the meson mass the larger is the effective conpling 
rise because heavier effective meson emissions and absorptions are eliminated earlier. 


Fermion-antifermion interaction 

The fermion-antifermion interaction order g'^ satishes a differential equation which is 
analogous to Eq. (3.40) but more terms need to be considered explicitly. The operator 
subscripts must distinguish fermions and antifermions and one has to include terms which 
result from the annihilation channels. 

The fermion-antifermion effective interaction term is written as 

= /[-Plvyj: ^ j [^^\{yp\9iili\^\p+Kaid\i-x)P-Kai^{^-y)P-p(^A^yP+p(T2- (3.73) 

f7icr2f73cr4 

Note the change of order of the spin numbering and momentnm assignments in comparison 
to Eq. (3.42). The new order results from the operator ordering dehned in Section 2.a. 
Momenta ki and k 2 are used for fermion and k^ and k^ for antifermion operators with even 
snbscripts for annihilation operators and odd snbscripts for creation operators. 
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There are three terms contributing to the derivative of ginw with respect to A: one due 
to the annihilation channel and two due to the exchange of the boson. One of the latter 
two contributions results from the emission of the boson by the fermion and absorption by 
the antifermion and the other one is for the emission by the antifermion and absorption 
by the fermion. In each of the terms, there are two similarity functions with different 
arguments. Signs in front of each of the three terms are determined by ordering of creation 
and annihilation operators. The result is 


dgiiiix 

d\ 

- I S2r2ir25i2r24r25u-i^^^-^ 
[ [y - x)P 

The inner similarity factors are 


SiguiVs gv4U2 rnrisruru^ 


+ S^r^2'r 3512 ^" 33^3534 


6{.x - y) 
{x - y)P+ 


guiU2 gv4Vs. 


(3.74) 


S^ 


fi4 


i2) 


^Ei2 


Ai7ji 



(3.75) 


Equation (3.75) is similar to Eq. (3.44) (the subscript A is skipped for clarity). The second 
subscript of the arguments of the similarity function / denotes the vertex, i.e. 1 stands for 
the vertex where the boson was annihilated and 2 stands for the vertex where the boson 
was created. In Eq. (3.74), the fermion regularization factors hrst subscript is the same as 
the corresponding inner similarity factor subscript (i.e. the subscript of S) and the second 
subscript is the same as the corresponding fermion momentum subscript. The boson reg¬ 
ularization factors are distinguished by the subscript 5 following the convention from Eqs. 
(2.50). Their first subscript is also the same as the corresponding inner similarity factor sub¬ 
script. Last two subscripts of the boson regularization factors equal subscripts of the fermion 
momenta from the vertex where the boson regularization factor originated. Arguments of 
the regularization factors have the same subscripts as the regularization factors themselves, 
i.e. Tj = r(eej/A^). The daughter energies in the arguments are calculated according to the 
rules given in Eqs. (2.47) to (2.50). We give the results below for completeness. The same 
arguments will appear in all theories of physical interest. 


Cll — 


Pm? 


X 


ei3 = 


Pm? 

1 — X 


ei4 — 


Pm? 

y 


ei2 = 


Pm? 

i-y 


K 


621 = 


±2 

212 


-|- m 


X212 


^4 + 

1 - a;2i2 


(3.76.a) 

(3.76.6) 
(3.76.C) 
(3.76.d) 
(3.77.a) 

(3.77.6) 
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62512 — 



X 

X -212 = — • 

y 

(3.77.C) 

^^^12 = - X 212 P ^ ■ 

(3.77.d) 

624 — 

2^234 

(3.77.e) 

K-^2 I ,,2 

^234 ^ P 

62534 — T 

1 - a ;234 

(3.77./) 

1 -y 

2^234 — 

1 — X 

(3.77.(7) 

^^34 ~ ~ P ^ + 2:234 • 

( 3 . 77 ./ 1 ) 

«^312 + 

632 — 

2:312 

(3.78.a) 

« 3 ll + 

63512 — T 

1 - 2:312 

(3.78.6) 

y 

2:312 = - • 

X 

(3.78.C) 

^312 ~ P 2:312ft: . 

(3.78.d) 

ft ;41 + rri ^ 

633 — 

2:334 

(3.78.e) 

ft;^! + P ^ 

63534 — T 

1 - 2:334 

(3.78./) 

1 — X 

2:334 — 

1-1/ 

(3.78.(7) 

^334 = + 2:334^^ • 

(3.78.6.) 


Arguments of the similarity functions and energy denominators which appear in Eq. (3.75) 
are calculated according to the rules given in Eqs. (2.12) to (2.19) and (2.24) to (2.26). The 
results are universal for one-particle-exchange two-particle interactions and are given below 
for completeness. 


AA4li — {ki + /ca)^ — (/cim + ^sm)^ 

= — Cii — ei3 . 


(3.79.a) 


EM 


2 _ 
11 ~ 


-AMl^ + 2 / 1 ^. 


(3.79.6) 



AEu = AMlJP^. 


(3.79.C) 


AA 4 I 2 — {k2m + kimY — (^2 + ^ 4 )^ 

= 614 H~ 612 — /i • (3.79.(7) 

T.MI 2 = AM \2 + 2 / 4 ^ . (3.79.e) 

AE 12 = AMl^/P^ . (3.79./) 

AA^21 = (^2534/i + /i^4m)^ — k^^ 

— C 2534 + C 24 ~ • (3.80.a) 

= AMl^ + 2m2. (3.80.6) 

AE 2 i = AMlJ{l-x)P+. (3.80.C) 

AA^22 = ^2m ~ (^2512^* + /i^lm)^ 

= rri^ — 62512 ~ 621 • (3.80.6?) 

T.MI 2 = -AMI 2 + 2/i2 . (3.80.6) 

AE 22 = AMl^/yP ^. (3.80./) 

AA^g;^ = (/i:35l2/i + k2mf — ^Im 

= 63512 + 632 — rn? . (3.81.a) 

= AM^^ + 2m^ (3.81.6) 

AE 31 = AMl^/xP^ . (3.81.e) 

^■^32 — ^4m ~ (^3534/i + ^3m)^ 

= m^ - 63534 - 633 . (3.81.C?) 

= -AMI 2 + 2/i^ . (3.81.6) 

AEs2 = AMl2/{l-y)P^- (3.81./) 
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In all cases, the arguments of the similarity functions are given by Eq. (2.18), i.e. Zi = 
AA4f/(SA4f + A^) for all subscripts appearing in Eq. (3.75). 

The same reasoning is used to integrate Eq. (3.74) as in the case of Eq. (3.65) for the 
fermion-fermion interaction. For the similarity function / approaching the ^-function with 
Mo = 1/4 in Eq. ( 2 . 20 ), we have f{zf) = 6 {\‘^ + 2 mf — \AMf\) with mf = /i^ in the hrst, 
and mf = w? in the second and third inner similarity factors in Eq. (3.74). 

Integration of Eq. (3.74) gives 


fl'illiA = guile + Cl guiv^gvAU2riiri^riiri2 

+ [C 2 r2ir25i2r24r2534 0{y -x) + Cs r 32 r 35 i 2 r 33 r 3534 6 {x - y)] guiU 2 gv^v^i , (3.82.a) 

where the coefficients are. 


6'i2-11^12-A , 6'ii_a6'ii_i2 

Cl = , , . .o I -h 


I AM 


121 


I AM 


111 


C2 = 


1/6*22-2i6*22-A , (1 — x)d 2 l-\d 2 l -22 


+ 


C3 = 


{y-x)\AMl2\ {y-x)\AMl]\ 

(1 — 1/) 6*32-316*32-A 2:6*31-^6*31-32 


+ 


(X 


(3.82.6) 


(3.82.C) 


(3.82.d) 


y)\AMl 2 \ {x-y)\AMli\' 

The symbols for ^-functions have the following meaning. 6 i_j = 6*(|AM^| — |AM||) and 
6 i_\ = 6 *(|AMi I — 2 mf — A^) with mf equal y? in ci and m? in C 2 and C 3 . 

The next step is the construction of the interaction Fx [GiniA] from ^ihia of Eq. (3.73) 
using Eqs. (2.8) and (2.9). 

Then, one has to hnd out if matrix elements of Fx [Gihia] between hnite free invariant 
mass states have a limit when e —0. Stated differently, one checks if the existence of the 
limit requires the initial value of gini^ to differ from zero to cancel potential divergences in 
the limit. Following the same steps as in the case of Eqs. (3.66) and (3.70), one can check 
that no divergences arise. Therefore, ginu = 0. 

The hnal answer for the effective fermion-antifermion interaction is 


1 

'P+ 


Hiiiix — Fx [Giiha] 

^ J [a:fi:][|/p] 5'iniA/(^iniA) ^XxP+nai ^X{l-x)P-Ka 3 ^X{l-y)P-pcj4 bxyP+pa 2 -i 

<Ji(72Cr3(J4 

(3.83.a) 


where 


fi'illiA = Cl guiV'i gv 4 U 2 + [c 2 0{y - x) + C 3 9{x - y)] guiU 2 gv^v^ . (3.83.6) 

In terms of the fermion momenta. 


with s = max(M 13 ,M 24 ), and 


Cl = 


6*(s — 3/4^ — A^) 
s- y? 


(3.83.C) 
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C29{y -x)+ c^eix -y) = ^ ^ ^ , (3.83.rf) 

a 0 

with a = /i^-gi 2 , b = 9a-b = ^-9b-a = 9{m^ya-mi_^i_yb), 9a-\ = 9 [m^ya -\x- y\{2m^ + A^)], 

9b-X = 9 [rrii-xi-yb — \x — y\ (2m^ + A^)], nixy = max{x, y) and rrii-xi-y = max{l — x, 1 — y). 

The argument of the outer similarity factor in Eq. (3.83.a), Zinix, is equal to Z 22 \ from Eq. 

(3.70.h). Note that Eqs. (3.83.a-d) provide the generalization of Eqs. (3.70.a-h) to the case 
of effective interactions of distinguishable fermions. 

When A^ is reduced below 4m^ — 3 the internal similarity factor in the annihilation 
term stays equal 1 independently of the value of A. The effective interaction term provides the 
full contribution of the annihilation channel to the fermion-antifermion scattering amplitude 
of order The fermion-antifermion-boson term in the effective hamiltonian which could 
contribute to the scattering acting twice, is zero for such low values of A^ because the mass 
gap between the boson and the fermion pair is larger than A allows. 

The internal similarity factor in the exchange term becomes equal 1 independently of A 
only when A^ becomes smaller than —2w? -|- 2m/i -|- The lower bound on A^ is —m? — 

{m -|- nY (see the discussion of Eqs. (2.18) and (2.20)). In the lower bound region, the 
effective boson emission and absorption vanish and the exchange interaction term provides 
the full scattering amplitude due to the one-boson exchange. The amplitude is equal to the 
standard result on-shell where the outer similarity factor equals 1. 

If the boson mass is much smaller than the fermion mass then, for small momenta, the 
fermion energies are quadratic functions of momentum while the boson energy is a linear func¬ 
tion of momentum. Therefore, for sufficiently small momenta, the boson energy is large in 
comparison to the fermion kinetic energies and their changes. Thus, the one-boson-exchange 
interaction is mediated by a relatively high energy intermediate state. Consequently, it is 
represented by a term in the effective hamiltonian. 

For small A, the effective hamiltonian contains potentials which are equal to standard 
scattering amplitudes in the Born approximation. The potentials differ from the Born am¬ 
plitudes off-shell in a unique way which is dictated by principles of the hamiltonian quantum 
mechanics and the similarity renormalization group: the outer similarity factor reduces the 
strength of the interaction off-energy-shell. In the light-front dynamics this is equivalent to 
the off-shellness in the free invariant mass. 

For A outside the lower bound region, the scattering amplitudes obtain also contributions 
from the effective interactions which change the number of bosons by one in the transition 
through the intermediate states. Analysis of Eq. (3.71) in application to the fermion- 
antifermion scattering follows the same steps as for the fermion-fermion scattering in the 
previous Section. The resulting on-shell scattering amplitude is independent of A. The 
amplitude is equal to the well known perturbative result in Yukawa theory to order g^. 


3.C QED 

This Section describes calculations of the effective mass squared term for photons, the 
effective mass squared term for electrons and the effective interaction between electrons and 
positrons in QED. The calculated terms are order e^. 

The initial expression which we use to calculate the renormalized hamiltonian of QED is 
obtained from the lagrangian C = -|- p — m)ilj by the procedure of evaluating 
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the energy-momentum tensor and integrating over the light-front, 
we use reads 


20 The formula 




QED 


= I dx d?x 




-d^^ + 


m 


2 id+ 


- -A^d A 


Ou 


(3. 

x +=0 

where 'ipm is a free fermion held with mass m and A^ is a free massless photon held with 
Ao+ = 0. 

We replace helds and Aq{x) for x"*" = 0 by the Fourier superpositions of creation 

and annihilation operators, we order the operators in all terms and we drop terms containing 
divergent integrals which result from contractions. This is done in the same way as in the 
Yukawa theory but more terms need to be considered. Then, we introduce the regularization 
factors. 

The ultraviolet regularization factors are already familiar and the same as in the Yukawa 
theory. The additional regularization is required due to the infrared singularities. Photons 
have diverging polarization vectors when their -|—momentum approaches 0. The correspond¬ 
ing seagull term, i.e. the 5th term in Eq. (3.84), is diverging too. We have to introduce the 
infrared regularization factors, e.g. (1 -|- 5/x)~^ described in Section 2.b. We also introduce 
a photon mass term = fi'j by adding it to —8-^“^ in the second term in Eq. (3.84). The 
mass term requires a few comments. 

Consider which is hnite and let e ^ 0. The leading correction order is small in 
comparison to ixj. The situation changes when ^ 0 while e is kept hxed. At some point, 
the correction order becomes larger than /i^ itself and one cannot consider the correction 
to be a small perturbation on the chargeless theory. 

Unfortunately, the limit /i^ —0 is not completely understood in the present approach. 
The second order calculations described in this paper are not sufficient. Higher order cal¬ 
culations need to be done to gain more understanding. In particular, it is not known what 
conditions correspond to gauge invariance identities used in the lagrangian approach to make 
the massless theory exist. It is also important to include the running coupling effects. 

The width dependence of the running coupling constant can be expected to matter in 
the analysis of the initial value problem for the photon mass term. The perturbative mass 
correction for photons diverges when the ultraviolet cutoff is being removed. The hnite part 
of the photon mass counterterm is not known. Therefore, we have to ht it to data. The 
divergent term is order e^. The hnite part must be of the same order. In QED, the running 
coupling is expected to grow with the hamiltonian width. Thus, one may expect that a 
very small mass term order in low energy ehective QED may correspond to a much larger 
photon mass term in the initial hamiltonian with A = cx) due to the much larger coupling 
constant at the scale 1/e. However, the small size of the ehective hne structure constant 
makes the issue rather academic. 

Introduction of the photon mass is relevant to the ehective infrared dynamics. The cutoh 
6 allows photons with x ~ 5 to appear in the ehective interactions of a hnite width A. Such 
a photon may contribute only a little to an intermediate state energy if it has sufficiently 
small transverse momentum. Therefore, emissions and absorptions of such photons are 
allowed even by the narrow similarity factors. The massless photons with small x introduce 
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dependence on the cutoff S into the effective hamiltonian. Moreover, large numbers of such 
photons can appear in the eigenstates. But large occupation numbers for bosons enhance 
the strength of the emission and absorption processes. This enhancement could, in principle, 
compete with the smallness of the coupling constant. A hnite photon mass suppresses these 
effects and, as a consequence, may become a major alteration of the theory. On the other 
hand, the mass may remain to be merely a stabilizing factor in the calculations if the effective 
dynamics turns out to be dominated by a limited range of photon wavelengths so that the 
momenta order S or smaller do not matter. This is certainly expected to happen in the hnite 
size neutral bound states. 

A hxed value of fis leads to the conclusion in perturbation theory that the photon eigen¬ 
states have masses equal to fis when the charge e approaches 0. Therefore, we will be forced 
to consider the limit —>• 0 in order to discuss physical photons to order e^. Also, the 
nonzero mass squared term for photons leads to additional divergences when <5 —0 and the 
limit of /i 5 —0 removes those. 

The infrared hniteness of QED suggests that physical results may be independent of fis 
when it is sufficiently small. But we do not provide an explanation of how this comes about 
in the light-front hamiltonian approach in general. We introduce the photon mass fis and 
investigate the limit /i^ —>• 0. The lowest order calculations in this paper lead to results 
which are independent of when it tends to zero. 


Photon mass squared 

The same procedure from Section 2 which led to Eqs. (3.4) and (3.5) in the case of 
mesons in Yukawa theory, leads in QED to 


Ql photon \ ^ ) 


+ h! t 


k+ 


^ka^kcr 


(3.85) 


No correction arises to the term because our regularization preserves the kinematical 

symmetries of light-front dynamics. 

A new feature in comparison to the Yukawa theory is the polarization of photons. With 
the kinematical symmetries explicitly preserved, only terms diagonal in the photon polar¬ 
ization emerge. For example, terms proportional to with cxi ^ <72 cannot appear 

because the regularization and similarity factors do not introduce dependence on the photon 
momentum. Note that such terms are allowed by the power counting. 

The net result of the photon self-interaction is an effective photon mass squared term 
which is independent of the photon momentum but varies with the effective hamiltonian 
width A. One obtains more complicated results for the effective photon free energy if the 
regularization or similarity factors violate kinematical symmetries of light-front dynamics. 

i 

The dependence of /i^ on A is determined to order by the equation 
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dX 


A^ai<72 ^ g2 


XK, 


dp{zl)Tr + m) ^kaAhm-m) 


dX 




Ps 


r^{x,K) 


(3.86) 


and the initial condition at A = cxd. However, it is also sufficient to know the effective mass 
squared at any single value of A to determine its value at other values of A using Eq. (3.86). 
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The initial condition at A = cx) is distinguished only because it provides connection with 
standard approaches based on the local lagrangian for electrodynamics. 

In Eq. (3.86), + —x), — and + so that 

2:a = —+ In the limit of Eq. (3.8), we have/^( 2 :f) = 6{\^ + 3n‘j — M‘^). 

In fact, Eq. (3.86) is free from infrared singularities and we could skip the introduction of 
/i| by letting it go to zero at this point. However, the systematic approach defines the 
hamiltonian of QED including the infrared regulator mass for photons and we can keep 
it here for illustration. The regularization factor r^{x,K) in Eq. (3.86) is the same as in 
Eq. (3.9) in Yukawa theory because only fermion regularization factors enter Eq. (3.86), 
according to Eqs. (2.47) to (2.49), and these factors are the same in both theories. 

Evaluation of the spin factor gives 


d\ 


= e 


dp{zl) 


r,(a;, k). 


(3.87) 


which is the QED analog of Eq. (3.6) from Yukawa theory. 

Integration of Eq. (3.87) is carried out through the same steps as in the Yukawa theory. 
We can use Eq. (3.22) to calculate the effective photon mass squared knowing its value 
/Xq at some value of A = Aq. 

The value of /Xq is found by requesting that the effective hamiltonian eigenvalues for 
photon states contain the physical photon mass /i, expected to be equal 0. However, solving 
the eigenvalue equation to second order in the coupling constant e through the same steps 
as in the case of mesons in Yukawa theory in Eqs. (3.23) to (3.25), leads to the physical 
photon mass fi = ns. fis is small but finite and it is considered, in terms of powers of e, to 
be of order e° = 1 when e —0. 

In the ^-function limit for the similarity function / with uq = 1/4, 


d'o — d'S 4- 
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Att 


dx / dn 


2M^ - 4k 


±2 


0 x(l — x)M^ — ju^ 


0(A^ + 3/i^-Af2)+o(e^), 


(3.88) 


Thus, at other values of A, the effective photon mass squared is 


+ S i' r x)m'^ - + o(e‘). (3.89) 

This result naturally depends on the infrared regularization parameter /i^ but no singularity 
appears when this parameter is set equal to zero. For A^ + 3/i| < 4m^, where 4m^ is the 
lowest possible free invariant mass squared for the two intermediate fermions, the photon 
mass is independent of the hamiltonian width A and it equals /i|. For larger values of 
A^, the effective photon mass grows with the width A so that its larger value compensates 
effects of the interactions which become active for the larger width. The net result is that 
the photon eigenstates have eigenvalues with masses squared equal /i| independently of A. 
Finally, the result favored by experimental data is obtained in the limit /i| —0 at the end 
of the calculation. 


Electron mass squared 
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Electron and positron self-interactions through emission and reabsorbtion of transverse 
photons lead to the fermion free energy terms of the form exactly the same in QED as in Eq. 
(3.27) in Yukawa theory. However, the effective mass of electrons and positrons depends on 
the width differently than in the case of Yukawa theory. Instead of Eq. (3.28), one obtains 
now 


dWj^ _ 2 

d\ ~ 

G 

where k = {k — p)o, = xk'^, p-^ = xk^ + k-*-, = {m? + k^)/x + (/i^ + — x), 

— + and z\ = AAf^/(SAf^ + A^). The regularization factor 

of Eq. (2.49) for the intermediate particles and the infrared regulator for the intermediate 
photon, as given at the end of Section 2.b, imply 
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(3.91) 


The sum over photon polarizations in Eq. (3.90) produces the well known expression 
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(3.92) 


and the spin factor in Eq. (3.90) is 
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(3.93) 


The new feature of this expression, in comparison to the Yukawa theory, is the divergence 
for a; —> 1, i.e. where the photon longitudinal momentum approaches 0. 

The rate of change of the electron mass term versus the effective hamiltonian width in 
the ^-function limit for the similarity function with mq = 1/4 is following. 
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(3.94.6) 


AI^ = . + ^ + ^ . (3.94.C) 

X 1 — x 

The divergence structure of the effective electron mass in Eq. (3.94.a) is obscured by 
the fact that the whole effective mass term is merely a number dependent on A while three 
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cutoff and regularization parameters appear in the integral: e, 6 and j^s- The only available 
condition is that the effective electron mass should have a limit when e —0. However, this 
condition has to be satisfied without generating divergences in the physical electron mass 
(i.e. in the electron eigenvalue energy) when we remove the infrared regularization. Since 
other contributions to the physical electron mass may diverge as 6 or fis tend to 0, and only 
the sum is hnite in the limit, one needs to keep track of the infrared structure in dehning 
the e-independent (i.e. ultraviolet hnite) part of the counterterm. 

The divergences due to e —0, <5 —0 and ^ 0, are not resolved in the single mass 
constant. Many elements of a complete analysis overlap in producing the hnal answer and 
many simplihcations are possible. We will proceed in this Section with a simplihed analysis. 
A more extended analysis will be required for other hamiltonian terms where the outcome 
of the procedure is not reduced to hnding only one number in the effective interaction. 
For example, in the electron-positron interaction term, the external momenta of fermions 
introduce a whole range of additional parameters. That case will be illustrated in the next 
Section. 

The (5-function under the integral on the right-hand side of Eq. (3.94.a) forces to be 
equal -|- 3m^. The smallest possible value of At is m + fis and the negative value of A^, 
namely, (m -|- fisY — 3w? = —rn? — {m + fisY + 2[(m -|- fisY — is required to reach this 
lower bound (see the comments about Eq. (2.20)). Below this bound the right-hand side 
of Eq. (3.94.a) vanishes, no emission or absorption of photons by electrons is possible and 
the effective electron mass stays constant. The smallest possible value of A^ allowed by Eq. 
(2.18) is —m? — {m + fisY- The difference between these bounds vanishes when the photon 
mass goes to zero. 

In the next Section, we will also consider the case of the inhnitesimal Uq (see the discussion 
below Eq. (2.20)), which leads to (5(A^ -|- — Ai"^) in Eq. (3.94.a), instead of 5{3m^ -|- A^ — 

M^). 

For A^ > (m -|- fisY ~ the right-hand side of Eq. (3.94.a) is positive. Therefore, 
the effective electron mass squared term grows together with the width of the hamiltonian. 
We will see below that this growth combines with the growing negative contributions of the 
corresponding effective transverse photon emission and reabsorption, so that the physical 
electron mass is independent of the hamiltonian width. 

For A^ close to the lower bound, u must be close to 0 and a; ~ 1 — fis/ijn + fis)- But 
Ai^ is limited and determined by the value of A. Quite generally, as long as Ai remains 
limited, u and /^^/(l — x) are limited. The invariant mass denominator in Eq. (3.94.a) 
equals A^ -|- 2im?. The denominator is small only when A^ is close to the lower bound. Then, 
the integration range is small too. The denominator, when expressed in terms of A, can 
be pulled out and put in front of the integral. The integration over u sets u = u{x) = 
A^ -|- 3m^ — m? jx — fi^/{I — x) and forces the condition u{x) > 0. This condition implies the 
following limits on the integration over x, provided A^ -|- 3im? > {m + fis)"^ since otherwise 
the integral is 0. 


xq — Ax < X < xq + Ax , 


(3.95.a) 



A^ -|- 3m‘^ j ’ 


(3.95.6) 
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(3.95.C) 


Within these limits, u{x) varies from the minimal value of 0 at the lower bound Xq — Ax 
through a maximum of + 3m^ — (m + fisY at a; = m/{m + ns) to the minimal value of 0 
again at the upper bound xq + Ax. In the case of infinitesimal mq, one replaces A^ + 2w? in 
these formulae by A^. 

Eq. (3.94.a) reads 
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The upper limit of integration over x for jx^ much smaller than and 2m^ + A^ , equals 
1 —/i5/(2m^ + A^) and approaches 1 when 0. For x close to 1, the factor (1 —a;)“^ in the 
square bracket of the integrand is large and leads to a logarithmic dependence of the integral 
on the upper integration limit. The logarithm would become inhnite for 0 if 5 were 

equal 0. Therefore, the limit /i^ —0 is sensitive to the presence of the regularization factor 
with 6^0. For the hnite fj,s, the region of a; —1 is regulated by e and can be considered an 
ultraviolet limit. A counterterm to the diverging e dependence could remove the divergence 
due to X ^ 1. Then, a separate cutoff parameter S would not be needed but the resulting 
terms would diverge for /i| —0. For finite /i|, the derivative of the electron mass with 
respect to A is hnite. For = 0, the integrand in the region a; ~ 1 is regulated solely by 
the infrared regularization factor [1 + 5/(1 — x)]~‘^ since the upper limit of integration over 
X is equal 1. 

For hnite A^, all three terms in i.e. u{x), w? jx and /^^/(l—a;), are limited. Therefore, 
for hnite one can take the limit e —0 in the integrand. The factor [1—5/(1—a;)]“^ remains 
and additionally cuts oh the integration region at a; ~ 1 — 5. The ratio of 5 to ixl/{2rn? + A^) 
determines the size of contributions obtained from the upper range of integration over x. 
For hnite hxed values of one can take the limit 5 —0 and log [(2m^ + A^)//^^] appears 
in the answer. 

The right-hand side of Eq. (3.96.a) contains terms which behave for large A as a constant, 
as A“^ with factors of logarithms of A and as functions vanishing faster than A“^. One 
integrates Eq. (3.96.a) over A from Aq to inhnity in order to express the ehective electron 
mass squared term at Aq, denoted ttIq, in terms of the initial m^. Clearly, the integration 
over A would diverge without the regularization factor which depends on e. The integration 
produces terms behaving as e~^, loge and terms convergent in the limit e —> 0. in 


67 



the initial hamiltonian must be supplied with a counterterm to subtract the diverging e- 
dependent terms in the effective hamiltonians. 

In summary, the infrared divergence due to 0 and 5 0 appears in the derivative 

of m\ with respect to A. Therefore, even if one requests that the electron mass term is hnite 
at some value of A, the effective masses of electrons in the neighbouring hamiltonians with 
even slightly different widths will diverge when fxs and S approach 0. We have to abandon 
the requirement that the effective electron mass term at any value of A remains hnite when 
the infrared regularization is removed. The effective masses diverge in the limit fis ^ 0 and 
5 —0. The only condition we can fulhll through the ultraviolet renormalization is that the 
effective electron masses for hnite A are independent of e. 

Mathematical details of the ehective electron mass term calculation are more compli¬ 
cated than in the Yukawa theory because the infrared regularization parameters are present. 
Otherwise, the calculation is essentially the same and we skip the description here. We only 
stress that the counterterm and the ehective masses of electrons and positrons depend on 
the infrared cutohs and they diverge when the cutohs are being removed. 

Thus, the ehective electron mass squared term in the limit e —0 is 
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(3.97) 

The hnite term m^ has a limit when e —0. Its dependence on the infrared regularization is 
not displayed, m^ is found from a suitable renormalization condition. 

The natural condition to be satished by m^ is that the ehective hamiltonian of some 
width A has the electron eigenstates with eigenvalues equal + m‘^)/p'^, where p denotes 
the electron momentum and m is the physical electron mass. The problems of measuring 
the physical electron mass and comparing the result with the eigenvalue of the hamiltonian, 
are not considered here. 

The eigenvalue equation for electrons can be solved in perturbation theory in the same 
way as for bosons in Yukawa theory in Eqs. (3.23) - (3.26) and fermions in Eqs. (3.35) - 
(3.39), or in QED for photons in Eqs. (3.88) - (3.89). 

One obtains the condition 
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and one can calculate mo from this condition. Consequently, 
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and m^ = m^ -|- o(e^). 

The physical electron mass is independent of the infrared regularization because the 
regularization dependent ml and the effective emission and absorption of photons combine 
to the regularization independent result. Note that the effective mass squared term for 
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electrons depends on the infrared regularization parameters and it diverges when the infrared 
regularization is removed. Higher order calculations are required to discuss contributions of 
the infrared photons in the experimental determination of the electron mass. 


Electron-positron interaction 

Calculation of the effective electron-positron interaction to order is of general interest 
as a fundamental force in electrodynamics - this interaction is responsible for the formation 
of positronium. Also, effective interactions between quarks and antiquarks in lowest order 
QCD effective hamiltonian will have similar features and the present calculation of the 
electron-positron interaction in QED provides an introduction to the QCD case. 

Generally speaking, the QED calculation of the effective electron-positron interaction 
proceeds in the same way as in the case of fermion-antifermion interaction in Yukawa theory 
except for three new elements. 

The first is that photons have polarization vectors which enter in the vertex factors and 
introduce additional dependence on the exchanged photon momentum. This dependence 
leads to infrared divergences for small longitudinal momenta of exchanged photons. 

The second feature is that the infrared divergences require additional regularization fac¬ 
tors. We use the parameter 6 and the photon mass squared jus ^ 0. The limits <5 —0 and 
/i 5 —0 are generally understood as to be taken at the end of a calculation of observables and 
not in the effective hamiltonian itself. However, it may also be possible to take the limits 
in the matrix elements of the hamiltonian between the states which do not induce infrared 
divergences, i.e. do not involve small x photons. 

The third feature is that the one-photon exchange interaction needs to be combined 
with the 5th term from Eq. (3.84) to obtain the standard results for the electron-positron 
scattering in the Born approximation. The 5th term from Eq. (3.84) provides the initial 
condition for the renormalization group flow of the effective hamiltonians. In order e^, this 
term is only supplied with the outer similarity factor by the operation Fx. It does not change 
in the flow otherwise because it is order itself. The initial condition provides an important 
contribution which leads to the Coulomb interaction. The one-photon exchange by itself is 
not sufficient to produce the effective Coulomb potential. 

This is a different situation than in Yukawa theory where no four-fermion seagull in¬ 
teractions appeared and the one-meson exchange interaction was sufficient to produce the 
Yukawa potential in the effective hamiltonians of small widths. 

The effective electron-positron interaction term has the analogous structure as the fermion- 
antifermion interaction in Eq. (3.73), i.e. 
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which is the QED analog of Eq. (3.74) from Yukawa theory. Notation is the same as in 
Eqs. (3.74) to (3.81) with the exception that is replaced by /i|. The new elements are 
the infrared regularization factors of Section 2.b, i.e. r^,/^ = r{k'^5/k'^), = r(A;^5/A;^), 

^ 1/5 = 5/k'^), r 4/5 = r{k'l^5/k'^), and the photon polarization vectors. The sum over the 

photon polarizations gives 
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The terms proportional to the four-vector k^Q can be rewritten using the Dirac equation 
for free fermions of mass m into a simpler form. For example, in the second term on the 
right-hand side of Eq. (3.101), we have 
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Using similar relations for all vertex factors involved one obtains 
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We use the notation Sij = {ki + kjY and qX = {ki — k^. 

Integration of Eq. (3.104) proceeds in the same way as in the case of Eq. (3.74). The 
initial condition at A = cx) includes the seagull term. 
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The result of the integration of Eq. (104) is following. 


(3.105.6) 
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(3.106) 


The coefficients Ci, C 2 and C 3 are given by universal Eqs. (3.82.b) to (3.82.d) which were 
derived already in Yukawa theory, with the replacement of the meson mass by the photon 
mass, = yj. 

This result will be now analysed term by term for illustration in the matrix elements 
between states of equal free energy. The matrix elements contribute to the electron-positron 
scattering in second order perturbation theory. 

The first term in Eq. (3.71) for the T-matrix calculated to order e^, has matrix elements 
equal to the matrix elements of the effective interaction from Eq. (3.106). When evaluat¬ 
ing the S'-matrix elements, one considers conhgurations where the free energy of incoming 
fermions equals the free energy of the outgoing fermions. This conhguration selects the 
energy-diagonal part of the effective interaction: S 13 = S 24 = s and ^43 = Q 34 = 5 ’^- Also, the 
external similarity factor which appears in the effective haniltonian as an additional factor 
to equals 1. Thus, fi'iniA in fhe energy diagonal part can be viewed as the scattering 

amplitude. It simplihes in the energy diagonal part to the following form. 


^ ^counterterm ^ 

fi'nllA = guile - CieUiYV3 eV4'y^U2 


+ —CiS -|- 1 ] 


eMi7’^n3 en47’^M2 

p+2 


C2r2512r2534r3/br2/b9{y -x) + C3r35i2r3334ri/3r4/39{x - y) 


eui'y^U 2 ev4'yf,V3 


{ C 2 q^ -h 1 r25i2r2534»^3/5^^2/56'(|/ “ x) 


+ 


C 3 q^ + 1 ] r 35 i 2 r 3534 ri/ 5 r 4 / 56 '(a; - ?/)} 


eMi7'''M2 en47’''n3 

{x — |/)2P+2 


(3.107) 


We have removed the regularization factors which equal 1 for fermions with hnite free energy. 

In the limit /i^ —0 and for the cutoff A close to —2m^, a number of simplihcations 
occur. We display the result for the case where the intermediate photon momentum fraction 
\x—y\ » {x, y, 1—x, l—y)5, i.e. when the photon momentum is not negligible in comparison 
to the fermion momenta. In this case, the infrared regularization factors for the photon which 
are still kept in Eq. (3.107), can be set equal 1. The result is following. 


^_ _ ^counterterm ^3 ^^^2 

yiiiix — yiiiie - 


+9 \q‘^\ — \x — y\{2m‘^ -|- \‘^)/max{x, y,l — x,l — y) 


eUiYU 2 eV4'^^eV3 
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(3.108) 


-e 


X — y\{2m^ + \^)/max{x, y, 1 


x,i-y)- 


eui'y^U2 ev4'y^V3 
{x — |/)2p+2 


The first term is the potentially necessary counterterm which we have not yet determined. 
Since the remaining part of the matrix element is not sensitive to e, the counterterm matrix 
element is equal zero. 

The second term is the well known expression for the electron-positron annihilation chan¬ 
nel scattering amplitude in the Born approximation. No limits on the fermion momenta 
appear because the invariant mass squared of two fermions is larger than 4m^ which is the 
minimal invariant mass difference in the transition between the two fermions and a one 
massless photon state. Since we assume the width to be small, the effective hamiltonian 
contains the full amplitude for transition through the intermediate photon state. 

The third term is equal to the standard second order expression for the electron-positron 
scattering amplitude via one-photon exchange, except for the 6'-function factor which forces 
the momentum transfer to be sufficiently large. 

The meaning of this restriction is visible in the T-matrix. The third term contributes to 
the e+e“-scattering amplitude through the hrst term in Eq. (3.71). The same contribution 
would originate from the second term in Eq. (3.71) if we were using the initial hamiltonian to 
calculate the T-matrix. In contrast, the effective hamiltonian with the small width A limits 
the effective photon emissions and absorptions to small momentum transfers and, therefore, 
it is not able to provide this contribution through the second term in Eq. (3.71). This 
contribution is then contained in the effective hamiltonian and comes through the first term 
in Eq. (3.71). 

The fourth term is unusual in the sense that it should not appear in the electron-positron 
scattering at all. The fourth term distinguishes the 2 :-axis in its structure and diverges when 
X ^ y. The ^-function factor in the fourth term is equal 1 where the 6'-function factor of the 
third term is equal 0. And vice versa, the fourth term ^-function equals 0 where the third 
term ^-function equals 1. 

The need for the fourth term becomes clear when one recalls that the second term in 
Eq. (3.71) also contributes to the electron-positron scattering amplitude. The relevant 
contribution comes through the effective one-photon exchange which results from the double 
action of Hjx- Hi\ is given by the operation Fx applied to the third term of the QED 
hamiltonian from Eq. (3.84), (see Eq. (2.11)). The operation Fx multiplies the photon 
emission and absorption vertices by the factor fx- This factor was set equal to a ^-function 
in the current example. Each interaction provides one factor of the 6'-function. The resulting 
factor in the second term of Eq. (3.71) is the same in QED as in Yukawa theory in Eq. (3.72), 
except for the antisymmetrization effect which leads to the 6'-function which stands in front 
of the fourth term in Eq. (3.108). 

Now, the second term in Eq. (3.71), contains the spin factor which is the same as in 
the last term in Eq. (3.104). The part complements the third term in Eq. (3.108) and 
produces the full well known one-photon exchange scattering amplitude which is free from 
the 6*-function factor. The remaining part provides the term which cancels the odd fourth 
term in Eq. (3.108). Thus, the effective hamiltonian calculated to second order in powers 
of the charge e contains an odd term and 6'-functions which are required to compensate for 
the odd contributions to the scattering amplitude from the effective hamiltonian order e 
acting twice. The above analysis of the energy diagonal part of the second order effective 
hamiltonian explains the role of different terms in Eqs. (3.106) to (3.108). 
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The analysis also suggests that apparently infrared diverging terms in the effective hamil- 
tonian may mutually compensate their diverging contributions in the scattering amplitude 
on energy-shell. In the current example, we see the interplay between the second-order 
seagull term and the double action of the hrst order emission and absorption of photons. 
The hrst order hamiltonian matrix elements diverge when the photon longitudinal momen¬ 
tum approaches zero. The second order seagull term compensates this divergence in the 
on-energy-shell T-matrix elements. 

The remaining point to make here is that the result of Eq. (3.106) with the counterterm 
equal 0 leads to the effective light-front hamiltonian version of the Coulomb force in the 
limit of small -|- 2m^ -C am?. The key elements in deriving this conclusion are the outer 
similarity factors and the smallness of a. The outline of the derivation is following (c.f. Ref. 

i)- 

If only the small energy transfers are allowed by the outer similarity factor, i.e. transfers 
much smaller than the electron mass, then, the wave functions of the lowest mass eigenstates 
of the effective hamiltonian are strongly peaked at small relative electron momenta and they 
fall off very rapidly as functions of the relative momentum. This is not true without the outer 
similarity factor because the function ^finiA alone is too large at the large energy transfers 
and it would produce singular contributions in the large relative momentum region, sensitive 
to the regularization cutoffs. 

Below the width scale, the wave functions fall off as dictated by the eigenvalue equation 
with small a. Above the width scale, the fall off is very fast due to the similarity fac¬ 
tor which justihes restriction to momenta much smaller than m. Then, the nonrelativistic 
approximation for all factors in Eq. (3.106) becomes accurate. 

In the nonrelativistic approximation, = Qu ~ ^0- (3.107) applies. Further, in 

Eq. (3.108), the ^-functions become effectively equal 1 and 0, respectively. The last term is 
not leading to important contribution despite its divergent longitudinal structure because it 
is canceled by the effective massless photon exchange as described earlier in this Section. 

The dominant contributions are provided by the second and the third terms from Eq. 
(3.108) which are well known to have the right nonrelativistic structure for predicting positro- 
nium properties in the Schrodinger equation. Now, the outer similarity factor becomes ir¬ 
relevant to the spectrum in the leading approximation because the coupling constant is very 
small. In the dominant region, the electron velocity is order a, the nonrelativistic approx¬ 
imation to the full dynamics produces wave functions with relative momenta order am (in 
the positronium ground state the wave function is (fc^ -|-a^m^/4)“^) and the outer similarity 
factor in the effective interaction can be replaced by 1. 

We can use the inhnitesimal mq in Eq- (2.20) and replace -|- 2m^ by A^ (see the 
discussion below Eq. (2.20)). When A is order am and x — y is order a then, the momentum 
transfer q‘^ is typically order which is much larger than 2{x — y)\^ in Eq. (3.108). 

Thus, the ^-function is equal 1 and the second term in Eq. (3.108) becomes equal to the 
standard Coulomb interaction with the well known Breit-Fermi structure of the spin factors. 
This step completes the derivation of the Coulomb potential. The derivation explains the 
effective nature of the Schrodinger equation with the Coulomb potential in the light-front 
hamiltonian formulation of QED. 

3.d QCD 
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The main conceptual complication in hamiltonian calculations in QCD is confinement, 
which is not understood. To order the manifestation of the confinement problem is the 
lack of a well defined initial condition for the renormalization group flow of the effective 
hamiltonians. It will require an extended research effort to find the class of acceptable initial 
conditions. 

For example, the renormalization conditions for the quark and gluon mass terms are 
missing in the second order calculations in this paper. We do not know the right values of 
the masses in the initial hamiltonian, and we may currently only speculate about their values 
in the small width hamiltonians (by comparison with the constituent quark model). So, we 
limit our presentation here to the differential equations themselves, without solutions. 

We stress the urgent need for the higher order calculations by describing some details 
of the second order calculation of the qq effective interaction. The calculation is similar to 
the one in QED above, except for the option for a differnt treatment of the last term in Eq. 
(3.108). Perry suggested that the long distance part of this term may remain uncanceled 
because of the gluon non-abelian gauge interactions. If the last term in Eq. (3.108) would 
not be canceled, it could be claimed to generate confinement in the light-front hamiltonian 
approach to QCD. We describe the structure of this term in the present approach since 
it is different than in Ref. P]. The differences result from the diferent definitions of the 
similarity transformation. 


Quark and gluon mass terms 

Results one obtains in QCD from Eq. (2.38), in second order in Q 2 \, can be illustrated 
by two equations for the hamiltonian width dependence of the effective masses of quarks 
and gluons. It is arbitrarily assumed that the equations are relevant to physics although 
their structure follows from the second order terms only. The range of widths and coupling 
constants for which these equations can pertain to the physics of quarks and gluons, are not 
known yet. 

The regularization factors are not written in detail. One can write them easily using 
results of Section 2 and previous examples in Section 3. In the case of effective masses a 
number of simplifications occur, as explained in the previous Sections. Let us consider an 
infinitesimal uo in Eq. (2.20) and simplify notation by replacing A by A itself. Then, we can 
write 


dX 


Gi\ — 


dp{zx/X^)/dX 

Pl2A- t;- z:, -P 2 IA 


Qix - E 


lA 


11 


Eix is the eigenvalue of Gw corresponding to the subscript 11. 
needed. Namely, 


(3.109) 

A set of arguments is 
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(3.110) 

(3.111) 

(3.112) 
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m\ and /xa are the effective quark and gluon masses, respectively. Then, 


dm\ 

~d\ 
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IxH 




dX 


k'^\2/x +A/{1 —xY]+2m\{l —xY/x TqgYx^n) (3.113) 
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^-i dp{z2/X^ 

dX 


+ m\ 

x{l — x) 


2k^ 


VqqYx, k) 


(3.114) 


The gluon couples to the quark-antiquark pairs and pairs of gluons while the quark couples 
only to the quark-gluon pairs. The running of the quark and gluon masses is more involved 
than in the previous examples. In Eq. (3.114), the number of colors is equal to 3 and the 
number of flavors to 1. 

These equations are not further studied here for two major reasons. The hrst one is 
that we do not know the initial conditions to use for such study. The second is that the 
equations involve two running couplings which are not known yet. The third and fourth 
order calculations are required to hnd them, at least approximately, before one will be able 
to make contact with data for hadrons using the effective hamiltonians of small widths. 

The importance of the effective mass issue for quarks and gluons is illustrated below by 
the calculation of the small energy transfer effective forces between quarks and antiquarks. 


Quark-antiquark effective interaction 

The second order results are similar to QED. For heavy quarkonia, one can directly look 
at Eq. (3.108). The only change required is the color SU{3) matrices sandwiched between 
color vectors of quarks and summed over colors of the exchanged gluons. The counterterm 
is set to 0. The second term gives the annihilation channel potential but the color matrix is 
traceless and this excludes the single gluon annihilation channel from the dynamics of color 
singlet QQ states. The third term leads to the color Coulomb potential with the well known 
Breit-Fermi spin factors. 

About the last term in Eq. (3.108), it was suggested P that a term of this kind may 
become a seed for conhning interactions if it is not canceled. The cancelation occurs when 
one evaluates the model interaction in the QQ sector in the explicit expansion in powers of g 
to second order using the operation R and Eq. (1.3) in the nonrelativistic approximation for 
the quark relative momentum, exactly the same way as for electrons in QED with massless 
photons. 

However, there are reasons for gluons to effectively acquire a mass, such as the unknown 
hnite part of the mass squared counterterm or the interaction energy due to the non-abelian 
terms which are absent in QED. Then, the cancelation of the seagull term can be questioned. 

One can assume that the gluon energies in the QQg sector may be lifted up so that the 
gluons cannot contribute to the model QQ interaction in the way the abelian massless pho¬ 
tons can in the model electron-positron interaction in QED. As a result of this assumption, 
one obtains the last term in Eq. (3.108) acting in the effective QQ sector. 
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In the nonrelativistic approximation, the Coulomb term and the term in question are 
(see Eq. (3.108)) 


X 


1/1 2A' 


g'^Arn? 


\x — y\ 2A^ — g" 


X 


yf 


(3.115) 


The regularization factors are the same in both terms and they are not displayed. The only 
effect of their presence which matters here is that \x — y\ is limited from below by about 
5/2. The initial infrared regularization parameter 5 comes in through the initial condition in 
the renormalization group flow of the seagull term. The flow is limited in the second order 
calculation to the dependence of the outer similarity factor on A. The factor 1/2 results from 
\x — y\/xquark being limited from below by 5 and the quarks having Xquark ~ 1/2. In fact, the 
lower bound on |x — y\ is given by 25max{x, 1 — x,y,l — y). This is different from Ref. [§ 
where instead of the ratios of the +-momentum fractions a separate scale for +-momentum 
is introduced. 

The third component of the exchanged gluon momentum is = {x — y)2m. Thus, we see 
that the uncanceled singular term is represented by the potential which is analogous to the 
Coulomb potential except for that the factor — 1/g^ is replaced by — l/g| and both terms 
have mutually excluding and complementary supports in the momentum transfer space. 

The seagull term ^-function can be rewritten as 6'[a;^ — (|g3|—a;)^—g-*“^], where a; = /2m. 

The support of this function is two spheres of radius to centered at g-*- = 0 and gs = iix. 
The spheres touch each other at the point g-*- = gs = 0. In this point, g| in the denominator 
produces a singularity. 

Let us initially consider both terms in Eq. (3.115) as the actual interaction in the model 
QQ sector, i.e. as if they were not affected by the operation R in Eq. (1.3). The Coulomb 
term works outside the two spheres in the ^-space and the singular seagull term works inside. 

In the region of the singularity, both g-*“ and are small in comparison to u. In the 
rough analysis, one can neglect the outer similarity factor 6{X‘^ — — fc'^|) since it is equal 

1 when q = k — k' approaches 0. k is the relative momentum of the created QQ pair and 
k' is the relative momentum of the annihilated QQ pair, uj = (A/m)A/2 A/2 and the 

spheres have the radius about A/m times smaller than the outer similarity factor width in 
the length of the quark relative momenta, i.e. the relative size of the spheres in comparison 
to the outer similarity factor support approaches 0 when A/m —> 0. 

Since g-*-^ is order gs in the singular region, the divergence when gs —>■ 0 is logarithmic. 
The lower limit of integration over |g3| for a given x is given by 2m 5 max{x, 1 — x). However, 
we assume x = 1/2 + o(g^) and we neglect terms of higher order than in the model QQ 
hamiltonian. 

The potential resulting from the uncanceled seagull term is given by the following ex¬ 
pression (c.f. 1^), 


V{r) 


d?q 

j2Xf 


exp {i qr) 


6'(2cc;|g3| - g^) 6'(|g3| - 2m5) 




(3.116) 


The sign ~ means that the diverging dependence on 5 is subtracted and the same coefficient 
stands in front of the integral as in the Coulomb potential term. The argument for the 
infrared subtraction goes as follows. 

If the gluons cannot cancel the last term in Eq. (3.108), they presumably cannot con¬ 
tribute to the model quark self-energies either, for the same reason. Because the size of 
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the quark mass in the effective hamiltonian is unknown, one may propose that its value is 
chosen in the second order calculation in the same way as for nucleons in the Yukawa theory 
in Eq. (3.37) or electrons in QED in Eq. (3.99): so that a would-be quark eigenstate has 
a finite constituent quark mass when the gluons are allowed to contribute at all momenta 
in the eigenvalue equation. This setting is equivalent to the solution Perry proposed for his 
coupling coherence condition for the quark self-energies. The argument also illustrates 
the urgency of questions concerning the initial conditions and higher order analysis. 

There is nothing wrong with the mass adjustment despite the infrared divergence. We 
have noticed in the previous Section that the arbitrary finite parts of the ultraviolet countert¬ 
erms can be infrared divergent. This time, however, the positive and infrared logarithmically 
divergent part of the effective quark mass term in the model eigenvalue equation for heavy 
quarkonia remains uncanceled when the transverse gluons with 2a; | ga | — > 0 are declared 

to be absent from the model dynamics. The uncanceled part of the effective quark mass term 
stands in the eigenvalue problem. The point is it can now cancel the diverging 6 dependence 
in the seagull term which is not canceled due to the missing of the gluon exchange below A. 
The cancelation occurs in the colorless states. It is analysed here in the nonrelativistic limit 
only. 

We describe the cancelation mechanism in the case of equal masses of quarks. The 
mechanism is similar but not identical to that in Ref. |P]. The infrared divergent mass 
squared term comes into the quarkonium eigenvalue equation divided by a;(l — x). But in 
the second order analysis the mass divergence appears only as a logarithmically divergent 
constant. The x-dependence is of higher order. The same diverging constant with the 
opposite sign is generated by the seagull term. 

The infrared divergent terms and their cancelation are not directly related to the ultra¬ 
violet renormalization procedure. They appear in the ultraviolet-finite effective small width 
hamiltonian dynamics. Note also that the introduction of the gluon mass fis in the regular¬ 
ization could matter for the lower bound on \x — y\ and it could even eliminate the whole 
contribution when the upper bound of uj meets the lower bound of fis- We assume here 
fis = 0 . 

The divergent part in Eq. (3.116) is independent of E and it is easily removed by sub¬ 
tracting 1 from exp zgF. Evaluation of the integral leads to the answer that for large r the 
seagull term produces a logarithmic potential of the form 

R(r) = ^ (3.117) 

71 

where a is equal 1 for the radial versor along the x-axis and it equals 2 when r is purely 
transverse. 

The resulting potential is confining. It is also boost invariant. But the asymmetry of 
the potential raises doubts. It suggests that an important piece of physics is missing in the 
reasoning used to derive it. The obvious sources of questions are the role of the operation 
R, the role of the nonrelativistic approximation, the size of the quark and gluon masses 
and the strong dependence of the term on the width A. The most urgent question is if the 
higher order calculations do lead to the uncanceled seagull term. The width dependence 
must disappear in the eigenvalues. The coupling constant dependence on A is necessary. 
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4. CONCLUSION 


We have defined and illustrated on a few perturbative examples a general method of 
calculating light-front hamiltonians which can be used for the relativistic description of 
interacting particles. The starting point in the calculation is a field theoretic expression 
for the bare hamiltonian density. This expression is multiply divergent in the physically 
interesting cases. Therefore, the hamiltonian theory requires renormalization. 

In the renormalization process, one calculates the whole family of effective hamiltonians 
as functions of the width parameter A which determines the range of the effective interactions 
on the energy scale. 

An effective hamiltonian of a small width A is much different from the initial bare hamil¬ 
tonian. It couples only those states whose masses differ by less than a prescribed amount. 
Thus, the effective theory contains only near-neighbour interactions on the energy scale. No 
scale is removed in the calculation but the correlations between dynamics at significantly 
different energy scales are integrated out. Therefore, in principle, the effective eigenvalue 
problem can be solved scale by scale using standard techniques for finite matrices which 
describe dynamics at a single scale. The remaining problems due to the small-a; massless 
particles need further investigation. 

Our formalism is based on the earlier work on renormalization of hamiltonians from 
Refs. and [Q where the hamiltonians are defined by their matrix elements in a given set 
of basis states. Wegner has developed similar equations for hamiltonian matrix elements in 
solid state physics. The present approach to renormalization of hamiltonians introduces 
the following features. 

Our equations are expressed in terms of creation and annihilation operators. Conse¬ 
quently, calculations of counterterms in perturbation theory can be performed without know¬ 
ing details of the specific Fock states which are needed to evaluate the matrix elements. This 
is useful because a large number of Fock states needs to be considered. The renormalization 
scheme is free from practical restrictions on the Fock space sectors. 

Expressing the effective hamiltonians in terms of the creation and annihilation operators 
and showing that the effective interactions are connected is a prerequisit to obtain the cluster 

m 


decomposition property. 


The effective interactions in our approach do not contain 


disconnected terms. The number of creation and annihilation operators in a single term is 
limited in perturbation theory by 2 + n{y — 2) where n is the order of a perturbation theory 
and V is the number of operators in the perturbing term. 

The physically motivated assumptions about the model space of effective states to be 
included in solving a particular problem are introduced after the effective hamiltonian is 
calculated. The interaction terms in the effective hamiltonian contain the similarity factors 
which diminish the dynamical significance of the Fock sectors with numbers of effective 
particles considerably different from the number of effective particles in the dominant sectors. 

The present operator formulation does not introduce spectator dependent interactions, 
even in the case where we include the sums of the invariant masses for incoming and outgoing 
particles in the similarity factors. The sums are useful for estimates of cutoff dependence in 
perturbation theory. 

The formalism explicitly preserves kinematical symmetries of the light-front frame. The 
structure of counterterms is constrained by these symmetries, including boost invariance. 
Hence, the number of possible terms is greately limited. Preserving boost invariance is 
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particulary important because it is expected to help in understanding the parton model and 
constituent quark model in QCD, simultaneously. 

It is essential to include the running of the coupling constants in the calculation of 
the small width dynamics. The examples of second order calculations we described in this 
article do not include the running coupling constant effects. Inclusion of these effects requires 
higher order calculations. Examples of related issues are following. The asymptotic freedom 
of QCD suggests that the effective coupling constant grows at small width. This rise can also 
enhance the initially small perturbative terms in the effective hamiltonians, such as masses 
of gluons or light quarks. In the meson-barion dynamics, the large couplings and narrow 
form factors are commonly used but the correlation between the size of the couplings and the 
form factor widths is out of the theoretical control. To gain some control on this correlation 
one needs to include the coupling constant dependence on the width. In QED, we expect 
the coupling constant to grow with the width. There, the hamiltonian calculation of the 
running coupling effects may help in understanding the ultraviolet structure of the theory 
and asymptotic nature of the perturbative expansion. 

Wegner’s equation can be adapted to building an operator approach similar to what we 
described in the present article. The initial equation which replaces our Eq. (2.29) when 
one uses the Wegner generator of the similarity transformation, is 

^ ^ 1|Wia,W2a],Wa] . (4.1) 

However, there is little flexibility left in the equation so that the widenning of the hamiltonian 
band is not readily available. 

Wilson and the present author have written a class of generalized equations for the flow 
of hamiltonian matrix elements and studied them in a simple numerical model. [?] These 
equations allow widening of the effective hamiltonian matrix at large energies. The gener¬ 
alized equations can also be adapted for the construction of the creation and annihilation 
operator calculus. Namely, 

^ = [F{W2a}, WaI . (4.2) 

These equations require detailed dehnitions of the similarity factors generated by the oper¬ 
ation F whose description goes beyond the scope of this article. 

In summary, the present formalism for renormalization of hamiltonians in the light-front 
Fock space provides a tool for working on a host of theoretical issues in particle dynamics. 
However, it remains to be verihed if the formalism can lead to quantitative improvements 
in our description of particles. The key problems are the lack of understanding of rotational 
symmetry beyond the nonrelativistic limit and the appearance of infrared singularities in 
gauge theories. Most urgent are the calculations of effective hamiltonians in the third and 
fourth order perturbation theory. 
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